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PREFACE 


This  report  is  the  first  draft  of  what  we  hope  will  eventually  be 
a  comprehensive  treatise  on  the  theory  and  calculation  of  EMP  coupling  to 
systems  located  in  the  source  region. 

There  lias  been  a  strong  tendency  for  many  years  to  rely  on  computer 
codes  for  EMP  coupling  calculations.  We  have  noticed  that  computer  codes 
built  in  the  absence  of  theoretical  understanding  almost  always  give  the 
wrong  answer  for  the  right  problem,  even  though  they  may  give  the  right 
answer  f o >'  the  Wrong  problem.  Source-region  coupling,  being  only  a  little 
more  difficult  subject  than  EMP  env 'ronments,  is  quite  amenable  to 
theoretical  analysis,  and  the  present  report  shows  how  such  analysis  can  be 
carried  out  for  some  important  examples. 

IVe  hope  to  add  to  this  report  over  the  next  few  years.  More 
examples  are  needed.  A  pioblem  rot  discussed  in  the  present  report  is  the 
effect  of  breakdown  in  air  rc.g.,  nuclear  lightning)  and  in  the  soil  on 
coupled  currents.  .'here  are  reasonable  prospects  that  sufficient  progress 
will  be  made  on  these  problems  in  the  next  year  or  so  that  they  can  be 
included.  Further,  experience  with  real  systems  such  as  MX  and  EoADS  may  show 
us  other  problems  that  need  analysis.  IVc  therefore  hope  that  the  copies  of 
tins  report  will  not  bo  bound  so  tightly  that  they  cannot  be  supplemented  by 
revisions  and  further  chapters. 
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CHAPTER  1 

INTRODUCTION  AND  BASIC  EQUATIONS 


1.1  INTRODUCTION 

A  nuclear  explosion  on  or  near  the  air-ground  interface  produces 
a  large  electromagnetic  pulse  (EMP) .  The  principal  source  of  the  EMP  is 
the  current  of  Compton  re-  oil  electrons  resulting  from  collisions  of  gamma 
rays  with  the  electrons  in  air  molecules.  The  Compton  current  is  signifi¬ 
cant  out  to  distances  of  several  kilometers  from  megaton  explosions.  Within 
this  source  region ,  the  air  conductivity,  associated  with  secondary 
ionization  produced  by  the  Compton  electrons,  has  a  strong  influence  on  the 
fields  generated.  The  presence  of  a  conducting  ground  also  has  a  strong 
influence . 


Calculations  of  the  coupling  of  electromagnetic  energy  into  systems 
located  within  the  source  region  must  take  into  account  the  existence  of  the 
gamma  rays,  the  Compton  current,  and  the  air  conductivity,  us  well  as  the 
EMP  fields.  Thus  source  region  coupling  is  more  coinp 1 icated  than  free-f i ell 
coupling ,  where  only  the  fields  need  to  be  considered.  Nevertheless,  a 
useful  approximate  theory  of  source  region  coupling  can  be  constructed,  and 
this  report  presents  the  theory  for  coupling  to  some  simple  but  practically 
relevant  system  geometries.  The  theory  will  hopefully  be  extended  to  other 
geometries  as  needs  arise. 

The  coupling  theory  presented  here  closely  parallels  the  theory 
of  source-region  EMP  environments  developed  previously  by  this  author 
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(References  1-2  to  1-4).  The  latter  theory  was  imoortant  in  that  it: 

•  gave  the  first  predictions  of  general  EMP  environments; 

•  showed  what  parameters  EMP  depends  on; 

•  showed  how  to  build  competent  computer  codes  for  more  detailed 
predictions; 

•  provided  accuracy  tests  for  the  codes. 

The  goals  and  uses  of  the  theory  of  source-region  coupling  are  similar: 

•  to  make  approximate  predictions  of  coupled  currents  and 
voltages,  especially  in  regimes  where  present  computer 
codes  are  not  valid; 

•  to  test  computer  codes  and  show  how  to  improve  them; 

•  to  provide  understanding  of  coupling  and  how  it  depends  on 
parameters ; 

•  to  allow  other  scientists  to  judge  the  correctness  of 
coupling  predictions. 

It  is  thus  hoped  that  this  report  will  be  useful  to  a  variety  of 
readers,  from  engineers  faced  with  the  task  of  making  predictions  for  actual 
systems  to  scientists  who  need  or  wish  to  judge  the  adequacy  of  our  under¬ 
standing  of  the  phenomena  and  of  methods  for  making  predictions.  In  this 
connection,  a  particular  reader  may  be  more  interested  in  some  sections  of 
this  report  and  less  interested  in  others.  We  have  tried,  however,  to  make 
all  of  the  report  readable  for  the  entire  spectrum  of  likely  readers. 

1.2  MAXWELL'S  EQUATIONS 

The  material  media  that  we  shall  be  dealing  with  most  commonly,  air 
and  soil,  are  essentially  non-magnet ic ;  that  is,  the  magnetic  permeability 
has  the  value  \i  appropriate  to  free  space.  Both  media  arc  generally 
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conductive,  and  the  soil  has  a  dielectric  permittivity  substantially  dif¬ 
ferent  from  Cq ,  the  free  space  value. 


The  two  time-dependent  Maxwell  equations  are  then 


38 

3t 


=  -  V  x  E  , 


(1-1) 


e 


3E 

at 


* 

j  + 


V  x  |  , 


(1-2) 


-v  2 

where  B  is  the  magnetic  field  (webers/m  ) ,  E  is  the  electric  field 
(volts/m),  and  J  is  the  current  density  (amps/m^) .  It  is  clear  that  these 
equations  are  sufficient  to  carry  the  fields  forward  in  time  if  initial 
values  are  given  and  if  J  is  specified.  From  these  initial  values  we  can 
evaluate  the  right-hand  sides  of  Equations  1-1  and  1-2,  which  then  tell  us 
how  E  and  B  will  change  in  the  next  infinitesimal  time  interval  fit. 

From  the  new  values  of  E  and  B  we  can  re-evaluate  the  right-hand  sides 
and  advance  the  fields  another  fit,  and  so  on.  This,  in  fact,  is  precisely 
how  numerical  solutions  of  Maxwell's  equations  are  obtained  (the  spatial 
derivatives  in  the  curl  operations  are  also  evaluated  in  finite  difference 
form) . 


Note  that  the  relation  of  cause  and  effect  in  this  way  of  looking 

at  Maxwell's  equations  is  different  from  what  most  of  us  were  taught, 

particularly  for  Equation  1-1.  The  picture  just  given  is  that  the  value  of 

V  x  E  determines  how  B  will  change  in  the  next  infinitesimal  time 

-> 

interval,  whereas  the  traditional  picture  is  that  a  changing  B  generates 

(inductively)  a  solenoidal  E ,  i.e.,  an  E  with  finite  curl.  Either 

picture  gives  the  same  result,  namely  that  the  right-  and  left-hand  sides 

are  equal,  and  we  do  not  actually  need  to  decide  which  side  causes  the 

other.  However,  the  new  picture,  which  is  the  one  generally  used  by  physicists, 

makes  it  easier  to  understand  how  time-dependent  solutions  evolve. 
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There  are  cwo  other  Maxwell  equations, 

V  •  B  =  0  ,  (1-3) 

V  •  (cE)  =  p  ,  (1-4) 

where  p  is  the  charge  density.  It  would  appear  from  the  foregoing  discus¬ 
sion  that  these  equations  are  not  needed  in  advancing  the  fields  in  time. 

Such  is  indeed  the  case,  for  taking  the  divergence  of  Equation  1-1  gives 

(V4)  =  -  V  •  (VxE)  =  0  .  (1-5) 

(The  divergence  of  the  curl  of  any  field  vanishes.)  This  equation  says 
that,  if  Equation  1-1  is  satisfied,  V  •  B  will  be  independent  of  time  at 
all  points  in  space.  Thus  if  V  •  B  vanishes  everywhere  initially,  then 
the  solution  of  Equation  1-1  will  have  V  •  1?  =  0  everywhere  at  all  times. 
Therefore  Equation  1-3  needs  only  to  be  imposed  as  a  condition  on  the  initial 
magnetic  field.  If  the  initial  magnetic  field  vanishes,  liquation  1-3  is 
satisfied . 


To  understand  the  role  of  liquation  1-4,  take  the  divergence  of 
liquation  1-2,  and  obtain 

gy  (V-rli)  =  -  V  •  J  .  (1-0) 

Now  the  conservation  of  charge,  which  is  a  well  verified  law  of  nature, 
states  that 


<}P 

Dt 


V  • 


■y 

.!  . 


Subtracting  liquation  1-7  from  liquation  l-(>  gives 

(V-cii-p)  =  0  . 


(1-7) 


(1-8) 


Thus  it  follows  from  Equation  1-2  that,  if  the  quantity  in  parentheses 
vanishes  everywhere  initially,  it  will  vanish  everywhere  at  all  times. 
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Therefore  Equation  1-4  also  needs  only  to  be  imposed  as  a  condition  on  the 
f  ,  ,  ->• 

initial  E  and  p.  If  E  and  p  both  vanish  initially,  Equation  1-4 
is  satisfied. 

If  B,  E  and  p  all  vanish  initially,  we  need  only  concern  ourselves 
with  Equations  1-1  and  1-2.  Note  that  these  equations  do  not  contain  p 
at  all;  p  need  not  be  calculated.  If  p  is  desired,  it  can  be  found  by 
time  integration  of  Equation  1-7. 

Note  that  Maxwell's  equations  and  the  conservation  of  charge, 

-y.  ..y  —y 

Equation  1-7,  arc  linear  in  the  variables  J,  p,  E  and  B.  Thus  if  current 
density  .J  produces  p. ,  p.  and  B.  ,  and  J,  produces  p9,  E7,  and  B~, 
then  current  density  ^  will  produce  p^  +  p2 ,  E^  +  and  +  B^ . 

We  have  assumed  here  that  e  (and  p„  of  course)  is  the  same  it  all  cases. 

,  ^  *>  -y 

This  linearity  is  somewhat  restricted  in  practice  when  .J  depends  on  E, 
as  we  shall  see. 

1.3  SOURCE  AND  CONDUCTION  CURRENTS 

In  EME  problems  the  current  density  is  made  up  of  two  parts.  First, 
there  is  the  source  current  of  Compton  recoil  electrons  produced  by 

the  flux  of  gamma  rays,  which  is  the  source  of  the  EMI’.  Second,  there  is 
the  conduction  current  ,J^  associated  with  the  flow  of  low-energy  electrons 
and  ions  induced  by  the  electric  field.  The  total  current  is  the  sum 

J  =  J  +  J  .  (1-0) 

sc 

The  Compton  current  is  formed  by  recoil  electrons  that  have  start¬ 
ing  energies  of  the  order  of  1  McV.  These  electrons  arc  stopped,  in 
material  media,  by  inelastic  collisions  with  the  media  atoms.  In  air,  the 
stopping  range  of  the  recoil  electrons  is  a  few  meters.  Therefore,  if  the 
electric  field  E  is  less  than  about  10‘  V/m,  the  effect  of  this  field  on 
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range  will  be  small  and  may  be  neglected.  If  the  field  were  10^  V/m,  the 
range  would  be  substantially  affected  by  the  field.  The  magnetic  field  B 
deflects  the  recoil  electrons.  The  deflection  will  be  small  if  the  Larmor 
radius  is  long  compared  with  the  stopping  range.  The  Larmor  radius  is  a  few 
meters  when  B  is  about  20  gauss  =  2  x  10  webers/m  .  Thus  for  magnetic 
fields  of  this  size  or  larger  the  deflection  will  be  substantial. 

In  many  applications  the  fields  are  less  than  the  critical  values 

just  given.  In  these  cases  we  may  assume  that  depends  only  on  the  gamma 

flux  and  is  independent  of  the  fields.  In  cases  where  the  fields  are  higher, 

we  shall  estimate  corrections  to  J  due  to  the  fields. 

s 

In  soil,  the  recoil  electron  range  is  only  a  few  millimeters  (soil 

3 

is  about  10  times  more  dense  than  air) .  Here  the  fields  are  never  large 
enough  to  affect  the  Compton  current.  Gamma  rays  are  attenuated  by  a  factor 
e  in  15  to  20  cm  of  soil.  Hence  the  Compton  current  is  appreciable  only  in 
the  top  meter  or  two  of  the  ground. 

The  conduction  current  is  generally  well  approximated,  in  both  air 
and  soil,  by  Ohm's  law, 

Jc  =  an  ,  (1-10) 

where  a (mhos/m)  is  the  conductivity.  In  air,  a  depends  somewhat  on  li , 
making  Maxwell's  equations  nonlinear.  Wc  can  usually  choose  an  li- independent 
value  of  o  which  over-estimates  coupling  effects.  Since  the  air  conductivity 
results  from  ionization  produced  by  the  Compton  recoil  electrons,  o  depends 
on  time  and  position.  In  the  ground,  a  is  independent  of  li,  except  at  very 
high  fields  where  breakdown  occurs.  It.  is  also  little  affected  by  ionization, 
except,  at  very  high  dose  rates.  Ft  may  be  assumed  independent  of  time  and 
position,  but  it  does  depend  on  the  frequency  of  the  driving  li-field  (as 
does  also  c)  .  Tlicse  points  will  be  discussed  in  detail  in  later  sections. 
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A  STANDARD  FORM  FOR  MAXWELL’S  EQUATIONS 


We  shall  write  the  dielectric  permittivity  in  terms  of  the  value 


er  for  free  space, 


£  =  £  £n  , 
r  0 


(1-11) 


where  e,  is  the  relative  permittivity.  If  we  also  make  use  of  Equations  1-9 
and  1-10,  the  Maxwell  Equation  1-2  becomes 


c)  B  I  -v 

e  en  —•  =  -  J  -  uE  +  —  V  x  B  . 
r  0  dt  s 


(1-12) 


It  is  convenient  to  replace  and  by  two  other  parameters,  namely 

the  speed  of  light  in  vacuum, 


c  =  — -----  w  3  x  IQ8  m/sec  , 


(1-13) 


and  the  impedance  of  Free  space, 


Z n  =  .  ---  a;  1 2 On  w  377  ohms 

0  Mc0 


(l-U) 


These  equations  can  be  solved  for  p^  and  r.y, 

l‘()  =  ^(/c  '  c()  =  l/cZo  • 


f  1  -IS) 


Inserting  these  expressions  in  liquation  12  gives 


t:  > 

rdf  -  y  >  > 

-  T,  =  -  Z_J  -  Z  oli  +  cV  x  B  . 
C  d  L  ()  s  0 


(1-16) 


livery  term  in  this  equation  now  has  the  dimensions  volt/m"';  note  that  cB 
is  the  electric  field  of  a  wave  in  vacuum,  as  follows  from  liquation  1-1. 


In  the  remainder  of  this  report  wje  shall  drop  the  subscript  r  on 
i  ;  i  will  always  mean  the  relative  permittivity.  Thus  our  standard  form 
for  Maxwell's  equations  is 
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We  note  here  also  the  relation  between  the  current  I (amps)  in  a  wire  and 
the  static  magnetic  field  Bq  encircling  it  at  radius  r, 

27rrBQ  =  u  I  =  ZQI/c  or  I  =  2TrrcB0/Zo  .  (1-3  9) 

Finally,  note  that 

Zg/2Tr  «  120'rr/2ir  =  60  ohms  .  (1-20) 
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CHAPTER  2 

THE  IMPEDANCE  OF  SOILS 


2.1  THE  RC  MODEL 

Scott2-1  measured  the  conductivity  and  permittivity  of  many 
samples  of  soil  over  the  frequency  range  10^  to  10  Hz.  He  noticed  that 
the  results  correlated  quite  well  with  the  water  content  of  the  soil.  He 
made  mathematical  fits  to  his  a  and  e  results  as  functions  of  frequency 
and  water  content.  In  making  these  fits,  he  made  no  attempt  to  ensure  that 
a(ai)  and  e  (w)  bear  the  relation  to  each  other  required  by  causality. 

Longmire  and  Longley2-2  noticed  that  Scott's  fits  could  be  refitted 
very  well  by  assuming  that  between  opposite  faces  of  the  soil  sample  there 
is  an  RC  network  of  the  type  shown  in  Figure  2-1.  In  this  network,  1 / R ^ 
represents  the  zero-frequency  conductivity,  C  represents  the  infinite. 
frequency  dielectric  constant,  and  the  other  branches  account  for  the  change 
in  o  and  e  with  frequency.  A  good  fit  was  obtained  with  one  such  branch 
for  each  decade  in  frequency  covered,  with  the  time  constant  IHC  of  the 
relevant  branch  chosen  equal  to  the  reciprocal  of  the  median  w  in  that 


Figure  2-1.  Network  representing  soil  impedance. 
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decade,  i. e.  ,  (R^CH  )  ]  =  Vio  2irt\,  where  f.  is  the  frequency  at  the  lower 
end  of  the  decade  in  question.  Thus  the  products  R^C^  were  chosen 
arbitrarily,  t.o  cover  the  frequency  range  uniformly.  The  ratios  Ih/C.  and 
were  then  chosen  to  fit  Scott's  e  curve.  Only  one  parameter,  R^, 
was  then  left  to  fit  the  a  curve,  but  it  was  found  that  a  good  fit  to  o 
was  obtained.  Furthermore ,  it  was  noticed  that  changing  the  fit  for  a 
soil  of  different  water  content  was  accomplished  by  scaling  all  of  the 
resistors,  except  ,  by  the  same  factor  and  leaving  the  capacitors  unchanged. 
Longmire  and  Smith2-3  used  these  results,  and  data  at  higher  frequencies, 
to  make  a  "universal  impedance  function"  of  soils  over  the  frequency  range 
102  to  108  Hz. 


2.2  THE  SOIL  ADMITTANCE 


The  Maxwell  Equation  1-18  for  fields  varying  as  e^llJt  (j2  =  ~  1) 


takes  the  form 


nii  =  -  z0^s  ♦  cv  x  ri  , 


(2-1) 


where 


r|  ZqO  +  j  (meters)  ^ 


(2-2) 


In  Equation  2-1  the  conduction  and  d isplaceinent  currents  have  been  combined 
into  tne  term  on  the  left.  If  we  define  the  admittance  Y (w)  of  unit 
volume  of  soil  by  the  relation  between  total  E -driven  current  5,,  and  E 


->  y 

•J,.  =  YE  , 


(2-3) 


then  obviously 


„  i) 

M  =  L  nY  ,  Y  =  ■=-  -  o 
0 


iwnr. 


0  • 


(2-4) 


The  dimensions  of  Y  arc  mhos/meter,  while  those  of  n  are  (meters) 
We  shall  call  rj  the  relative  admittance. 


IB 


The  admittance  of  the  RC  network  is 


I 


jioC 


^  R  +  •'  U>l"m  ^  i  *  ju)R  C 
U  n  •  n  n 


(2-5) 


The  real  and  imaginary  parts  of  Y  are  related  to  o  and  r.  by  liquation 


2-4.  Defining  RC  rates  (3  by 


Bn  2  <W 


-1 


(2-0) 


Reference  3  fits  Scott's  data  by  the  formula 

1 3  i  <or. 


i  =  o.  +  not  e  +  }  a 

0  ■  0  n  1  ^ 

n- 1 


0_ 

+  iui/(3 . 


(2-7) 


Here  is  the  zero-frequency  conductivity,  tv>)  is  the  infinite-frequency 

relative  permittivity,  r,  is  the  permittivity  of  free  space  in  MRS  units, 

the  b  are  a  fixed  set  of  rates, 
n 

13  =  2 a  (  10) n_  1  sec  '1  ,  (2-K) 

and  the  a  are  a  set  of  dimensionless  fit  coefficients,  for  soil  contain 
n 

ing  10  percent  water  by  volume  the  fit  parameters  are  given  in  Table  2-1. 


for  this  fit  the  relative  permittivity  and  conductivity  are 


13 
+  V 


ir-- 1  1  +  (.„/|<) 


13  a  i 

ii  0  n 

+  L  . . 


0 


(2-0) 


(2-10) 


n  -  1  1  +  { '•'/l'11 ) 


It  can  be  seen  that  >  decreases  while  11  increases  with  increasing  frequem 
Craphs  of  >  and  o-u  as  functions  of  frequency  arc  given  in  figures  2-2 
and  2-3.  figure  2-1  shows  how  o  varies  with  water  content,  and  also 
gives  the  scale  factor  f  bv  which  the  rates  !•  must  be  multiplied  for 
different  water  content. 
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from  Equation  2-9,  with  Equation  2-8  and  Table  2-1 .  Smooth  curve  is  drawn 
through  decade  points. 


Figure  ?-4.  Frequency  scale  factor  F  and  zero  frequency 
conductivity  n  from  Scott's  results. 


Conductivity  o0 (mho/meter) 


Table  2-1.  Fit  parameters  for  soil  containing  10  percent 
water  by  volume. 


ao 

=  8x10"' 

*  mho/m 

e 

oo 

=  5  (relative) 

n 

an 

n 

an 

n 

an 

1 

3.40(6) 

6 

1.33(2) 

11 

9.80(-l) 

2 

2.74(5) 

7 

2.72(1 ) 

12 

3 . 92 ( - 1 ) 

3 

2.58(4) 

8 

1.25(1) 

13 

1  73 ( - 1 ) 

4 

3.38(3) 

9 

4.80(0) 

5 

5.26(2) 

10 

2.17(0) 

This  fit  is  expected  to  be  good  for  frequencies  between  10^  and  10^  H 
for  a  wide  range  of  water  contents.  The  author  has  never  seen  any  data  that 
cannot  be  fitted  reasonably  well  by  this  model  by  adjusting  only  the  assumed 
water  content  and  the  value  of  aQ  (to  a  value  that  may  be  different  from 
that  indicated  by  Figure  2-4) . 


The  fit  for  the  relative  admittance  is,  according  to  Equations  2-4 

and  2-7, 


n 


zoao 


13 


£ 

n=l 


jh)/c 

n  1  +  ju)/3 

n 


(2-11) 


2.3  CAUSALITY  AND  REALITY 


The  requirement  of  causality  is  that  the  current  must  vanish  until 
a  field  is  applied.  For  example,  let  E(t)  be 


E(t)  =  0  ,  t  <  0  , 

-Yt 


(2-12) 


E(t)  -  Enc 


t  £0  , 


( Y  =  real ,  positive) 


The  Fourier  transform  of  I:  is 


i 

/ 


H(uJ  =  /  K(t)e  •'  dt  = 


j  cot 


0 

Y  +  ju) 


(2-  n  j 


The  current  density  in  the  frequency  domain  is  then 

0, 


J(w)  =  Y  (to) 


_ j0_ 

Y  +  jw 


and  in  the  time  domain 


J(t) 


27  / V(»> 


j(0t  . 

C  Uii!  . 


(2-14.) 


(2-15) 


For  t  <  0,  the  integration  contour  can  be  extended  to  enclose  the  negative 
imaginary  half  plane.  The  factor  l/(Y+.i<‘>)  has  a  pole  at  «  =  jY.  If 
Y (u>)  lias  no  poles  in  the  negative  imaginary  half  plane,  the  integral  will 
vanish  as  required,  for  t  <  0.  Inspection  of  liquation  2-7  shows  that  the 
poles  of  Y  arc  at  id  =  j 3^ .  in  the  positive  imaginary  half  plane.  Thus 
causality  is  satisfied.  The  generally  required  relation  between  a  physical 
o  and  r  inn  is  that  they  must  form  the  real  and  imaginary  parts  of  a 
complex  function  which,  when  analytically  continued  from  the  real  id  axis 
into  the  negative  imaginary  half  plane,  lias  no  poles  there.  Any  RC  network 
provides  this  property.  Resonances  in  *’  and  >.  could  be  accommodated  by 
adding  inductances,  but  it  appears  that  none  are  needed. 


The  fa(  i  that  the  electric  field  H(t)  and  the  current  density 
J(t)  are  real  functions  of  time  places  another  condition  on  the  admittance 
Ypn)  and  the  relative  admittance  n(o).  lor  general  real  lift)  ,  liquation  2-15 
shows  that  the  complex  conjugate  H*('n)  is  related  to  H(w)  by 

l;.*(o)  =  li  (-<■•)  .  (2- lb) 

The  same  relation  holds  between  .)*(■.  )  and  .Jpu).  Since 

Y  (  a )  =  ■)  (id)  /H  (>o)  ,  (2-17) 
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it  follows  that  Y  and  n  also  obey  the  real i ty  condition 


Y*  (in)  =  YC-w)  , 
n*(uo;  =  n(-fa)  . 

From  liquations  2-2  and  2-4,  it  then  folloi  that 
o  ( -to)  r  o(uj)  ,  e(-w)  =  e((jj)  . 

2.4  EXPONENTIALLY  RISING  FIELD 


(2-18) 


(2-19) 


In  the  early  part  of  the  I.MP  the  electric  field  i  isos  approximately 
exponentially  to  a  level  fairly  near  the  peak  field, 

h (t)  w  F.0cat  .  (2-20) 


Hence  it  is  useful  to  evaluate  n  for  the  case  in  v.’hich  jio  is  replaced 
by  ' x .  One  obtains  the  real  expression, 


■1 


zoao 


1 3 


L 

n=  1 


_ a/c _ 

dn  1  +  a/F 
n 


(2-21) 


A  graph  of  t;  as  a  function  of  u  for  the  10  percent  water  soil 
is  shown  in  Figure  2-5.  The  same  figure  shows  the  real  and  imaginary  parts 
of  1 1  as  a  function  of  <o  for  the  oscillatory  case.  Note  that  while 
there  is  no  simple  relation  between  n(u)  and  n^po)  and  t|^((o),  except 
that  contained  in  liquation  2-11,  q(o.)  is  not  far  from  the  sum  of  n  (w) 
and  q.  (oj)  for  w  =  a.  (Ac'ually,  n(u)  is  a  little  less  than  the  sum.) 

We  shall  call  the  case  graphed  in  Figure  2-5  our  standard  soil. 
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TIME  DOMAIN  TREATMENT 


Reference  2  showed  how  to  treat  frequency-dependent  media  in  the 
time  domain.  Write  the  total  E-driven  current  that  flows  into  the  network 
as  the  sum  of  the  currents  in  its  branches, 


where 


t r  -  c  if  +  i  +  V  3  , 

E  °°  3t  R_  n  ' 

0  n 


J  *  (E  -  —  f  S  dt)/R 
n  LnJ  n  n 

-oo 

t. 

■  B  (C  ft  -  /  J  dt)  . 
nv  n  J  n 


(2-22) 


(2-23) 


(2-24) 


Comparison  of  Equation  2-5  with  the  fit  formula  (2-7)  establishes  the  relations 

a 

r 

7 

\) 

Thus  liquations  2-22  and  2-23  can  be  written  as 

c 


Cco  =  ,:wEo  *  R0  “  °0  »  Cn  anc0  3  cZ,' 


•>  C»  ^  y 

Z0',li  =  "ir  3t  +  V’o1'  +  Z0  53  ,J„ 

n 

J  =  3  (  3-1.  it  _  f, j  dt)  . 
n  n\cZ„  J  n  / 


(2-25) 


(2 -2b) 


Inserting  Z  i,  for  pit  in  Equation  2-1),  we  can  take  that  equation  back 
to  the  time  domain,  with  the  result. 


(f j  r) I •  - y  ■  «  — v 

i-  st-  -  VWLV  *  cVx  15  ■ 


(2-27) 


This  is  the  time-domain  form  of  the  Maxwell  equation  (1-18)  for  t lie  frequency- 
dependent  medium.  The  are  to  be  obtained  from  Equation  2-2(i,  which  can 
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be  converted  into  a  differential  equation,  if  desired,  by  differentiating 
it  with  respect  to  t. 


The  fact  that  J  is  a  vector  means  that  Equation  2-26  must  be 
solved  for  each  non-vanishing  component  of  E.  In  a  stratified  medium,  the 
parameters  u  ,  em,  and  the  u  could  have  different  values  for  different- 
directions  . 
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CHAPTER  3 

SOURCE  REGION  ENVIRONMENTS 


3.1  INTRODUCTION 

A  complete  discussion  of  liMP  environments  is  beyond  the  scope  of 
this  report.  Theoretical  discussions  are  given  in  References  1-2  to  1-4,  and 
many  detailed  computer-based  calculations  have  been  made.  Access  to  much  of 
the  available  information  is  controlled  by  such  U.S.  Government  agencies  as 
the  Defense  Nuclear  Agency  and  the  Air  Force  Weapons  Laboratory.  The 
Goverivnent  normally  provides  HMD  environment  specifications  for  systoms  it 
sponsors . 


For  the  purposes  of  this  report,  it  will  be  necessary  to  know  only 
the  gone -a  1  features  of  the  liMP  environment,  such  as  the  order  of  magnitude 
of  rise  times,  amplitudes,  and  decay  times,  and  approximate  relations  between 
Compton  current,  air  conductivity,  and  the  fields.  These  will  lie  presented 
in  this  chapter.  Predictions  of  coupling  to  actual  systems  should  use 
environments  supplied  by  the  sponsoring  agency. 

While  it  should  not  be  assumed  that  the  environments  hypothesized 
in  this  chapter  are  precisely  correct  for  any  particular  real  case,  they  are 
representative,  in  crude  approximation,  of  those  that  might  be  observed  at  a 
point  on  or  near  the  ground  surface  at  a  distance  of  500  meters  from  a  1 
megaton  explosion  at  the  ground  surface  or  within  a  few  tens  of  meters  above 
the  surface. 
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3.2 


THE  GAMMA  FLUX 


Figure  3-1  shows  a  gamma  flux  as  a  function  of  time,  which  will  be 
used  as  an  example  in  later  sections  of  this  report.  The  quantity  graphed 
is  actually  the  dose  rate  delivered  to  air,  and  the  relation  between  dose 
rate  and  the  actual  flux  F  of  gamma  energy  is 


n  fV~MeV 

Y  L  2 
in  sec 


Y 


■)  «  2  x  1013  D  (£^-) 
sec 


(3-1) 


The  average  energy  of  the  gammas  is  about  2  MeV  per  photon,  although  the 
total  spectrum  covers  the  range  from  a  fraction  of  1  MeV  to  many  MeV. 


The  time  indicated  in  Figure  3-1  is  actually  retarded  time,  i.o., 
the  time  origin  is  set  when  gammas  first  begin  to  arrive  at  the  observer's 
position.  The  time  for  gammas  to  travel  500  meters  is  1.67  microseconds, 
since  the  speed  of  light  is  300  m/ps . 


In  a  crude  approximation,  the  gamma  flux  is  collimated  in  the 
radial  direction  from  the  burst  point  (the  point  of  the  nuclear  explosion) . 
In  the  case  assumed  here,  the  flux  is  approximately  horizontal.  Due  to 
scattering  and  finite  source  size,  the  actual  angular  distribution  of  the 
gammas  covers  several  tens  of  degrees  around  the  radial  direction. 


Figure  3-1  indicates  what  sources  arc  responsible  for  various 
parts  of  the  gamma  flux.  The  prompt  gammas  arc  emitted  by  the  nuclear 
device  itself.  Air  inelastic  gammas  are  made  in  inelastic  collisions  of 
energetic  neutrons,  emitted  by  the  device,  with  the  nuclei  of  air  atoms. 
Ground  capture  gammas  arc  produced  when  neutrons  are  captured  in  the  ground 
near  the  burst  point.  Air  capture  gammas  arc  produced  later  when  neutrons 
arc  captured  in  the  air.  Fission  product  gammas  arc  emitted  over  long  times 
by  the  nuclear  fragments  resulting  from  fission  of  uranium  or  plutonium. 

Of  these  sources,  only  the  prompt  gammas  have  effectively  a  point  source; 
the  others  originate  in  volumes  with  dimensions  of  the  order  of  a  few 
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hundred  meters.  Figure  3-1  has  been  drawn  as  if  each  of  the  sources,  except 
the  prompt  and  fission  product  gammas,  were  a  decaying  exponential.  Note 
that  the  doses  from  each  of  the  sources  arc  about  equal,  since  the  sources 
with  lower  dose  rates  last  longer  in  time. 

The  rise  of  the  prompt  gamma  flux  has  been  chosen  exponential 
in  time,  i . e . , 

I)  =  Aeat  ,  (3-2) 

where  A  is  a  constant  and  a  has  been  taken  as 

a  =  2  x  1()8  sec"1  .  13-3) 

The  exponent  Lai  form  is  crudely  representative,  and  is  convenient  for  cal¬ 
culations.  The  value  of  a  chosen  is  in  t ho  correct  range,  but  should  not 
be  taken  as  either  an  upper  or  a  lower  bound  on  actual  values. 

3.3  THE  COMPTON  CURRENT 


(lamina  rays  traveling  through  matter  collide  occasionally  with  atomic 
electrons,  knocking  the  electron  generally  forward  and  scattering  the  photon. 
The  mean  free  path  of  the  gammas  in  air  for  those  Compton  scatter lnjg  collisions 
is  a  few  hundred  meters.  The  recoil  electrons,  which  have  initial  energies 
of  about  1  MeV,  move  forward  an  average  of  a  few  meters  before  stopping,. 

Thus  a  steady  flux  of  gammas  will  produce  a  steady  flux  F  of  recoil 
electrons,  in  the  same  direction,  of  about  1  percent'  of  the  gamma  flux.  The 
relation 


:  w  )(,i:  /l-MoV 

c  \  m'"sec  '  \m“sec 


> 


(3 -'ll 


holds  approximately  in  air  (and  in  other  media  of  low  atomic  number,  such  as 
soil)  over  the  gamma  energy  range  of  interest  hero.  From  liquations  3-4  and 
3-1,  and  the  charge  of  an  electron,  one  can  deduce  an  approximate  relation 
for  the  Compton  current  density  .J^  (the  source  of  1-MI’), 
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The  lifetime  of  a  recoil  electron  (before  stopping)  is  ubout  10  '  second 

-9 

in  real  time,  and  about  10  second  in  retarded  time  (the  electron  moves 
forward  at  ubout  0.9  of  the  speed  of  light).  Thus  the  steady-state  relation 
(3-5)  is  valid  when  the  changes  in  D  in  periods  of  10  second  are  small 
compared  with  I).  This  condition  is  fairly  well  satisfied  by  the  dose  rate 
in  figure  3-1.  The  Compton  current  density  graphed  in  figure  3-2  is  obtained 
from  liquation  3-5. 

3.4  THE  AIK  CONDUCTIVITY 


f.ach  Compton  recoil  electron,  in  slowing  down,  produces  about 
3  x  10  pairs  of  secondary  electrons  and  positive  ions,  which  make  the  air 
electrically  conducting.  The  rate  of  production  of  ionization  is  directly 
proportional  to  the  dose  rate, 


s(.L°»  1U  2  x 

m  “'sec  ' 


(3-0) 


The  free  electrons,  because  of  their  small  mass,  respond  more  quickly  than 
ions  to  applied  electric  fields,  and  are  the  dominant  contributors  to  the 
air  conductivity  at  early  times.  However,  the  electrons  gradually  attach 
themselves  to  (),  molecules,  forming  the  negative  ion  0,.  The  rate  a 
of  attachment  per  electron  Is  about 

a  f-  1  x  lo^sec  '  ,  (3-7) 

in  sea-level  air.  Thus  the  density  N  of  free  electrons  satisfies  the 
equal-  i  on 

;)N  ^ 

.y:°  =  S  -  aN  .  1  (3-8) 

<)t  o 
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If  the  dose  rate  rises  exponentially,  as  in  liquation  3-2,  the  solution  is 

electrons 

3 
m 

On  the  other  hand,  if  the  dose  rate  changes  little  in  periods  of  10 
second , 

N  *t  -  *  2  x  107  D  .  (3-10) 

e  a 

This  equation  is  obtained  from  (3-9)  by  setting  a  =  0. 


)- 


2  x  10 


15  D 


a  +  a 


(3-9) 


In  the  presence  of  an  electric  field  H,  the  free  electrons  drift 
through  the  air  at  an  average  speed  v  which  is  roughly  proportional  to  Ii, 

v  =  Meli  •  (3-11) 

The  electron  mobi  1  ity  |J  is  of  the  order  of  magnitude 


U  w  0.3 
c 


_m  .volt: 
sec7  m 


in  sea-level  air.  Actually,  p 
.Joule  heating  of  the  electrons, 
3  x  io3  <  |i  <  3  x  K)5,  is 


(3-12) 

depends  significantly  on  li ,  because  of 
and  a  better  expression,  over  the  range 


P 


e 


li*  =  3  x 


104  V/m  . 


(3-13) 


The  electrical  conductivity  a  is  the  ratio  of  the  conduction 
current  density  .J  to  Ii.  from  the  equations  above  and  the  electron 
charge  o,  the  electronic  conductivity  can  be  deduced  as 

o  =  N  e|i  0.8  x  j ()  ^  — -  mho/m  .  (3-14) 

n  e  e  a  +  ay  I-. 

After  the  peak  of  the  dose  rate,  a  should  be  set  equal  to  zero  here. 

At  late  times,  the  air  conductivity  is  dominated  by  positive  and 
negative  ions,  because  they  disappear  more  slowly  than  electrons.  The 
equation  governing  the  positive  ion  density  N+  is 
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bN  N 
+ 


Here  N  us  the  density  of  negative  ions  (0o  etc.),  and  b  is  the  mutual 
neutralization  coefficient  for  positive  and  negative  ions.  The  value  of  b 
is 


10 


12 


m'Vscc 


(3-10) 


At  late  times,  most  of  the  electrons  that  have  been  produced  are  attached 
to  0-,,  so  that  N  w  N  .  The  solution  of  liquation  3-15  is  such  t.ha<  the 
two  terms  on  the  right-hand  side  nearly  balance,  and  3N  /dt  is  sma  i 
compared  with  either  of  these  terms,  so  that 


N  re  N 

+ 


(3-17) 


The  mobility  of  the  ions  is  about 


P 


i 


_in _ /V 

sec' m 


The  ion  conductivity  can  then  be  deduced  as 

o.  j-  2N+e|i.  w  2 . 5  '  10  ^  \[i)  mho/m  . 
Comparison  of  this  result  with  liquation  3-14  with  a  -•  0, 

.  7 

that  o.  and  u  are  equal  when  I)  u  10  rads/sec. 


(3-18) 


(3-19) 
H*,  shows 


The  electronic  and  ionic  conductivities  and  the  total  conductivity 
are  graphed  in  figure  3-3.  for  the  electronic  conductivity,  li  has  been  set 
equal  to  li*,  so  that  the  result  is  indicative  rather  than  precise, 
especially  at  early  times  when  li  will  be  larger  than  li*. 
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Time  (seconds) 

Figure  3-3.  Air  conductivity  for  the  dose  rate  of  Figure  3 
a  and  or.  are  shown  separately,  along  with  the 


3.5  SCALING  WITH  DISTANCE! 

The  gamma  flux  of  Figure  3-1  is  crudely  representative  for  an 
observer  at  500  meters  from  a  1  megaton  burst  near  the  ground  surface.  For 
other  yields  „>j  and  distances  r,  the  gamma  flux  scales  roughly  as 

-r/A 

D  ~  yy  .  (3-20) 

r 

Here  A  is  the  effective  mean  free  path  of  the  gammas  in  air;  a  representative 
value  in  sea- level  air  is 

A  w  300  meters  .  (3-21) 

According  to  Section  3.3,  the  Compton  current  density  has  the  same 

scaling  as  D.  The  scaling  of  the  air  conductivity  is  less  simple; 

scales  as  D,  while  a.  scales  as  \/!T.  At  early  times,  a  is  dominant  at 

X  c 

most  distances  of  interest  in  this  report. 

3.6  THE  RADIAL  E  FOR  SPHERICAL  SYMMETRY 

The  Compton  current  J  and  the  air  conductivity  o  are  approxi¬ 
mately  independent  of  angle  about  the  burst  point.  The  presence  of  the 
ground  destroys  complete  spherical  symmetry,  of  course.  However,  for 
observers  above  the  ground,  e.g.,  at  power  line  heights,  the  effect  of  the 
ground  on  the  fields  will  not  occur  immediately,  but  will  be  delayed  by  the 
finite  speed  of  light  and,  more  importantly,  by  the  diffusion  time  of  fields 
through  the  conducting  air  between  ground  and  elevated  observer.  It  is 

therefore  useful  to  examine  the  solution  of  Maxwell's  Equations  1-17  and  1-18 

- y 

for  the  case  of  radial  J  and  spherical  symmetry  in  J  and  a.  Note 
that  c  =  1  in  the  air. 

The  fields  f!  and  6  start  from  zero.  Integrating  liquation  1-18 
over  a  small  time  interval  will  give,  by  integration  of  ,  a  radial  and 
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spherically  symmetric  E.  The  curl  of  such  an  E  vanishes,  so  that 
integration  of  Equation  1-17  leaves  B  =  0.  Thus  the  V  x  $  term  in 
Equation  1-18  remains  equal  to  zero,  and  E  remains  radial  and  spherically 
symmetrical.  The  vector  signs  may  be  dropped  and  Equation  1-18  becomes 


-  ||  =  -  ZnJ  -  ZnoE 
c  3t  Os  0 


(3-22) 


Note  that  this  equation  contains  no  spatial  derivatives:  E  is  determined 
at  each  point  by  the  local  J  and  a. 

At  sufficiently  early  times  J  and  a  are  small  and  E  will  be 
small,  and  oE  will  be  negligible  compared  with  J..  In  this  time  frame 


E  «  cZ, 


(3-23) 


If  Jg  rises  as  exp(at),  E  will  also,  and 


0  » - -  J  . 

a  s 


(3-24) 


In  this  time  frame,  it  can  be  said  that  is  charging  up  the  capacitance 

of  space  (cZ„  =  l/en) . 


Eventually  al!  will  become  comparable  with  J  ,  if  the  dose  rate 
is  large  enough.  In  this  case  the  91i/3t  term  can  be  neglected  in  Equation 
3-22,  giving  the  approximate  solution 


n  =  -  r  E  . 
a  s 


(.3-25 


This  equation  defines  the  saturated  field  1:  ,  which  is  such  that  conduction 
current  cancels  Compton  current.  Since  .J  and  a  are  both  proportional 
to  b  (at  early  times),  they  tend  to  rise  and  fall  together,  so  that  E  is 
almost  constant  after  saturation  is  reached.  Thus  DE/Bt  is  indeed  small. 
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The  value  of  1-  can  be  determined  from  liquations  3-5,  3-14,  and 
3-25.  Th°  result  comes  out  directly  as 


li  w  2.5  x  10 
s 


t 


+_a 
'  8 


10 


s 

rf* 


With  the  value  of  li*  given  in  liquation  3-13,  this  result  becomes 

o 


S 
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(3-2<>) 


Thus,  if  saturation  occurs  during  the  exponential  rise  in  the  sample  gamma 

O 

flux  with  a  -  2  *  10  ,  the  peak  li  will  be 


1.8  x  K)  V/m  . 


(3-27) 


However,  after  the  peak  in  the  gamma  flux  <x  -  1),  and 
li.  k.  2  ■<  l O’1  V/m  . 


(3  ~2,V) 


At  late  times,  when  the  ion  conductivity  is  dominant,  li .  falls  rough lv  as 

Vft.  However,  by  this  time  the  effect  of  the  ground  asymmetry  will  usually 

be  Felt,  figure  3-4  shows  lift)  for  the  sample  case,  neglecting  the  ground 

effect  altogether.  Note  that,  the  peak  1!  occurs  before  the  peaks  in  D  and 

,)  . 

s 


The  question  as  to  whether  or  when  li  reaches  the  saturated  value 
can  be  answered  by  comparing  the  capacitively-l  united  field  of  liquation  .3-23 
or  3-24  with  l.i  .  Thus  saturation  will  occur  during  the  exponential  rise  if 
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(3-21)) 


in  the  example.  figure  3-3  shows  that  this  occurs  well  before  the  peak  o. 
Saturation  is  much  easier  to  reach  after  the  peak  of  the  gamma  flux. 
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The  diffusion  time  of  ground- induced  fields  to  elevated  observers 
will  be  discussed  in  Section  3.9. 

3.7  FIELDS  GENERATED  BY  AIR-GROUND  ASYMMETRY 

The  radial  li  of  Section  3.6  is  gem  rated  everywhere  in  the  air 
in  the  beginning.  However,  the  ground,  Doing  usually  (but  not  alwuys)  a 
better  conductor  than  the  air,  shorts  out  tne  radial  li  applied  to  it.  A 
current  flows  in  the  ground,  and  this  current,  induces  a  horizontal  magnetic 
field,  in  the  direction  perpendicular  to  the  applied  radial  H  (i.c.,  in  the 
azimuthal  direction  around  the  burst  point),  and  horizontal  and  vertical 
electric  fields.  The  induced  horizontal  H-field  approximately  cancels  the 
applied  radial  li  at  the  ground  surface.  In  the  usual  system  of  spherical 
polar  coordinates  with  r  measured  from  the  burst  point,  0  measured  from 
the  vertical,  and  <J>  the  azimuthal  angle,  the  field  components  present  are 
liv,  Hq  ,  and  The  induced  fields  propagate  away  from  their  point  of 

generation,  through  the  conducting  air  and  soil. 

A  detailed  exposition  of  the  ground- induced  fields  is  possible  (see 

References  1-2  to  1-4)  hut  lengthy,  and  will  not  bo  given  here.  Qualitative 

understanding  can  be  had  by  recognizing  three  phases.  The  wave  phase  occurs 

at  early  times  when  the  air  conductivity  is  negligible  or  small.  In  this 

phase  the  ground- induced  fields  propagate  through  the  air  rather  freely, 

hut  with  some  attenuation  duo  to  the  small  conductivity.  Because  the 

radial  I!  applied  to  the  ground  appears  to  move  outward  with  the  speed  of 

light  (the  speed  of  the  gamma  flux),  the  ground  induced  fields  radiate  chiefl) 

outward  in  a  small  angular  zone  (in  0)  just  above  the  ground.  The  fields 

are  predominantly  If,  and  B,  ,  and  thev  have  the  relation  li  ,  ^  cB, 

0  <|>  '  0  <|> 

appropriate  for  u  vertically  polarized  transverse  wave  propagating  along  the 
ground,  li  is  small  at  the  ground,  but  rises  to  the  radial  li  of  Section 
3.6  at  the  upper  edge  of  the  small,  angular  zone.  For  a  gamma  flux  rising 
as  exp(at),  li  and  B  rise  first  as  exp  (ext),  then  as  exp(u't)  where 


a'  is  a  small  fraction  of  a,  and  then  as  exp(at/2),  the  differences  being 
due  to  the  effects  of  air  conductivity  between  the  burst  point  and  the 
observer.  The  last  time  dependence  is  the  most  important,  since  it  goes  with 
the  largest  field  amplitudes. 

The  wave  phase  ends  and  the  diffusion  phase  begins  when  a  reaches 
the  value  indicated  in  liquation  3-29.  This  is  the  time  when  the  radial  D 
saturates  (without  ground  effects)  and  ulso  the  time  when  the  conduction 
current  exceeds  the  displacement  current.  In  the  diffusion  phase  the  term 
Sli/cDt  can  be  dropped  from  the  Maxwell  liquation  1-18.  The  two  liquations  1-17 
and  1-18  then  define  a  diffusion  problem  like  that  in  the  well-known  skin 
effect,  as  will  be  discussed  in  Section  3.9.  True  wav<'  propagation  ceases, 
due  to  conductivity.  11^  continues  to  increase  as  exp  (at/ 2)  and  li^ 
decreases  as  exp  (-at/2),  15  is  limited,  at  the  ground  surface,  by  the 

finite  ground  conductivity,  and  Is  nowhere  greutor  than  li ^ .  The  ground- 
induced  fields  diffuse  up  into  the  air  until  they  can  go  no  further,  i.e., 
until  the  skin  depth  in  the  air  is  comparable  with  the  radius  r  from  the 
burst . 


When  the  diffusion  is  complete,  the  quasi.stnt.ic  phase  begins.  In 
this  phase,  the  deposition  of  charge  by  the  Compton  current  is  balanced  by 
removal  of  charge  by  the  conduction  current  driven  by  the  statue  electric 
field.  The  electric  field  is  thus  derivable  from  a  potential  <!>, 

H  =  _  V<!>  ( V  x  it  -  0 )  ,  (3-30) 

and  the  conservation  of  charge,  liquation  1-7,  becomes 

V  ■  (I.  -  oV<l>)  =  0  .  (3-31) 

With  .t  and  o  given,  this  equation  determines  <I>,  from  which  H  can  be 
computed.  The  magnetic  field  is  then  determined  by  the  static  form  of 
liquation  1-18, 
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The  fields  in  the  qua s i s t at i c  phase  are  not  exactly  static,  of  course,  but 
the  changes  are  so  slow  that  the  time  derivative  terms  in  Maxwell's  equations 
are  small  compared  with  other  terms.  The  correct  fields  at  each  time  are 
near  the  static  solution  for  the  instantaneous  .5^  and  o.  The  approximate 
solution  of  liquations  3-31  and  3-32  will  be  discussed  in  Section  3.12. 

One  additional  point  must  be  made  regarding  the  diffusion  phase. 

At  positions  sufficiently  close  to  the  burst,  the  peak  air  conductivity  o 

il 

exceeds  the  ground  conductivity  o  .  Since  soil  conductivities  are  usually 
not  much  greater  than  about  10  "  mho/m,  this  is  true  in  the  ease  of  figure 
3-3.  This  condition  i|iodifies  the  diffusion  problem  to  some  extent:.  When 
lVu  °o '  *•).  '•  |:,s  t'igl'l  down  to  the  ground  surface.  In  addition,  the 
Compton  current  in  the  top  layer  of  the  ground  (the  top  half  meter  or  so) 
becomes  an  important  source  of  fields,  as  will  be  discussed  in  Section  3-11, 

figure  3-5  shows  the  fields  in  the  air  just  above  the  ground  as 
functions  ol  time.  These  fields  are  consistent  with  the  sources  in  the 
example  discussed  in  this  chapter.  Although  they  have  not  been  obtained  in 
detailed  calculations,  they  will  suffice  for  our  development  of  source-region 
coupling  theory. 

In  the  remaining  sections  of  this  chapter,  somewhat,  more  detailed 
discussions  of  the  phases  are  given.  In  reading  these  sections,  it  will  he 
useful  to  refer  to  figure  3-5  for  illustration  of  the  features  deduced  or 
stated . 
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3.8 


THE  WAVE  PHASE 


At  oarly  times  the  effect  of  the  air  conductivity  is  either 
negligible  or  small.  At  these  times  the  ground-induced  fields  are  well 
represented  by  outgoing  spherical  waves.  While  both  outgoing  and  ingoing 
wuves  arc  generated  by  the  Compton  current  (in  the  presence  of  the  ground 
asymmetry) ,  the  fact  that  the  Compton  current  movos  outward  with  the  speed 
of  light,  maintaining  approximate  phase  with  outgoing  liM  waves,  causes 
the  outgoing  waves  to  be  built  up  to  amplitudes  much  larger  than  those  of 
the  ingoing  waves.  The  effect  of  the  air  conductivity  is  to  attenuate  the 
outgoing  waves  to  some  degree.  In  the  wave  phase,  the  air  conductivity  is 
generally  small  compared  with  the  ground  conductivity,  und  in  first  order 
the  ground  may  be  regarded  as  u  perfect  conductor.  The  finiteness  of  the 
ground  conductivity  loads  to  some  additional  attonuution  of  the  outgoing 
waves . 


The  relation  between  it  und  it  and  the  attenuation  of  the  out¬ 
going  waves  can  be  understood  by  considering  plane  waves  in  a  conducting 
medium.  Un  the  assumption  that  li  and  B  have  the  forms 


lt.lt  =  Ot(),lt())o-I(wt"k"‘10 


(3-33) 


y  y  .  ,  , 

where  are  constants  denoting  amplitude  and  polarization,  liquations 

1-17  and  1-18  become,  for  freely  propagating  waves  (.)=()), 


k  x 


(3-34) 


(  -  JV'l! 


o'  '"o  '  ck  x  B(> 


(3-33) 


Crossing  k  into  (3-35)  and  using  (3-34.)  to  eliminate  k  x  li(1  leads  to: 


a)  j  i :oi 

c  '  c 


:o0)l5(, 


k  x  (kx|t()) 

y  y  y  .■>  v  -  y 

k ( k • B^ )  +  ( k • k ) 


(3-3(,) 
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The  initial  condition  (1-3)  becomes,  for  the  present  case, 


k  •  !$()  =>  0 


(It  is  necessary  here  to  use  the  initial  condition  since  the  assumed  fields 
are  oscillatory  at  all  times,  rather  than  starting  from  zero.)  Therefore 


k  ■  k  ■  ko“-j  • 


(3-37) 


where  k(J  is  defined  as  the  propagation  constant  in  the  absence  of  conductivity, 


kn  s  VT  W  . 

0  c 


(3-38) 


liquation  3-37  indicates  that  for  non-vanishing  a,  k  will  be  complex,  with 

■  ►  y 

real  und  imaginary  parts  and  k^, 

k  «  kr  +  jk L  .  (3-39) 

As  will  lie  seen  in  Section  3.10,  liquation  3-37  does  not  require  that  the 
directions  of  k^  and  k,  lie  the  same.  If  they  are  arbitrarily  chosen 
to  bo  the  same ,  so  that 

'£  =  nk  ,  (3-40) 

where  n  is  a  real  unit  vector  and  k  is  a  complex  number,  then  liquation 
3-37  allows  n  to  be  arbitrary  but  determines  k, 


k  b  k. 


i  •  2°<5- 

VTk 


(3-41) 


When  the  magnitude  of  the  imaginary  term  in  the  radical  is  small  compared 
with  unity,  this  solution  is  approximately 


Z  o 

k  w  k0  “  J  2VT 


(3-42) 


, n  the  air,  c  =  1,  and  the  attennat ion  length  l  is  independent  of  frequency, 


(il'ICWLf. 


(air) 


(3-43) 


(>  =  2/Z()o  . 

-y  -y  -y  y 

The  directions  of  li  ) » U  ^  and  n  arc  mutually  orthogonal,  with  f ^  *  B() 
in  the  direction  of  n.  from  liquation  3-34  it  follows  that  the  complex 


amplitudes  of  £  and  I5()  are  related  by 

C»0  YT  (l-J  B0  . 


(3-44) 


Since  the  imaginary  term  here  lias  been  assumed  to  be  small  compared  with 
unity,  the  magnitudes  of  li  and  Bn  are  related  approximately  by 

W  U 


y.  S 


(3-45) 


.in  an*  . 


in  the  earliest  part  of  the  wave  phase,  o  is  negligible  every¬ 
where.  In  this  case,  the  fields  will  rise  us  eut  if  the  Compton  current 
,) ^  does;  hence  the  name  a  w_a v e  j)lia^S£  for  this  regime.  The  ground-induced 
fields  are  I!,,  and  11,  in  the  spherical  coordinates  indicated  in  figure 

1)  <|J 

3-(>,  and  these  fields  are  produced  by  the  shorting  out  of  li  at  the 
ground  surface. 

The  exponential  rise  of  the  fields  ceases  when  the  attenuation 
becomes  important,  anywhere.  Since  o  is  largest  very  near  the  nuclear 
device,  attenuation  first  becomes  important  at  r  =  r(),  the  radius  of  the 
device  itself.  The  transition  occurs  when  the  attenuation  length  « 
becomes  loss  than  r  ,  i.e.,  when 


0(V 


(3-4(0 
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z 


Figure  3-6.  Cartesian  and  spherical  coordinates.  The  x,y 
plane  is  the  air-ground  interface. 


This  condition  is  reached  very  early  in  the  rise,  and  it  is  difficult  to 
detect  the  a  wave  phase  at  appreciable  distances  from  the  explosion. 


Let  an  observer  be  located  at  radius  r  »  r^.  If  r  is  not 
too  large,  attenuation  will  eventually  become  important  at  this  distance. 

At  most  distances  of  interest,  the  dominant  variation  of  o  with  distance 
(at  constant  retarded  time)  comes  from  the  factor  e  1 ^  in  Equation  3^20. 
The  time  of  A-saturat 1 on  at  the  observer  is  defined  as  that  when  & 
becomes  less  than  A,  i.e.,  when 


a(r) 


--v-  fj  2  x  10  °  mho/m  at  r  . 

V 


(3-47 ) 


The  condition  (3-46)  is  also  callci  A -saturation  at  r  .  Between  the 
retarded  times  at  which  the  conditions  (3-46)  and  (3-47)  are  reached,  the 

01  *  t 

fields  at  r  rise  roughly  as  c  ,  where  a'  is  a  small  fraction  of  ot. 


After  A-saturat  ion 

at/2  ...  . 

ri.se  as  e  ,  provided  ,J 

continues  until  a  reaches 

called  cx-saturation,  and  is 

displacement  current. 


at  the  observer,  the  fields  IL  and  B, 
is  still  rising  as  c  .  This  dependence 
the  value  indicated  by  liquation  3-29,  which  is 
the  condition  that  conduction  current  exceed 
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The  end  of  the  wave  phase  occurs,  for  a  given  observer,  at  the 
time  of  a-saturation  at  his  location.  During  the  entire  wave  phase,  the 
relation  (3-45)  holds  between  Eg  and  B^.  These  fields  are  confined 
chiefly  to  a  layer  of  air  just  above  the  ground,  with  thickness  5  of  the 
order  of  a  few  meters.  At  the  onset  of  a-saturation 

»  6  w  c/  ■  (3-48) 

and  the  fields  are 

“  >'-e  “  -  Bs  '  (,-49> 

Here  E  is  the  saturated  field  defined  in  Section  3.6. 
s 

3.9  THE  DIFFUSION  PHASE 


The  diffusion  phase  begins  when  the  air  conductivity  reaches  the  value 

given  by  liquation  3-29,  and  in  it  the  displacement  current  is  negligible  in  the 

air.  The  dominant  fields  arc  li  and  B.  .  li  is  near  the  saturated  field 

r  i|>  r 

Es  except  in  a  layer  just  above  the  ground,  and  B^  is  appreciable  only  in 
tliis  layer.  In  this  thin  layer  it  is  convenient  to  use  the  Cartesian  coordin¬ 
ates  indicated  in  Figure  3-6.  I'or  an  observer  located  on  or  near  the  y  axis, 
the  y-dircction  is  approximately  the  same  as  the  r-direction  and  the  x-direction 
is  approximately  the  same  as  the  -^-direction.  The  distance  above  the  ground 
is  z.  Maxwell's  Equations  1-17  and  1-18  become,  for  this  case 


4  B  ,  Ah 

_.i  = 

At  Ac 


(3-30) 


’Oh  -•  -  Z  ■)  -  c 

Or  Or 


Substituting  l;  from  the  second  equation  into  the  first  gives 


OB  / . J  \  .  .  i) B . 

.A  =  ?_  In  +  .)  l  a 

At  be  '  o  I  ZiS  Az  o  Ac 


(3-52  ) 
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This  is  a  type  of  diffusion  equation  for  the  magnetic  field.  The  first 
term  on  the  right-hand  side  is  the  source,  for  without  it  =  0  would 

be  a  solution.  Both  and  a  are  approximately  independent  of  z  in 

the  air,  so  the  source  exists  only  at  the  ground  surface,  where  a  changes 
from  air  to  ground,  and  in  the  top  10  to  30  cm  of  the  ground,  where 
falls  rapidly  due  to  attenuation  of  the  gamma  flux  in  the  ground.  If  the 
ground  conductivity  is  large  compared  with  the  air  conductivity,  most  of  the 
source  occurs  at  the  ground  surface. 


The  total  magnetic  flux  ,‘f  per  radial  meter. 


=  J  B^dz  ,  (3-53) 

produced  by  the  source  can  be  found  by  integrating  Equation  3-52  over  z. 
Since  B,  and  3B,( 

4*  (, 

deep  in  the  ground, 

~  =  (— J  =  -  E  .  (3-54) 

3t  \  a  ' 


Since  B^  and  SB^/Bz  vanish  deep  in  the  air  and  ground  and  Jr/o  vanishes 


ai  r 


Thus 


•Xt)  =  .'f 
s 


L 

•:  -  J!  - '  ■ 


(3-55) 


where  -'f  is  the  flux  at  the  time  t  of  the  beginning  of  the  diffusion 
s  s 

phase  (u-saturation)  .  (Note  that  and  .>  are  negative,  i.e.,  B({)  runs 


clockwise  around  the  burst  point.) 


of  z . 


In  the  air  the  source  term  vanishes  and  o  is  (roughly)  independent 
Thus  Equation  3-52  becomes 


3B. 

_<iL 

at 


3“B, 
_c_  f 

“0°  3  z“ 


(3-56) 
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This  equation  governs  the  diffusion  of  up  into  the  air.  The  solution 

is  qualitatively  different  during  the  exponential  rise  and  after  the  peak 
of  the  gamma  flux.  During  the  rise,  let  us  try  the  assumption  that 

~  eBt  ,  (3-57) 

where  3  is  a  constant  to  be  determined,  liquation  3-56  then  determines 
the  z-dependence, 

B.  ~  c  ,  ( o - 5 .  ) 

V 

where  the  skin  depth  6  is 


Z,.o6  \  Pno3 


(3-59) 


Since  a  increases  as  c  ,  6  decreases  during  the  rise  as  c 
amplitude  of  can  be  estimated  from  5  and  the  total  flux 

B,3  &  .f  ,  or  B.  » .'S/6  . 

(|)  .)> 


-ut/2  .... 

I  he 


(3-60) 


Since  li  is  constant  during  the  exponential  rise  of  the  gamma  flux, 
liquation  3-55  indicates  that  ( - )./  increases  only  linearly  with  time,  or 
slowly  compared  with  the  exponential  increase  of  1/6.  Thus  approximately 

B ,  -  e’lt/J  .  (3-61) 


Comparison  with  liquation  3-5"’  shows  that 


vX’ 


(3-62! 


Note  that  in  writing  liquation  3-60,  we  neglected  the  flux  in  the 
ground.  This  is  permissible  if  the  skin  depth  in  the  ground  is  small  compared 
with  that  in  the  air  (so  that  the  ground  contains  little  flux),  or  if  *iie 
ground  conductivity  is  large  compared  with  the  air  conductivity.  The  case 
in  which  this  condition  does  not  hold  will  be  discussed  in  Section  3.11. 
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After  the  peak  in  the  gamma  flux,  o  falls  and  t he  skin  depth 
increases.  Note  that  liquation  3-56  can  be  brought  to  simpler  appearance 
l)>'  changing  the  time  variable  to  T  defined  by 


dT 


cdt 

V 


cdt 

V 


(3-63) 


where  t. 

P 


is  the  time  of  the  peak,  liquation  3-56  then  becomes 
3T  "  3, 2  ' 


(3-64) 


Solutions  of  this  equation  can  be  found  as  functions  of  the  similarity 
variable  z/VT.  Thus  the  skin  depth  is 


(S  re  VT  = 


(3-65) 


The  amplitude  of 


can  again  be  estimated  from  liquations  3-55  and  3-60. 


figure  3-4  shows  that  I:  ,  after  falling  about  one  decade  from  its 
peak,  is  then  almost  constant  for  several  decades  in  time.  During  most  of 
this  interval,  liquation  3-55  becomes 

./( t )  re  -  !■;  t  .  (3-66  ) 
s 


Figure  3-3  shows  that,  to  a  crude  approximation,  o  can  be  written  over  the 
same  time  interval  as 


a(t ) 


(3-67) 


where  the  subscript  p  indicates  peak  values.  With  this  approximation,  the 
skin  depth  becomes 


is 


From  Equation  3-60,  the  estimate  of 


<l> 


(3-69) 


Thus  B^  is  roughly  constant  in  time  in  the  diffusion  phase  after  the  peak 
of  the  gamma  flux. 


As  stated  before,  13  tends  to  be  small  at  the  ground  surface  and 
rises  to  E  at  heights  of  a  few  skin  depths.  If  the  ground  conductivity 
is  not  very  large  compared  with  the  air  conductivity  o&,  then  E^  at 
the  surface  is  given  approximately  by 

VT“ 

0  w  .  (3-70) 

1  Vo  +  V~0~ 

a  g 

This  formulae  comes  from  considering  the  impedances  of  air  and  soil  within 
one  skin  depth  from  the  surface.  The  return  conduction  current  that  would 
flow  in  one  skin  depth  in  the  air,  if  the  ground  were  not  present,  is  shared 
with  one  skin  depth  in  the  ground. 

The  field  Eg  would  be  small  in  the  diffusion  phase  if  it  were 
not  for  the  effect  of  the  magnetic  field  on  the  Compton  current.  The 

Compton  electrons  are  turned  upwards,  away  from  the  ground  by  the  magnetic 

•  -  3  2 

force  on  them.  For  greater  than  about  10  '  Weber/m  =  10  gauss,  the 

resulting  Jg  is  comparable  with  J^.  A  roughly  static  field  Eg  then 

arises,  of  sufficient  magnitude  to  drive  a  conduction  current  cancelling 

Jg.  Thus  E0  is  comparable  with  Es ,  provided  is  as  large  as  indicated 

above . 


3.10 


FIELDS  IN  THE  GROUND  AT  EARLY  TIMES 


Cables  and  other  components  of  systems  are  often  buried  at  depths 
of  one  to  a  few  meters  in  the  ground.  It  is  therefore  important  to  see  how 
fields  propagate  in  the  ground.  The  assumption  of  oscillatory  fields  of 
the  form  (3-33)  leads  to  Equation  3-37,  which  can  also  be  wrr  ^en  as 


The  factor  in  parentheses  here  is  -jn,  where  n  is  the  relative  admittance 
of  the  soil  defined  in  Chapter  2.  figure  2-5  shows  that  at  the  higher 
frequencies  of  interest,  neither  the  real  or  imaginary  part  of  r]  is 
negligible . 


The  fields  in  the  ground  can  be  related  to  the  magnetic  field 
at  the  surface,  discussed  in  previous  sections.  Over  distances  of  only 
several  meters,  B^  at  the  surface  can  be  regarded  as  a  function  only  of 
t  -  i.c.,  except  for  time  delay,  B^  is  the  same  at  different  r.  Thus 

the  radial  phase  velocity  of  all  Fourier  components  of  B^  is  c,  so  that 
the  component  with  frequency  w  has  radial  wave  number 

kr  =  oj/c  .  (3-72) 


Since  k  has  the  two  components  k  and  k 


and  since  k  in  the  ground 
r  1 


must  match  that  of  B.  at  the  surface,  Equation  3-71  determines  kr 


. ,  oo 
±k  =  — 
z  c 


as 


(e-1)  -  j 


zo°c 

(0 


(3-73) 


Thus  k^  is  complex,  so  that  the  waves  attenuate  as  they  propagate  downward 
in  the  ground. 


The  fact  that  k  is  complex  while  k  is  real  means  that  the 
z  1  r 

real  and  imaginary  parts  of  k  do  not  have  the  same  direction.  The 
dispersion  relation  (3-71)  docs  not  force  real  and  imaginary  parts  to  be 
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parallel,  i.c.,  phase  planes  and  amplitude  planes  need  not  be  parallel, 
liquation  3-71  determines  one  Cartesian  component  ot  k  if  the  other 
components  are  specified. 


The  real  and  .imaginary  parts  of  arc  plotted  as  the  points  in 

Figure  3-7  as  a  function  of  tn  for  our  standard  soil,  for  which  t ho  relative 
admittance  was  graphed  in  Figure  2-5.  Messier  (Reference  3-1)  noticed  that 
a  remarkably  good  fit  to  soil  propagation  constants  is  obtained  by  the 
simple  formula 


(3-74) 


The  curves  in  Figure  3-6  represent  the  real  and  imaginary  parts  of  this 
formula,  which  obviously  has  the  same  limits  for  low  and  high  <n  as 
liquation  3-73.  The  values  of  f.  and  o,.  used  in  the  fit  are 

cu)  =  6.5  ,  o(|  =  H  x  11)"^  mho/m  .  (3-75) 


Note  that  o',  lias  the  same  value  as  in  Table  2-1,  whereas  i:  has  been 

u  7)  ,  9  - 1 

adjusted  slightly  to  give  a  better  fit  over  the  range  10  10  sec 


if  the  magnetic  field  B^(t)  at  the  ground  surface  is  represented 
by  its  Fourier  transform  B  (<n) ,  then  B,  in  the  ground  at  depth  c  (taken 

tl  (p 

positive)  is  obtained  by  propagating  each  frequency  component  with  its  k_, 


B.U.t)  =  •“•  /  B  too )  exp  |  j  (wt-k  -/.)  |  dm 

cp  Z 'I  J  u  u 


(3-76) 


If  Messier's  approximation  for  k_  is  used,  this  equation  can  he  written  as 


B|)U,t)  =  2~-J  B;i fo>) exp |  jiut '  - 


j  (»T,  jdi.i  , 


(3-77) 


where 
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u 

O) 

00 


3 
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Figure  3-7.  Propagation  constants  in  standard  soil.  Points:  real  and 
imaginary  parts  computed  from  exact  Equation  3-73.  Solid 
curves:  Messier's  approximation.  Equation  3-74.  Dashed 
line:  assumed  k  =  cj/c. 


n..nuki-i,')i,u..i  i,ni;i.'j.  Tprrn^rrmTwrrr: 


and 


Ti  =  -  1  z/c  ■ 


5  Vo2^ 


(3-78) 


(3-79) 


The  appearance  of  t'  in  the  integral  has  the  result  that  the  onset  of  the 

field  at  depth  z  is  delayed  by  the  time  after  onset  at  z  ~  0.  T. 

is  the  time  for  the  highest  frequencies  to  propagate  a  distance  z  in  the 

vertical  direction.  (The  actual  phase  propagation  direction  is  not  vertical, 

but  is  in  the  direction  corresponding  to  the  components  k  ,  k  .) 

r  2 1 


The  radical  in  the  exponential  in  liquation  3-77  gives  u  diffusive 

spreading  in  time,  in  addition  to  the  delay.  This  can  be  seen  if  B  (t) 

u 

is  taken  as  an  impulse  function  at  t  ~  0,  for  which  B^fuQ  =  1.  Evaluation 
of  the  integral  then  gives 


l^U.t) 


exp(-T2/4t ' ) 


(3-80) 


It  can  be  seen  that  the  time  integral  of  this  function  is  unity,  independent 

2 

of  *1*2 .  A  graph  of  is  presented  in  Figure  3-8.  Since  T9  ~  z“, 

the  impulse  function  at  z  =  0  is  broadened  into  a  longer  and  lower  pulse 
with  increasing  z.  For  our  standard  soil 

T  =  1.01  x  10"8  z2sce  .  (.3-81) 

_  ^ 

Thus  at  1  meter  depth,  the  impulse  response  is  a  pulse  of  approximately  10 
second  duration.  For  arbitrary  B() (t)  ,  the  response  at  depth  z  can  be 
obtained  by  folding  13 (t)  with  the  impulse  response. 


it  is  also  instructive  to  examine  the  z-depcndencc  of  the  impulse 
response  (3-80)  for  fixed  t.  For  this  purpose  it  is  convenient  to  define 
a  skin  depth 
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v. 


j 


T-B.(z,t) 


fi(t)  = 


(3-82) 


and  the  maximum  depth  reachable  by  waves  in  the  ground  in  time  t 


z  (t)  =  ct/VT  -  1 

nr  •  oo 


(3-83) 


With  these  definitions ,■ Equation  3-80  can  be  written  as 


B,.(z,t)  =  - - 

<|)  Vvt  (1  -  — ) 

z 

111 


z  /  5 

'z~372  exP 


(z/6)" 

.  (1 


(3-84) 


This  formula  contains  two  characteristic  lengths,  and  so  is  not  graphublc  as 


a  single  curve.  The  two  lengths  cS  and  z  arc  equal  when 


4(c  -1)  o 

ct  =  — ~ - w  7.3  m  ,  or  v  »  2.43  x  10  sec  .  (3-85) 

V() 


For  smaller  t,  z  <  6.  and  for  larger  t,  z  >  6.  iiic  numerical  values  here 
’in  h  m 


arc  from  liquation  3-75  for  Messier's  fit.  Figure  3-9  shows  tB, (z,t)  for 


1  T  _  O 

z,n  =  2  6  >  5*  and  26,  corresponding  to  t  =  O.bl,  2.43,  and  9.72  *  10 


second.  Also  shown  is  the  limiting  case  for  z  >>  6.  In  this  presentation 


evidence  of  propagation,  as  contrasted  with  diffusion,  practically  disappears 
by  the  time  z  >26. 


If  instead  of  B<j>,  li  at  the  surface  (or  I!  )  is  specified  as 


H  (t)  or  li  (a>)  the  same  formulae  can  be  used  to  obtain  li  (or  ii  ) 


at  depth.  Thus  these  fields  arc  also  broadened  in  time  and  decreased  in 
amplitude  with  increasing  depth. 


The  determination  of  B.  ,  li  and  li  at  the  surface  must  be 

<(>  r  z 


accomplished  by  solving  Maxwell's  equations  in  the  air  and  ground 
simultaneously.  The  analysis  of  this  section  does  not  avoid  that  work, 
but  only  explains  the  relation  of  the  fields  in  the  ground  to  those  at  the 
surface.  The  analytical  method  developed  here  could  be  used  to  determine 


the  ground  fields  if  the  ground  conductivity  were  always  large  compared 

with  the  air  conductivity.  In  this  case,  B  (t)  could  lie  determined  by 

a 

solving  Maxwell's  equations  in  the  air  over  perfectly  conducting  ground. 

Next,  could  be  found  in  the  ground  by  the  methods  of  this  section, 

and  li  and  ti  determined  from  B,.  It  is  more  convenient  to  obtain 

r  z  <p 

the  fields  by  use  of  the  finite-difference  codes  LFiMP-SUBL. 

3.11  EFFECT  OF  COMPTON  CURRENT  IN  THE  GROUND 

The  attenuation  length  of  the  gamma  rays  in  soil  is  about 

\  ,  fa  0.2  m  .  (3-86) 

The  relation  of  this  length  to  that  in  air,  liquation  3-21,  is  determined 
.  3 

by  the  density  of  soil,  about  2  gm/enf  ,  as  compared  with  that  of  air,  about 
-  3  3 

1.23  x  10  '  gm/enf  .  It  i.s  clear  that  any  gammas  observed  in  the  ground  at 

appreciable  distances  from  the  burst  must  have  traveled  mostly  through  the 

air  and  entered  the  ground  only  near  the  point  of  observation.  If  the 

burst  were  on  a  flat  surface  and  there  were  no  air  scattering,  the  prompt 

gamma  flux  would  drop  to  zero  at  the  air-ground  interface.  Gammas  scattered 

in  the  air  can  enter  the  ground,  but  arrive  with  a  time  delay  corresponding 

“  8 

to  their  longer  path.  Since  the  prompt  pulse  is  only  a  few  times  10 
second  in  width,  scattered  paths  that  are  longer  by  more  than  about  10  meters 
than  the  direct  pat.li  do  not  contribute  to  the  prompt  pulse  at  distance. 

Only  gammas  scattered  through  small  angles  can  contribute,  and  the  number 
of  these  is  only  about  5  percent  of  the  unscattered  gammas  during  the 
prompt  pulse.  Thus  for  a  burst  on  a  flat  surface,  the  gamma  flux  drops  by 
a  factor  of  about  20  at  the  air-ground  interface  and  decays  further  to 
negligible  values  in  depths  of  the  order  of  10  cm.  Since  the  ratio  of  Compton 
current,  density  to  gamma  flux  is  approximately  independent  of  material  and 
density  (see  liquation  3-4),  the  total  Compton  current  in  the  ground  is  equal 
to  that  in  only  about  10/20  =  0.5  cm  of  air  above  the  surface.  The  Compton 
current  in  the  ground  is  negligible  during  the  prompt  pulse  for  a  burst  on 
a  flat,  surface. 
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If  the  burst  is  above  the  surface,  or  if  the  surface  is  curved  and 
is  exposed  to  line  of  sight  from  the  burst  at  the  observer,  unscattered 
gammas  can  enter  the  ground.  If  the  angle  between  the  line  of  sight  and 
the  surface  tangent  is  vp ,  as  in  Figure  3-10,  then  the  Compton  current 
density  in  the  ground  attenuates  with  depth  approximately  as 


-> 

J 


s 


JsOe 


-  z  /d 


(3-87) 


where 

d  =  A  sinip  =  0.2  sini|>  m  .  (3-88) 

8 

Mere  is  the  Compton  current  density  in  the  air-ground  interface.  The 

current  density  is  (approximately)  continuous  across  the  interface,  but 
decays  in  the  ground  in  a  depth  d . 


The  Compton  current  in  the  ground  is  important  when  the  air 
conductivity  exceeds  the  grou  d  conductivity.  At  distances  of  interest,  this 
happens  only  during  the  prompt  pulse.  In  this  case,  the  electric  field  in 
the  air  is  limited  by 


Figure  3-10.  Geometry  and  coordinates  for  ground  field 
analysis . 
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2  x  io"  V/m  . 


In  the  ground,  if  E  reached  the  value  Js/a^,  it  would  be  larger  than  that 
in  the  air.  The  field  in  the  ground  does  not  generally  reach  this  value, 
because  the  return  conduction  current  flows  over  a  thickness  of  one  skin 
depth  6  of  soil,  and  6  is  usually  larger  than  d. 

The  fields  produced  in  the  ground  can  be  calculated  exactly  for 
a  gamma  flux  rising  as  exp(at),  as  is  done  in  Reference  3-2,  under  the 
assumption  that  a  >>  a  so  that  the  air  can  be  regarded  as  a  perfect 
conductor.  The  results  of  that  calculation  are  summarized  briefly  here. 

The  relative  admittance  nC°0  0"  *)  has  been  defined  in  Section 
2.2,  and  Figure  2-5  shows  n(ot)  for  our  standard  soil.  The  a  ual 
admittance  Y  =  n/Z^  (mho/m) .  The  skin  depth  in  the  soil  is 


*•'  0.4m  . 


(3-90) 


The  numerical  value  here  is  for  the  example  presented  in  this  chapter. 

If  the  conduction  current  returned  on  the  same  paths  followed  by 
the  Compton  current,  the  electric  field  in  the  ground  would  be 


E  0  -  1 

v  n  s 


(3-91) 


Instead  of  this  relation,  Reference  3-2  shows  that  the  maximum  horizontal 
component  E  in  the  ground  is 


0  sOy  fi2  ’ 


(3-92) 


where  .J^  is  the  y  component  of  the  Compton  current  density  at  the  surface 
This  value  is  reached  at  depth  ~d  in  the  ground;  F.  =  0  at  the  surface  in 


L,.;  •*.«  ■  t\ j&iifas- i-Yf Liistii 


this  model  (large  air  conductivity)  because  the  field  produced  in  the 

■>  i 

air  was  neglected.  (One  could  add  1;  (air).)  The  factor  d“/iS“  in 
liquation  3-92  comes  from  two  sources.  One  factor  d/<S  comes  from  the 
fact  that  the  return  current  flows  over  depth  6  while  the  Compton  current 
flows  over  depth  d  <  <S.  The  other  factor  comes  from  the  proximity  of 
the  Compton  current  to  the  highly  conducting  air.  By  using  liquation  3-89 
to  relate  ii  to  the  saturated  field  !i  in  the  air  and  liquations  3-88 

y  s 

and  3-90,  liquation  3-92  can  be  written 

a  a  ■> 

li  fe  ii  n  7  -  A  sin  ilicosil/  , 

ym  s  a  0  c  g 

.  2 

w  5s m  (jicosi^li^  (example) 

In  the  example,  the  peak  a,  was  used  from  Figure  3-3.  The  maximum  value  of 

2  .  '* 
sin  i(icosi|j  is  about  0.3b,  at  T  =  63°. 


The  maximum  value  of  li„  in  the  ground  occurs  just  below  the 
surface,  and  is 


•7 


1! 

s 


s  i  mji  , 


r-  lDsiniJili  (example) 


(3-94) 


Skin  depth  spreading  does  not  reduce  the  vertical  conduction  current  density 
or  li  .  This  large  ii.,  decays  in  the  depth  d. 


The  magnetic  field  B  (=-B. )  is  reversed  from  the  usual  direction 

x  <!> 

near  the  surface,  because  most  of  the  Compton  current  in  the  ground  returns 
as  conduction  current  in  the  air.  The  value  of  15  at  the  surface  is 

'\x0  “  C'  lUsy()  =  "oV'sO^n^cos'F  >  j 

(3-95) 

-  5  *  10 "°s  i  nificosij;  Weber/m“  (example)  .  ' 
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in  the  example  the  peak  current  density  was  used  from  Figure  3-2.  This  result 
may  be  compared  with  the  field  pre  .cod  l.y  the  air  current  at  the  ground 
surface , 


*\x()a  ~  ‘'o^a^syO  ~  U0 


.COSIp  , 


-2.2  x  10  cosif  Webcr/nf 


sO 

i  example ) 


(3-96) 


The  field  11  decays  in  depth  d  in  the  ground,  and  changes  sign  again 
due  to  return  conduction  current  flowing  below  the  Compton  current. 


When  the  air  conductivity  falls  below  the  ground  conductivity, 

after  the  peak,  6  soon  becomes  larger  than  6  ,  and  the  ground  current 
1  g 

and  fields  are  again  dominated  by  sources  in  the  air. 


Two  cautions  should  lie  noted  regarding  the  results  of  this  section. 
First,  the  peak  II  _  estimated  in  liquation  3-94  is  in  a  range  that  might 
lead  to  breakdown  in  the  soil.  Second,  the  peak  dose  rate  in  the  example 
is  si.if  fi  c  i  ent  ly  high  that  it  might  increase  the  conductivity  of  the  top 
lay*.  'f  the  ground.  Not  much  can  be  said  with  certainty  about  the 
probability  of  occurrence  of  either  of  these  effects;  experiments  with 
relevant  soil  samples  arc  needed.  Both  effects  could  be  expected  to  result 
in  a  decrease  in  the  fields  at  greater  depths  in  the  ground,  but  rely  ng  on 
this  expectation  might  be  risky. 


Attention  was  first  called  to  the  importance  of  the  ground 
Compton  currents  by  R.  R.  Schaefer  md  W.  R.  Cruivm. 

3.12  THE  QUASiSTATIC  PHASE 


As  stated  in  Section  3.,’,  r.  h.  quns  istat  i  c  phase  begins  when  the 
skin  depth  in  the  air  becomes  as  1  ;(/<;:  as  allowed  by  the  spherical  geometry, 
i.e.,  when  '  exceeds  the  distance  r  from  the  burst  point, 
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5  >  r . 


(3-97) 


The  approximation  (3-67)  for  o(t)  is  representative  if  the  time  origin  is 
suitably  chosen,  as  in  Figure  3-3.  In  that  case  (our  usual  example) 

0^  =  0.5  mho/m  ,  t  6  *  10  ^  sec  .  (3-98) 

Equation  3-68  then  gives  the  skin  depth 

6  ^  3.6  x  106  t  (example)  .  (3-99) 

With  r  fa  500  meters,  the  quasistatic  phase  begins  at 

t  =  1.4  x  lo'4  sec  .  (3-100) 


Figure  3-3  shows  that  ion  conductivity  is  beginning  to  be  significant  at 
this  time.  Over  most  of  the  quasistatic  phase,  ion  conductivity  is 
dominant . 

It  was  shown  in  Reference  1-2  that  a  fair  approximation  to  the 
solution  of  the  governing  Equation  3-31  for  the  quasistatic  phase  is  obtained 
by  taking  the  potential  function  <J>  to  be  a  function  of  0  alone.  Thus 
the  predominant  electric  field  is  E„;  the  electric  field  lines  are  approxi¬ 
mately  circular  about  the  burst  point  as  center.  The  fields  Eg  and 
can  be  deduced  from  this  model  and  the  conservation  of  charge,  Equation 
3-31.  If  we  assume  F,  <<  Eg,  this  i.  (nation  becomes  (on  writing  the 
divergence  operator  in  spherical  coordinates) 


_1 _ 

rs  i  n0 


_9 

96 


sinOoEg 


V  •  3 

s 


(3-101) 


Since  n  is  approximately  independent  of  0  and  since  the  distance  scaling 
of  .1  is  approximately  that  in  Equation  3-20,  this  equation  can  be  written 


a  9_ 
rsinO  9 


s  inOH 


n 


.1 

s 

"A  ’ 


This  relationship  is  not  accurately  valid,  but  one  may  compute  the  following 
numbers : 

x  i  r/A  =  0.5,  1,  2,  4 ,  5, 

xe'x/2  =  0.39,  0.61,  0.74,  0.54,  0.41  . 

It  can  be  said  that  the  relation  is  valid  to  ±30  percent  over  the  range  from 
r  =  0.5A  to  5A .  In  much  of  the  quasistatic  phase,  the  hard  gammas  from 
air  capture  arc  strong  contributors,  for  which  the  effective  A  (including 
build  up)  is  about  400  meters.  Thus  the  model  is  justified  over  the  range 
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r  =  0.^  to  2  km.  Hardening  of  the  gamma  spectrum  with  distance  make  the 
relation  3-104  more  accurately  satisfied,  since  X  increases  with  distance. 
Thus  Equation  3-103  gives  reasonably  good  values  for  Li^  over  the  range 
indicated . 


Photographs  of  large  yield  nuclear  explosions  show  lightning-like 
discharges  developing  in  the  time  frame  from  10  '  to  3  x  lu  *•  second 
(Reference  3-5).  The  discharges  rose  from  sharp  objects  (antennas)  on  the 
ground,  and  followed  quite  closely  the  0-dircction.  The  inference  that  they 
were  driven  by  the  quasistatic  electric  field  is  hard  to  resist.  The 
growth  rate  of  the  discharges  contains  information  on  the  magnitude  of 
the  electric  field,  but  analysis  of  the  data  is  not  complete  at  this  writing. 

The  magnetic  field  can  be  found  from  liquation  3-32.  The  r-component 
of  this  equation  is 


— • — rT  sinOB.  =  —  .J 

rsinO  30  (pcs 


(3-105) 


This  equation  has  the  same  form  as  liquation  3-103,  and  the  same  integration 
procedure  leads  to 


13  =--  r.l  tan  f  , 

(\)  e  s  2 


C  T  =  l'„> 


(5-107) 


The  time  dependence  of  1:^  and  13(1  can  be  deduced  from  liquations 
3-103  and  3-107.  lb  is  independent  oi'  time  until  the  conductivity  becomes 


ionic,  then  falls  as  VJ^.  13^  falls  as  throughout  the  entire 

quasi  static  phas  . 


The  radial  field  li^  at  the  ground  can  be  found  from  the  condition 
that  the  entire  Compton  current  passing  out  through  the  hemisphere  at  radius 
r  must  return  in  the  ground.  The  radial  dependence  of  I-  is  complicated 
by  a  rather  complex  current  flow  pattern  in  the  ground,  and  the  time  dependence 


if  affected  by  the  fact  that  diffusion  persists  longer  in  the  ground  than 
in  the  air.  Once  diffusion  is  completed,  the  time  dependence  of  ii  __  at 

t  7 

the  ground  is  the  same  as  that  of  J  .  Diffusion  takes  about  3.5  x  lo 

s 

second  at  r  =  500  meters. 
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CHAPTER  4 

COUPLING  TO  SHORT  BURIED  CABLES 


4.1  INTRODUCTION 

Ground-based  systems  hardened  to  blast  often  have  electrical 
conductors  buried  at  depths  of  one  to  a  few  meters  in  the  ground  for 
mechanical  protection.  Burial  also  affords  some  protection  from  BMP,  but 
by  no  means  total  protection  since  the  fields  penetrate  to  these  depths 
without  strong  attenuation,  especially  in  the  lower  frequency  components. 

It  is  useful  to  distinguish  two  categories  of  buried  cables.  The  category 
considered  in  this  chapter  includes  cables  shorter  than  a  few  hundred  meters, 
i.e.,  one  gamma-ray  absorption  length  in  air.  For  such  cables  the  EMP 
fields  may  be  assumed  to  have  roughly  constant  amplitude  along  the  length 
of  the  cable.  The  phase  is  not  constant,  for  in  general  the  BMP  will  sweep 
over  the  cable  with  a  speed  determined  by  the  angle  x  between  the  cable 
run  and  the  radial  direction  from  the  nuclear  burst.  (See  Figure  4-1.) 

The  phase  speed  v  is 


Figure  4-1.  Geometry  of  burst  and  cable  run. 
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(4-1) 


v  =  c/cosx  ,  cosx  =  cosxecosxa  . 

In  case  the  burst  is  above  the  surface,  the  angle  x  is  made  up  from  two 
angles,  the  elevation  angle  xe  of  the  burst  point  as  seen  from  the  cable, 
and  the  azimuthal  angle  x  in  the  ground  plane.  The  horizontal  LiMP 

electric  field  E^  at  the  wire,  under  the  assumptions  of  constant  amplitude 
and  constant  phase  speed,  has  the  form 


Eh  ~  Eh(t 


(4-2) 


where  z  is  the  distance  along  the  wire  (not  the  vertical  coordinate) .  The 
field  component  of  interest  is  that  parallel  to  the  cable,  which  is  related 
to  the  horizontal  field 
burst  by 


E  in  the  projected  radial  direction  from  the 
pr  1 


Eh  =  C0SXaKpr  •  (4~3) 

The  cable  is  most  strongly  driven  when  x  =  0>  and  since  the  location  of  the 

cl 

burst  will  not  usually  be  predictable,  this  case  should  be  assumed. 


If  E^  is  Fourier  analyzed  in  terms  of  waves  of  the  form 
exp[j  Iwt-kz) ] ,  then  Equation  4-2  leads  to  the  result  that  k  is  deteimined 
by  U! , 

k  =  ~  =  ~  cosy  •  (4-4) 


The  fact  that  v  s:  c  means  that  the  EMP  sweeps  over  the  cables 
faster  than  free  signals  can  propagate  along  it.  Since  also  of  the  soil 

is  considerably  greater  than  unity,  the  EMP  sweep  speed  is  considerably 
greater  than  the  free  signal  speed.  In  particular,  signals  arising  from 
effects  at  the  cable  ends  arrive  at  most  points  along  the  cable  with 
significant  delay  after  the  EMP  arrives.  It  is  therefore  useful  to  calculate 
first  the  response  of  the  cable  ignoring  end  effects,  i.e.,  as  if  the  cable 
were  infinitely  long.  This  problem  is  taken  up  in  Section  4.2. 


The  other  category  of  cables  includes  those  that  are  so  long  that 
the  BMP  environment  cannot  be  assumed  to  have  constant  amplitude  over  the 
length  of  the  cable.  This  category  is  discussed  in  Chapter  5. 

Some  of  the  buried  cables  of  interest  will  be  of  multi-wire  type. 

However,  these  will  generally  have  an  outer  conducting  sheath  to  shield  the 

wires  from  BMP,  lightning  or  other  electrical  interference.  The  conducting 

sheath  will  generally  be  covered  by  an  insulating  sheath  for  protection 

against  corrosion  of  the  (metallic)  conducting  sheath.  In  this  chapter,  only 

the  insulator  and  the  outermost  conductor  will  be  considered.  Transfer 

coupling  from  the  outer  conductor  to  internal  wires,  if  any,  is  a  separable 

problem,  at  least  approximately.  The  cable  geometry  is  defined  by  Figure 

4-2.  The  conductor  is  drawn  as  a  hollow  cylindrical  shell  of  thickness  d, 

but  it  could  be  a  solid  cylinder.  The  electrical  conductivity  of  the 

conductor  is  generally  much  larger,  by  a  factor  of  the  order  of  10^°,  than 

that  of  the  soil.  Thus  the  conductor  has  the  same  resistance  per  unit 

4  5 

length  as  a  cylinder  of  soil  with  radius  of  the  order  10  to  10  times  large)’ 
than  a^,  depending  on  the  thickness  d.  The  relative  permittivity  of  the 


Figure  4-2.  Definition  of  cable  parameters  and 
cylindrical  parameters. 
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insulator  is  a  little  larger  than  unity,  r..  r;  2  being  typical.  The 
conductivity  of  the  insulator,  typically  very  small  compared  with  that  of 
soil,  may  be  enhanced  by  the  gamma  radiation  penetrating  to  the  burial 
depth.  The  peak  gamma  induced  conductivity  may  be  as  high  as  10  ^  mho/m, 
and  is  time  dependent  (proportional  to  the  dose  rate) .  The  analysis  of 
this  chapter  will  treat  a  as  constant,  but  will  estimate  the  importance 
of  gamma-induced  conductivity.  The  soil  parameters  c?  and  will  be 

treated  as  independent  of  time  but  frequency  dependent. 

The  presence  of  the  air,  with  its  time-varying  conductivity,  affects 
the  coupling  to  the  cable.  The  analysis  here  will  first  assume  that  the 
soil  extends  to  infinity  in  all  directions,  and  then  show  how  the  presence 
of  the  air  can  be  taken  into  account  in  an  approximate  way.  In  addition, 
the  wire  initially  will  be  assumed  to  have  perfect  conductivity. 

4.2  PERFECTLY  CONDUCTING  WIRE  IN  INFINITE  SOIL 

HMP  environments  arc  calculated  without  wires  present,  bet  the 
Fourier  component  of  the  environmental  electric  field  parallel  to  the  wire 
axis,  at  the  position  of  the  axis,  with  frequency  ir  be 

lih(w)  =  Fq (oj)  exp |  j  ('.nt-kz)  ]  .  (4-5) 

Here  k  is  specified  by  Equation  4-4.  The  total  parallel  electric  field 
must  vanish  at  the  wire  surface.  A  current  1  flows  in  the  wire,  such  that  the 
additional  fields  produced  by  I  have  a  parallel  electric  field  canceling 
the  applied  iiMP  field. 

The  HMP  electric  field  has  components  both  parallel  and  perpendicular 
to  the  wire  axis.  The  perpendicular  component  causes  the  wire  to  polarize, 
with  positive  charge  appearing  on  one  side  and  equal  negative  charge  on  the 
opposite  side.  The  currents  involved  in  this  polarization  are  small  when 
the  wire  diameter  is  small  compared  with  the  wavelengths  in  the  HMP,  a 
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condition  which  is  generally  well  satisfied.  Further,  these  currents  do  not 
flow  along  the  wire,  so  are  of  no  consequence  for  equipment  connected 
to  the  wire.  The  EMP  field  of  importance  is  the  component  of  the  electric 
field  parallel  to  the  wire  axis,  and  the  variation  of  this  field  around  the 
circumference  of  the  wire  is  negligible. 


The  usual  cylindrical  coordinates  r,  0,  z,  right  handed  in  that 
order  and  indicated  in  Figure  4-2,  are  convenient  for  the  problem  at  hand. 
The  field  components  associated  with  the  wire  current  are  E  ,  E^  and  Bq, 
and  they  are  all  independent  of  0.  Maxwell's  Equations  1-17  and  1-18 
become,  for  fields  with  t  and  z  dependence  given  by  Equation  4-5, 


>B0  =  jkEr  +  |:Ez  , 
nEr  =  jkcB0  , 


(4-6) 

(4-7) 


HE 


-  £  JL  rR 

r  9r  0 


(4-8) 


i 


I 


I 


where  n(w)  is  the  relative  admittance  introduced  in  Section  2.2.  Equation 
4-7  can  be  used  to  express  E  in  terms  of  B0 , 


H 

r 


=  US£  B 
n  0 


(4-9) 


Use  of  this  result  In  Equation  4-6  allows  expression  of  B„ 


B0  = 


.  3E 

n  _i _ z 

c  2  9r 

K 


in  terms  of 


(4-10) 


where  k  is  defined  by 

k 2  (to)  5  -  (A£  +  k2)  .  (4-11) 

From  the  definition  of  n.  Equation  2-2,  this  becomes 
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(4-12) 


K=  £  (f  -  jZ0<»  “k2]1/2  * 

=  [  (e-cos^x)  (~)2  -  Jz0a  •  (4-13) 

In  the  second  line  k  has  been  evaluated  by  Equation  4-4.  The  choice  of 
sign  of  the  square  root  here  will  always  be  taken  such  that  K  would  be 
real  and  positive  if  a  vanished.  Thus  the  imaginary  part  of  <  is 
negative  for  oj  >  0. 


Using  Equation  4-10  in  Equation  4-8  gives  a  differential  equation 


for  Ez, 


3r  3r 


k2E  . 
z 


4-14) 


The  solutions  of  this  equation  are  Bessel  functions  of  order  zero,  Jq  (Kr)  , 
Yg(Kr)  .  (See  Reference  4-1.)  The  value  of  k  is  different  in  insulator 
and  soil.  Let  apply  to  the  insulator,  to  the  soil.  In  the  soil, 
the  field  E^  produced  by  the  wire  current  should  approach  zero  at  large 
r.  This  requires  the  combination  of  J0  and  Y  , 


I'z  =  MJ0(|<2r)  _  jYo(K2rj]  =  Ali0(K2r)  ’ 


(4-15) 


where  A  is  an  arbitrary  constant.  Equation  4-10  then  gives  Bg  in  the 
soil , 


Be  ■  zq  -  jW>>  ;  sj'W’  ■ 


(4-16) 


In  the  insulator,  the  solution  contains  two  arbitrary  constants  B  and 


Hz  =  BJ0(Klr)  +  CY0(Klr)  , 


(4-17) 


B0  =  &T  l BJ  i  (<1^)  +  CYj  (Kjr)  ] 


(4-18) 
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The  constants  A,  B  and  C  arc  determined  by  the  requirements  that  Ei. 


cancel  the  applied  field  at  r  =  a.  and  that  H  and  B„  both  be 

11  0  1  z  0 


continuous  at  r  =  a. 


if  the  two  continuity  conditions  are  written  out  first,  they  can 
be  solved  for  the  ratios  B/A  and  C/A,  with  the  results: 


B  Vl  -  G,,1YC, 


C 

A 


‘Vl  -  (:"lJ0 


(4-19) 


In  these  expressions,  !1(1  and  11^  are  evaluated  at  K2a2’  anc‘  G0’  ''o’  G1 


and  are  evaluated  at  k ^ a 0 .  The  factor  G  is 


G  = 


Kln2 


K2n, 


(4-20) 


and  IV  is  the  Wronskiun 


w  2  Jovi  -  Vo 


Va2 


(4-21  ) 


The  condition  that  I!  cancel  Ih  at  r  =  a  can  be  written  as 

z  0 


li 


0  B  ,  ,  .  C 

A  A  '  oW  +  A  \)  K1  ;*1  • 


(4-22) 


Since  B/A  and  C/A  are  known,  this  equation  determines  A,  and  liquation 
4-19  then  determines  B  and  C,  which  completes  the  solution. 


The  total  current  in  the  wire  is  related  to  the  value  of  B 


0 


at  r  =  a^  by  liquation  1-19,  which  becomes 


2  n  a  n , 


Vi 


l  WP  ♦  s  VVi’ 


(4-23) 


The  impedance  Z  of  the  wire  is  defined  as 


I  (a))  li,, ('» ) / 1  (ui)  ohms/meter  . 
w  0 


(4-24) 


By  use  of  the  result  obtained  above,  2 


cun  be  found  to  be 


I'  I  i  HhESssggiigagggTOjBfia 


z 

w 


“0K1 
2  ira^ 


(VrG1IrVYo  -  (11oYrG11iYo)'Io 


'(11oJrG1!iJ(,)Yi  +  (lloYrGlliYo1Ji 


(4-25 J 


Hero  11^  and  11^ 
and  the  asterisks 
The  functions  II 


arc  evaluated  at  K0a?;  .J^,  .)  ,  and  at  KjU.,; 

indicate  evaluation  at  K^a^  ,  not  complex  conjugation, 
and  llj  are  defined  in  liquations  4-15  and  4-16. 


4.3  SMALL  RADIUS  APPROXIMATION  FOR  WIRE  IMPEDANCE 


The  exact  expression  for  the  wire  impedance,  liquation  4-25,  is 
difficult  to  deal  with,  although  it  can  be  evaluated  numerically.  Reference 
4-1  gives  quite  accurate  polynomial  approximations  for  the  various  Res.sel 
functions  useful  for  numerical  work.  However,  it  is  fortunate  that  the 
arguments  of  all  of  the  Bessel  functions  are  small  in  cases  of  interest  in 

p 

this  report.  The  largest  argument  is  K^a.,.  l-'or  <o  =  1(V  ,  c  r„  10,  c  0.04 

(see  Figures  2-2  and  2-5),  liquation  4-13  gives  | k0 |  »  11  m  .  If  the 

radius  a0  =  0.01  m,  then  lK-,a0|  w  0.11.  For  to  =  10  ,  |  k  0 1  m  0.013,  and 

for  smaller  to,  the  argument  is  even  smaller.  1'hus  it  is  reasonable  to  expect 

that  the  smul 1 -argument  expansions  of  the  Bessel  functions  may  yield  a  suf¬ 
ficiently  accurate  evaluation  of  2  .  These  expansions  are 

'  wr  1 


•»0(x)  =  1  -  +  Oix4)  ,  .)  j  (x)  =  * 

V()(x)  -  I  V.n  ^  (1  -  V>  ♦  V+  O(x') 


0(x5) 


(4-26) 


Y)(x, 


0(  x' 


where  y  is  the  F.uler-Maseheroni  constant, 

Y  =  1.781  ...  .  (4-27 

By  using  these  expansions  and  taking  care  to  collect  till  of  the  terms  of 
lowest  order,  liquation  4-25  can  be  reduced  to 


jSj;,  y&t**  *+****«■**•  *"•  ■ 


77 


Z  (to)  & 


2tt  n 


2 

K2 

—  (An 


'YK2a2 


- 

n  J 


K1  a2 
—  An  — 


(4-28) 


This  formula  holds  for  any  to,  real  or  complex,  provided  |w| 
is  not  so  large  that  the  arguments  of  the  Bessel  functions  are  no  longer 
small.  We  shall  test  the  accuracy  of  the  formula  ir.  the  Laplace  domain 
rather  than  the  Fourier  domain,  i.e.,  for  exponentially  rising  applied 
fields  of  the  form 


Ej  (a)  =  Eg  (a)  exp |at  -  TTz ]  . 


(4-29) 


The  assumption  of  constant  amplitude  and  phase  speed,  Equation  4-2,  now 
implies 


c-  a  a 
k  =  -  =  —  cosy  . 

VC 


(4-30) 


Comparison  of  Equations  4-29  and  4-30  with  Equations  4-8  and  4-4  respcctivel 
shows  that  the  formulae  can  be  written  in  the  Laplace  domain  by  making  the 


replacements 


j«  .  j 


4-31) 


Then  the  media  parameters  become 


0  =  Z()o  (a)  +  --  l(u)  , 


(4-32,: 


(4-33) 


Thus  the  replacement 


(3-34) 


should  also  be  made  in  the  formulae,  liquations  4-12  and  4-13  become 
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K 


(4-35) 


■  [f  f  *  V>  -  ■ 

=  [(e-cos2X)(|)2  +  Z Qc  |]1/2  .  (4-36) 


The  exact  formu’a  for  1^,  Equation  4-25,  could  now  be  rewritten  in  terms 
of  the  nrdifieci  ..^el  functions  (Bessel  functions  of  imaginary  argument)  , 
hut  ’*'e  shall  leave  it  as  it  stands.  It  is  convenient,  however,  to  evaluate 
the  small  radius  approximation  in  terms  of  the  Laplace  domain  parameters. 
Making  the  replacements  indicated  above  in  Equation  4-28  leads  to  the  totally 
real  result 


Z  (a) 
w 


2lr  ln2 


£n 


2 


YK2a2 


(4-37) 


To  test  the  accuracy  of  the  sinai  1  radius  approximation,  we  have 
numerically  evaluated  both  Equations  4-25  and  4-37  for  the  following  cable 
and  soil  parameters: 

a  =  0.5  cm  ,  a  =  1.0  cm  ,  cosy  =  1  , 

1  “  (4-38) 

£  1  =  2  ,  CJj  =  0  ,  e?  =  10  ,  c?2  =  0.01  mho/m  . 

The  computed  results  are  shown  in  Table  4-1.  The  accuracy  is  quite  adequate 
up  to  the  highest  frequencies  of  interest. 


Tabl e  4-1 . 

Comparison  of 

exact  formula  and 

smal  1 

radius  approxi- 

mat ion  for  Z  ( 

a),  for  perfectly 

conducting  wire. 

a 

Exact 

Approximation 

Error 

Skin  Depth 

108  sec"1 

89.21  ohms/m 

89.17  ohms/m 

0.04  % 

0.74  m 

3  x  1  08 

213.1 

212.6 

0.2 

0.31 

109 

506.7 

498.8 

1  .6 

0.16 
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For  '>  *c  cal  application  of  our  results,  note  that  an  BMP  field 

8-1  c 

that  rises  with  i  =  10  sec  to  a  peak  of  about  10  V/m  would  drive  a 

3 

current  of  about  10  A  in  this  wire  at  the  peak  of  the  LIMP .  The  current 
will  continue  to  rise  after  the  peak,  for  the  impedance  is  mostly  inductive, 
as  will  be  seen  in  later  sections  of  this  report.  Note  that  wire  resistance 
is  negligible  compared  with  the  computed  impedance,  at  these  frequencies. 

In  use  of  the  Laplace  domain  results,  it  should  be  noted  that 
a  (a)  and  c  (a)  in  Equation  4-32  are  not  the  same  as  c(oo)  and  e  (w) , 
unless  o  and  a  are  independent  of  frequency.  Equation  2-21  gives 
H(a)  directly  for  universal  soils,  and  Figure  2-5  gives  n(aj  for  our 
standard  soil  (10  percent  water  contentj. 

It  should  also  be  noted  that  the  proximity  of  the  conducting  air 
to  the  cable  has  not  yet  been  included  in  the  analysis.  The  effect  of  the 
air  pro.  mity  is  not  large  for  a  .v  10  sec  and  burial  depths  of  1  meter 
or  more.  After  the  peak,  the  air  proximity  will  have  a  larger  effect, 
enhancing  the  current  during  that  period  in  which  the  air  conductivity  is 
larger  than  the  ground  admittance  q/Z()  • 

The  quantities  2/yi<2  In  Equation  4-28  and  2/yk'2  in  liquation 
4-37  arc  the  skin  depths  in  the  soil;  they  give  a  measure  of  the  depth  of 
penetration  of  the  wire-induced  Fields  into  the  scil.  The  skin  depth  is 
complex  in  Fourier  domain,  in  which  case  its  magnitude  indicates  the  depth 
of  penetration.  It  is  real  in  the  Laplace  domain,  and  values  arc  given 
i  n  Table  4-  . 
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4.4 


THE  TRANSFER  FACTORS 


In  addition  to  the  wire  current,  it  is  useful  to  be  able  to  predict 
the  radial  fields  E^  in  the  insulator  and  the  voltage  across  the  insulator. 
Formulae  for  these  quantities  can  be  obtained  from  Equations  4-9  and  the 
fact  that,  when  the  small -argument  approximation  of  the  Bessel  functions 
is  valid, 


_v_ 

2irrc 


(4-39) 


This  is  the  static  approximation.  Equation  1-19,  and  is  valid  near  the 

wire.  The  result  for  E  is,  in  the  Fourier  domain, 

r 


E»  =  2^rH0M  5  T^VW)  *  (4'40) 

w 

and  in  the  Laplace  domain, 

Br<«>  *  m  r  E0<«>  5  I(0)E»(“)  •  <4-4l> 

w 

The  name  transfer  factor  will  be  used  for  the  factors  T  multiplying  F. 
in  these  equations.  They  are  dimensionless.  According  to  Equations  4-4  and 
4-30  for  the  sweeping  EMP,  they  arc 


T  (oi) 

j  (a)/ c  1  cosy  0 

2li  rq  Z  ’ 

w 

(4-42) 

T  (a) 

(a/c)cosx 

2irrri  Z 

w 

(4-43) 

In  an  insulating  material,  o  =  0,  both  of  these  formulae 

simplify  to 

T  = 

Z_ 

cos^  0 

2iirr:  Z 

(4-44  ) 

In  all  of  the  expressions  for  T,  the  appropriate  ri>  and  e  is  to  be 
used,  i.e.,  n Coo)  or  n(o0,  etc. 

Table  4-2  gives  the  values  of  T  for  the  example  (4-38)  and  for 
points  r  in  the  insulator  just  outside  r  -  a^(T^),  and  in  the  soil  just 
outside  r  =  a^CT^).  It  is  seen  that  E^  will  be  much  larger  than  the 
applied  Eq,  and  that  the  transfer  factors  are  larger  at  lower  frequency. 
Indeed,  the  radial  fields  are  in  the  breakdown  range. 


T  from 


The  voltage  V  across  the  insulator  is  obtained  by  integrating 
a^  to  a^ .  Thus 

V  =  TVE0  ,  (4-45) 


where 


T 


v 


-  C0SX  \  „  *2 

2tr  e  Zv  ‘ 


(4-46) 


The  dimensions  of  T  are  meters.  Ty  is  also  given  in  Table  4-2  for 
the  same  example. 


The  field  transfer  factor  T  is  smaller  for  wires  of  larger  radius, 

as  indicated  by  liquations  4-42  to  4-44.  Z  decreases  as  the  radius 

w 

increases,  but  only  logari.thmi.cly . 


Table  4-2.  Transfer  factors  for  the  example  (4-38). 


a 

Ti 

To 

T 

2 

v 

108 

sec""' 

67.2 

3.2 

0.233  meters 

3  x 

108 

28.2 

2.1 

0.098 

109 

11  .8 

1  .1 

0.041 
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4.5 


R,  L  EQUIVALENT  CIRCUIT  OF  THE  WIRE  IMPEDANCE 


During  the  rise  of  the  EMP,  the  Laplace  domain  form  of  the  wire 
impedance  can  be  used  directly  to  calculate  the  wire  current  and  radial 
electric  field.  After  the  peak  in  the  EMP,  this  simple  procedure  is  not 
applicable.  Fourier  transform  techniques  are  applicable:  transform  the 
applied  field  to  the  a)  domain,  use  Zw((o)  to  calculate  1(a)),  and  invert 
to  find  I(t).  While  this  procedure  is  straightforward,  it  is  time 
consuming,  and  the  analyst  tends  to  lose  contact  with  the  numbers  and 
confidence  in  the  results.  This  occurs  especially  since  the  whole  process 
is  usually  left  in  the  hands  of  computational  technicians,  who  may  have 
little  feel  for  electromagnetics. 

An  alternative  approach  is  to  continue  to  use  analytical  techniques 
to  construct  a  simple,  approximate  method  that  can  be  applied  directly  in 
the  time  domain.  In  this  section  an  11, L  circuit  will  be  devised,  which 
has  the  same  impedance  Z ,  to  good  accuracy,  as  the  actual  wire. 

At  low  frequencies  the  impedance  Z  ,  calculated  for  a  perfectly 
conducting  wire,  becomes  so  small  that  the  finite  resistance  of  the  wire  is 
not  negligible.  We  therefore  add  to  Z  the  resistance  It  ^  (ohms/meter) 
of  the  outer  conductor  of  the  cable.  11  may  depend  on  frequency,  as  it 
will  if  the  thickness  d  of  the  outer  conductor  is  more  than  a  skin  depth 
in  it.  Most  commonly  this  will  not  be  the  case  at  those  low  frequencies 
for  which  II  ,,  is  significant  compared  with  Z  .  We  shall  treat  II  ,, 
as  independent,  of  frequency,  and  take  it.  to  be  the  d.c.  resistance  of  the 
outer  conductor. 

A  convenient  equivalent  circuit  would  have  the  form  shown  in  s  1 
Figure  4-3,  for  then  the  current  in  each  branch  could  be  computed  sepanj^^y 
and  the  total  cable  current  obtained  by  adding  the  branch  currents.  Since 
the  admittance  of  the  network  is  most  easily  written,  it  is  desired  to  make 
the  fit  to  the  wire  admittance  Y  , 


H3 


•v 
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Figure  4-3.  Equivalent  R,L  circuit  to  fit  cable 
impedance. 


- - - =  y 

R  .  +  Z  w 
wO  w 


Ml  1 

^  R  +  iwl. 
n=l  n  J  n 


(4-47) 


This  equation  can  also  be  written 


in  g 

"  E  TTJ, 

n=l  J 


i“/6„  ’ 


(4-48) 


where 


g  =  1/R 
^n  n 


3  =  R  /L 
n  n  n 


(4-49) 


The  fit  can  be  made  by  first  choosing  arbitrarily  a  set  of  Bn's,  spaced 
one  decade  apart,  say.  The  g  are  then  determined  by  requiring  liquation 
4-48  to  hold  exactly  at  a  set  of  ai's,  c.g.  ,  a)  =  8^ ,  37,  3^,  •.*3m»  which  gives 
a  set  of  m  linear  equations  to  solve  for  the  m  quantities  g  . 


The  fit  can  be  made  either  in  the  hup  lace  or  Fourier  domain.  In 
the  Laplace  domain,  both  sides  of  liquation  4-48  are  real.  In  the  Fourier 
domain,  one  can  fit,  say,  the  imaginary  (inductive)  part  of  Y  (w) .  The 
real  (resistive)  part  should  then  fit  automatically,  since  both  sides  of 
liquation  4-48  are  analytic  functions.  The  fit  will  not  be  exact,  of  course, 
at  frequencies  between  the  fit  points. 


Lavery  (Reference  4-2)  used  this  technique  to  fit  the  admittance 
of  a  cable  with  parameters 

a.  =  0.5  cm  ,  a  =  1.0  cm  ,  d  =  0.05  cm  ,  cosX  =  1,  ) 

/  (4-50) 

e1  =  2  ,  a  =  0  ,  R  =  1.129  X  10  ohm/m  ,  ) 

imbedded  in  our  standard  soil.  The  3  chosen  were 

n 

1  +  n 

3  =  5  X  10  ,  II  =  1  to  7  .  (4-51) 

The  real  and  imaginary  parts  of  Y  (<n)  and  of  the  fit  are  graphed  in 

2  q  W 

Figure  4-4  for  10  s:  w  £  10‘  sec  ^ .  It  is  seen  that  the  fit  is  quite 

g 

accurate,  except  for  the  real  part  at  u  10  .  This  discrepancy  could  be 

removed  by  adding  higher  (3^ 1  s ,  but  we  shall  not  depend  on  the  fit  at  these 

high  frequencies.  The  values  of  g  ,  R  and  L  for  the  fit  are  given  in 

n  n  n  ^ 

Table  4-3.  In  this  table,  an  entry  a(b)  means  ax  io  .  The  zero 
frequency  conductance  and  the  infinite  frequency  inductance  Lw  of 

the  network  are 

Table  4-3.  Fit  parameters  for  the  example  (4-50)  in 
standard  soi 1 . 


n 

gn (m/ohm) 

Rn  (ohm/m) 

L  (henry /m) 

1 

8.439(2) 

1 .1 85(-3 ) 

2 . 370( -6) 

2 

2.137(1  ) 

4.679(-2) 

9. 359( -6) 

3 

1 .306(0) 

7 . 657  (-1  ) 

1  .531 (-5) 

4 

1  . 972 (-  1  ) 

5.071 (0) 

1  . 01  4  ( -5 ) 

5 

2 . 684 ( -  2 ) 

3.726(1  ) 

7 .452 (-6 ) 

6 

5.002(-3) 

1 .999(2) 

3 . 990( -6) 

7 

1 .258 (-3) 

7.949(2) 

1  . 5  90 ( - 6 ) 
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Figure  4- 


Real  and  imaginary  parts  Ywr  and  Yw-j  of  cable  admittance 
as  a  function  of  w.  The  curves  are  computed  from  the 
definition  of  Yw,  Equation  4-47,  with  Zw  given  by  Equation 
4-28.  The  points  are  computed  from  the  R,L  circuit  fit.  The 
fit  was  made  to  Yw-j  at  the  points  u  =  10n+2,  n  =  1  to  7. 
The  R,L  model  values  for  Yw~  are  also  very  close  to  tne 
true  values,  except  at  m  ^  10  .  This  discrepancy  could  be 
removed  by  adding  R,L  branches  with  higher  rates  p  . 


HG 


(4-52) 


G  =  8  =  8.668  x  lCf'  mho-m  *=»  1/R  ,  \ 

0  *-*•  n  wO  ) 

n  ( 

L  =  ( gn0n)  *  =  5.862  x  io  7  henry/m  . 


The  admittance  in  the  .Laplace  domain  is 


'  6*. 

Y  (a)  <*  T  - -- 

w  1  +  a 


(4-53) 


— ,  ^  ■  a/8 
n=l  a 

Figure  4-5  compares  Y  (a )  with  the  fit  for  the  example  (4-50). 

t  h 

The  current  in  the  nx  branch  of  the  network  can  be  computed  by 
solving  the  differential  equation 

dl 


L  ~  +  R  I  =  E,,(t)  , 
ii  dt  n  n  0  ' 


or 


Tr+  <Vn  " 


The  solution  of  this  equation  is 


-8  (t-t  )  -8  t  r  8  t‘ 

Jntt3  '  0  W  '  B,.fne  J  0  ‘0 

t0 


(t')dt' 


(4-54) 


(4-58) 


This  solution  allows  for  an  arbitrary  initial  current  I^ft^)  at  the 
starting  time  ,  in  case  it  is  convenient  to  approximate  ii^ft)  by 
different  analytical  forms  in  different  time  periods.  For  example,  a 
somewhat  crude  but  useful  approximation  to  the  early-time  part  of  the 
horizontal  electric  field  (see  Figure  4-6)  is 


Ii  ft)  =  ii  eat  for  t  <  0  , 

0  m 

=  Ii  c"Yt  for  t  >  0  . 
m 


(4-56) 
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Here  E  is  the  maximum  value,  achieved  at  t  =  0..  If  E  is  chosen 
m  m 

somewhat  larger  than  the  actual  peak  field,  the  form  (4-56)  will  bound  the 
actual  field,  and  the  computed  current  will  therefore  bound  the  actual 
current. 


For  the  form  (4-56)  ,  the  current  for  t  <  0  can  be  calculated 


directly  from  Y  (a) , 


I(t)  =  Y  (a)E  e  u  for  t  <  0  . 
v  •  w  m 


(4-57) 


For  this  purpose,  either  the  original  Y^Ca)  or  the  fit  to  it  may  be  used. 
At  t  -  0  =  t  ,  there  is  already  a  current  flowing  in  the  wire,  and  the 

v  * 

part  of  this  current  flowing  in  the  n  branch  is 


g  E 

r  ( m  -  n  111 

V03  '  1  +  a/6. 


(4-58) 


From  Equation  4-55,  the  nu'  current  for  t  >  0  is  easily  calculated  to  be 


-St  gSE  -6t 

In(t)  =  e  n  In(0)  +  (e'Yt-e  n  )  . 


(4-59) 


The  apparent  singularity  in  the  second  term  here  if  y  should  approach 

3  actually  does  not  occur,  us  the  exponentials  also  cancel  in  the  limit 

y  -*■  6  .  Summing  I  (t)  over  n  gives  the  result,  for  t  >  0, 
n  n 


g  6  e 

I  (t)  =  (a+Y)Ii  T 


m  ^  (a+3n)(Y-3n)  wl  m 


(4-60) 


Here  Y  (-y)  is  Y  (a)  evaluated  at  a  =  -Y* 
w  w 

The  form  (4-56)  has  a  discontinuity  in  slope  at  its  peak,  which 

gives  it  more  high-frequency  content  than  the  true  liMP.  A  function  with 

a  discontinuity  in  slope  has  a  Fourier  transform  which  falls  off  no  faster 
2 

than  1/w  at  large  w,  whereas  the  transform  of  the  true  EM!’ 


u,.'..  ii  rrl  ■>  taVa u.-ihdJ  iJ  j.JL-ii  ■  i 


(which  has  all  derivatives  continuous)  falls  off  faster  than  any  finite  power 
of  co.  The  form  can  be  improved  easily  by  adding  other  exponentials.  For 
example,  the  form 


c  ...  c  at  2at. 
E0(t)  =  Emf2e  'e  ) 


«  E  (  —  e’Yt 
m\2a-y 


,  t  <  0  , 

Y  -2at 
2a-Y  £ 


t  >  0 


(4-61) 


has  maximum  value  H  at  t  =  0  and  zero  slope  there,  as  is  easily  verified. 

m  r 

This  function  is  graphed  in  Figure  4-6  for  the  case  Y  =  a/2.  The  function 

3 

has  discontinuous  second  derivative,  and  its  Fourier  transform  falls  as  l/oo 

at  large  w.  It  is  possible,  with  the  rise  and  decay  rates  a  and  y 

fixed  in  the  first  exponentials  on  each  line,  to  vary  the  constants  in  the 

last  two  exponentials  and  the  coefficients  of  all  terms  in  such  a  way  that 

t  h 

continuity  is  maintained  through  the  fourth  derivative.  If  the  N  derivative 

N+ 1 

shows  the  first  discontinuity,  the  Fourier  transform  falls  as  l/oo 
asymptotically.  The  algebra  determining  the  coefficients  becomes  quite 
complicated  for  the  very  smooth  forms.  Use  of  the  simple  forms  bounds  the 
high-frequency  content. 


The  current  for  the  form  (4-61)  can  be  written  down  by  applying 
liquations  4-57  and  4-59  to  each  of  the  exponential  terms.  In  this  way, 
all  but  the  final  summing  of  the  terms  can  be  done  analytically.  Alternatively, 
the  differential  liquation  4-54  or  the  integral  in  liquation  4-55  can  be  solved 
or  evaluated  numerically.  Hither  approach  gives  a  fast  and  quite  accurate 
way  of  calculating  the  cable  current. 


Note  that  the  transfer  factor  for  the  insulating  layer,  liquation 
4-44,  displays  no  frequency  dependence  except  in  the  factor  Z  .  Therefore, 
on  including  the  d.c.  wire  resistance,  it  is  possible  to  write 
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.  _  zoc°s*  Ko 

r  2irre  R  .  +  Z  ’ 
wO  w 

ZqCosx 

=  “ 2TTre—  1  ’  (insulator)  (4-62) 

and  the  la  line  here  holds  in  either  frequency  or  time  domains.  Calcula¬ 
tion  of  I  therefore  imediately  yields  E  in  the  insulator. 
Unfortunately,  this  is  not  true  for  11  in  the  soil.  At  lower 
frequencies,  where  n  is  approximately  constant, 


i!  «  i  |> . 

r  2'irroQ  c  it 


(.soil,  near  wire) 


(4-63) 


At  high  frequencies  where  c  rather  than  0  dominates  n  and  c.  does 
not  vary  strongly  with  frequency,  liquation  4-62  applies  approximately  in  the 
soil  near  the  wire. 


4.6  THE  TIME-VARYING  INDUCTANCE  MODEL 


The  method  developed  in  Section  4.5,  while  transparent,  and  accurate, 
still  requires  a  considerable  amount  of  calculation.  A  simpler  method  is 
desirable,  even  »f  it  is  loss  accurate.  Such  a  method  exists  for  the  type 
of  applied  fields  found  in  liMP  environments. 


The  expression  (4-28)  for 


can  be  simplified  by  milking  some 


approximations,  liquation  4  11  shows  that,  if  the  k"  term  is  dropped 


k  .  U) 

n  r’"  ■' 


c 


(4-64  ) 


Mow,  the  k'  term  contributes  the  term-cos"/  in  liquation  4-13.  In  the 
2 

soil,  cos  ,<  is  no  more  than  about  .10  percent  of  r,  whore  t  •  10.  In  the 
insulator,  dronping  co';  'y  makes  a  bigger  percentage  change.  The  logarithm 
factor  tf.n(a,/,'i|)  is  usually  much  smaller  than  the  logaiithm  for  the  soil 


lldl'yiiailiiiiH  ill  *  I  i’1 
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term,  so  that  the  insulator  term  is  not  very  important.  However,  if  a  is 
.  .  2 

negligible  in  the  insulator,  k  /n  is  still  simple,  and  can  be  retained. 
The  approximate  form  of  Z  is 


Z 

w 


2'fTC 


2 

cos  X 


•jfi-n 


(4-65) 


(Note  that  ftn(-j)  =  -jfr/2.)  This  expression  gives  an  impedance  which  is 
too  large,  but  not  by  more  than  about  10  percent. 


If  the  logarithm  in  Equation  4-65  were  independent  of  frequency 
and  real,  it  would  represent  the  impedance  of  a  pure  inductance.  If 
is  written  in  terms  of  its  magnitude  and  phase, 


(4-66) 


it  is  clear  from  Equation  4-13  that  tj»  varies  between  zero  and  it/ 4 ; 

<|)  approaches  zero  when  the  dielectric  term  is  dominant  and  approaches 
ir/4  when  the  conductivity  term  is  dominant.  On  separating  real  and 
imaginary  parts,  Z becomes 


Z 

w 


j"zo 

2irc 


+ 


2 

cos  X 
c 


1 


2 'ire 


(4-67) 


The  absolute  value  on  the  factor  w  in  the  real  term  is  demanded,  for 
negative  u>,  by  the  reality  condition,  Z*(to)  =  Zw(-u>),  and  can  also  of  course 
be  derived  from  the  properties  of  the  Bessel  functions  by  carrying  through 
the  analysis  for  negative  w.  The  second  term  is  resistive  and  leads  to 
energy  dissipation.  If  the  soil  is  a  good  dielectric,  cp  -+  0,  and  the 
energy  is  radiated  away.  If  the  soil  is  a  good  conductor,  <|>  rn/4,  and 

the  energy  is  dissipated  in  Joule  heating  of  the  soil. 
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The  magnitude  of  the  resistive  term  is  generally  fairly  small 


compared  with  the  first  or  inductive  term.  The  first  logarithm  is  typically 


between  3  and  12  in  cases  of  interest,  whereas  (—  -  <J>)  ranges  between 


ir/2  and  tt/4.  Dropping  the  resistive  term  should  increase  the  calculated 
current  and  provide  an  upper  bound.  We  shall  drop  it,  and  compare  calculated 
results  with  those  obtained  from  the  more  accurate  equivalent  circuit  method. 


In  terms  of  the  inductance  L  of  the  wire, 


jtoL  , 


(4-68) 


and 


27TC 


£n 


/  2  \ 
W1  kYI  a  2  7 


1  - 


2 

cos  x 


a2 
Jin  ~ 


(  -69) 


Note  th  t 


^0  -7 

=  2  x  10  Henry/meter  . 


2trc  2tt 


(4-70) 


The  inductance  depends  on  l^l'  hut  only  quite  slowly  because  of  the 


logarithm.  Thus  an  approximate  fit  to 
to  Equation  4-64, 


would  be  adequate.  According 


cT  'V 


(4-71) 


A  good  fit  to  VhJ  is  given  by 


V'n2  1  »  V 

1  10 

e  —  + 

oo  £ 

(4-72) 

For  our  standard  soil,  the 

yield  the  results: 

exact 

evaluation  of  ^ |  q9  | 

and 

this  approximation 

to 

It 

3 

104 

io5  io6  io7 

io8 

109  ) 

exact  =  1.76 

1.77 

1 . 79  1.82  2 . 04 

2.90 

( 

5.94  {  (4-73) 

approximate  =  1.74 

1  .73 

1.78  1.87  2.15 

3.03 

5.82  ' 
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exact  result  over  the  entire 


The  approximation  is  within  5  perce.  ..  ^ 
range.  Therefore,  a  good  fit  to  2/‘» 


6 


yT^T  >.  (v^;  *  vi^) 


(4-74) 


The  quantity  6  defined  here  is  the  skin  depth  in  the  soil  for  the  present 
cylindrical  problem.  In  terms  of  <5,  the  inductance  is 


L 


(4-75) 


Because  of  the  logarithm,  the  inductance  changes  only  slowly  with 


frequency.  For 

the 

example 

(4-50), 

Equation  4-75 

gives  the 

results : 

w  =  103 

104 

q 

10 

io6 

107 

108 

109 

sec”1  ) 

6  =  353 

111 

34.5 

10.4 

2.86 

0.641 

0.105 

( 

meters / (4  -  76 ) 

L  =  2.16 

1.93 

1  .70 

i.46 

1.20 

0.90 

0.54 

pH/ m.  ' 

The  relation  between  the  applied  Eq(u>)  and  the  current  1(a))  , 
including  now  the  wire  resistance  ^wq>  can  be  written 

jojfLo')1^]  1  00(oj)  -  RwQI  .  (4-77) 

The  exceedingly  slow  variation  of  L(w)  suggests  that  it  might,  be  a 
reasorably  good  approx  .mation  to  regard  L(w)  as  constant  in  inverting  this 
equation  to  the  time  domain,  it  would  appear  appropriate  to  use  the  value 
of  I.  for  that  range  of  (o  which  gives  the  dominant  contribution  to  the 
integral  op  I;  (oj) exp ( jut)  over  w,  i.e.,  to  the  inverse  Fourier  transform 
of  liquation  4-77.  Since  the  exponential  here  oscillates  rapidly  with  w 
when  tut  »  I,  the  appropriate  range  of  w  is 

uJ  «  1/t  ,  (4-78) 

unless  li  (w)  varies  rapidly  with  u>.  The  latter  possibility  depends  on 
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the  choice  of  the  time  origin.  The  F.MP  E  (t)  varies  rapidly  with  t 
initially,  and  then  more  and  more  slowly  at  later  and  later  times.  If  the 
time  origin  is  chosen  at  that  time  when  E^ft)  ps  changing  most  rapidly, 
then  Eq(w)  will  not  vary  rapidly  with  oo.  For  example,  the  form  (4-56) 
for  E 0  C t )  has  the  Fourier  transform 


nn(w)  -  li 


m  (Y+jw) (a-jw) 


(4-79) 


However,  if  the  time  origin  were  shifted  to  time  t(),  then  E ^ (to)  would 
acquire  a  multiplicative  factor  exp(-jwt0).  Thus,  rapid  variation  is 
avoided  by  choosing  the  time  origin  as  stated.  We  agree  to  make  this 
choice,  and  use  liquation  4-78  (otherwise  we  would  choose  ui  &  l/(t-t  )). 


The  suggested  approximate  time  domain  equation  is  therefore 

[L(t)I  (t) ]  =  R0(t)  -  Rw0l  ,  (4-80) 

where  Lft)  is  given  by  liquation  4-75  with  6  evaluated  from  liquation 
4-74  and  4-78,  L .c.  , 


1 . 1 2ct 


(4-81) 


This  prescription  would  have  difficulty  if  used  at  t  =  0.  However,  the 
EMI’  lias  a  finite  rise  time  t  ,  which  for  the  exponential  rise  model  i.s 
t  1/a.  Thus  in  using  liquation  4-81,  t  should  be  set  equal  to 

t  =  tr  =  1/a  if  t  <  t  ,  (4-82) 


and  the  time  origin  should  be  chosen  so  that  t  =  t  when  the  rise  rate 
begins  to  fall  significantly  below  a.  Obviously,  this  prescription  is 
somewhat  imprecise  at  early  times.  Note,  however  that  quite  accurate 
evaluation  of  the  current  during  the  rise  of  the  applied  field  can  be  had 
by  using  the  form  (4-56)  or,  better,  the  form  (4-61)  and  the  Laplace  domain 
impedance.  For  the  latter  form,  the  current  is 


1 


,  .  r  / _ 2  at  1 

^  ni  \  Z  (a)  e  Z  (2a) 
'  w  w 


0  .  (4-83) 


Lavery  (Reference  4-2)  has  tested  the  approximate  equations  against 
the  equivalent  circuit  model  cf  Section  4.5  for  the  applied  field 


EQ(t)  =  0  for  t  <  0 


(4-84) 


=  Ea( 


-Yt  -at.  . 
e  -e  )  tor  t 


for  which  the  equivalent  circuit  equations  can  be  solved  analytically.  The 

8  “1 

cable  example  (4-50)  was  used,  with  a  =  10  sec  .  In  Equation  4-81,  t 

-  8 

was  set  equal  to  the  larger  of  t  and  10  seconds.  Equation  4-80  was 
integrated  numerically  for  several  values  of  y  as  indicated  in  Figure  4-7, 
which  compares  the  currents  computed  by  the  two  methods.  The  comparison 
shows  that  the  varying  inductance  approximation  i.s  quite  good. 


Although  Equation  4-80  was  integrated  numerically  in  this  example 
to  show  how  good  the  approximation  is,  the  current  can  obviously  be  estimated 
by  crude  integration  of  the  equation  over  blocks  of  time  in  each  of  which 
li^(t)  and  L(t)  are  regarded  as  constant.  Note  that  the  decay  time  L/R^ 
is  of  the  order  milliseconds. 


Our  choice  of  putting  L(t)  inside  the  time  derivative  in  Equation 
4-80  instead  of  outside  was  somewhat  arbitrary,  mathematically.  The  choice 
made  gives  better  agreement  with  the  accurately  calculated  currents  for 
the  shorter  driving  pulses.  If  b  wore  put  outside,  the  current  would  not 
decrease  immediately  after  the  short  pulse,  but  would  decay  only  on  the  time 

„  Z  ... 

scale  L/R  ,,  rj  1 0  ‘  second.  Putting  L  under  the  time  derivative  gives  hack 
some  of  the  dissipation  associated  with  the  real  part  of  Y  (to)  ,  through  the 
term  Idh/dt.  Since  db/dt  is  positive  this  term  has  the  effect  of  a 
res i stance . 
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Comparison  of  currents  calculated  from  R,L  equivalent  circuit 
(curves)  and  from  time-varying  inductance  model  (points), 
for  the  applied  field  of  Equation  4-84.  The  values  of  y  are 
i ndi cated . 
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4.7 


FREELY  PROPAGATING  SOLUTIONS 


The  analysis  thus  far  has  considered  only  the  particular  solution 
of  Maxwell's  equations,  or  that  part  driven  by  the  applied  field.  The 
equations  also  allow  freely  propagating  solutions,  and  these  are  needed  to 
satisfy  conditions  at  the  end  of  the  cable.  The  freely  propagating  solutions 
~  exp(j  (u)t-kz))  occur  for  w  and  k  such  that 

Zw  -  0  ,  (4-85) 

so  that  a  current  I  car.  exist  without  an  applied  field.  Thus  k  is  no 
longer  determined  by  liquation  1-4,  but  for  each  it  is  chosen  to  make  Z  ^ 
vanish . 


Z  is  given  by  Equation  4-28,  and  Equation  4-11  gives 

2  o 

K  .  00  k 

_  - _  4  —  +  ~ 

n  c  n 


(4-86) 


for  arbitrary  k  and  u).  Ignoring  the  dependence  of  the  logarithm  on 
k  (through  <2)  allows  solution  of  Equation  4-85  for  kz,  with  the  result 


,2  .0) 

k  -  -  J  -  H; 


!bn ((6/ £n(a2/a^) 


n 


(4-87) 


1 


lln(a2/a1)  +  2,n(6/a,,) 


where  the  skin  depth  6'  for  propagation  is  complex, 


d  =  2/yjK2  .  (4-88) 

2 

Since  k  and  6  depend  on  k  (sec  Equation  4-12),  Equation  4-87  has  to 
be  solved  by  iteration  for  accurate  results,  in  general. 


for  u)  <  10^  sec  ^  and  when  the  conductivity  of  the  insulator 
is  indeed  small  (true  except  in  case  of  high  radiation  exposure),  the  first 
step  in  the  iteration  gives  approximately  correct  results.  The  first  step 
puts  k“  =  0  in  the  calculation  of  and  6.  for  to  <  i()^  the  dielcctri 

L.  ^ 

part  of  p  curi  be  neglected  and 
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V-j^0O0a)/c  * 


so  that 


Ihi 


V) 


6  V- j 
m 


(4-89) 


(4-90) 


Y  ',"oao“/c 


Here 


is  the  magnitude  of  .6.  '’  Thus 

\ 

'6 


An  I 


=  An 


(a  -  j 


(4-91) 


The  ratio  of  the  i] 1  s  is 

j  G  Jjj/  C 

n2  Z0°Q 


for  the  example  (4-50)  at  u>  =  10(\  ^n(6  /a  )  «7.3,  Ania^/a^  =  0.695, 
and  n, /n0  ~  j/450.  The  second  term  in  the  denominator  of  liquation  4-87 
can  he  neglected,  with  the  result. 

2  An  (6  /a  )  -  i  •; 

,2  _  /  u>  m  1J  -4 


An  (a7/ »j) 


(4-91) 


This  result  verifies  that  it  was  proper  to  neglect,  k"  in  calculating  kv 
Since  An(6  /a^)  is  considerably  larger  ihan  ir/4,  the  square  root  can  he 
calculated  approximately,  giving  the  formula  for  k, 


H6,  i/'V 


i/r/8 


k  "  ’  C  yjl n(a2/ai)'  V  '  AiU^/a  ")/  ‘ 

The  imaginary  part  of  k  gives  the  attenuation  of  the  propagating  wave. 


■> 

Tor  large  ui,  the  solution  for  k",  liquation  4-87,  should  ho 

7 

iterated  by  putting  this  value  of  k“  Sack  into  the  formula  (liquation  4-86) 
for  Kq ,  and  so  on.  Ter  accuracy,  the  full  n  should  of  course  be  used, 
instead  oi  Z  0^.  figure  4--8  shows  the  real  and  imaginary  parts  of  k  for 
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Figure  4-8.  Cable  propagation  constants.  Solid  curves,  real  and  imagi 
parts  kr  and  k-j  as  functions  of  m.  Short  dashed  curves,  t 
same  from  the  approximate  Equation  4-92.  Points,  K(i)\  nc 
that  this  is  virtually  the  same  as  kr(io).  Long  dashed  cur 
ratio  of  phase  speed  to  speed  of  light. 


the  example  (4-50)  ,  along  with  results  from  the  approximate  formula  Equation 
4-92. 


(4-99) 


tZ0  (a2\ 

v'a>  ■  ■ 

Note  that  the  sign  of  V  depends  on  the  sign  of  k  or  Tc,  i.e.,  it  is 
different  for  solutions  propagating  to  right  and  left.  The  factors 
multiplying  1  in  these  equations  are  analogous  to  the  characteristic 
impedance  of  a  coaxial  transmission  line.  In  the  present  case,  the 

outer  conductor  (the  soil)  does  not  have  perfect  conductivity,  and  there 
are  significant  electric  fields  in  it  out  to  radii  of  the  order  of  6. 

4.8  THE  END  CONDITIONS:  OPEN  CIRCUIT 

The  freely  propagating  solutions  are  added  to  the  driven  solution 
to  satisfy  conditions  at  the  end  of  the  wire.  The  easiest  case  to  analyze 
is  that  in  which  the  end  of  the  wire  is  insulated  from  the  soil,  i.e.,  the 
case  of  open  circuit.  If  the  driven  and  freely  propagating  currents  are 
denoted  by  I ^  and  I  respectively,  the  end  condition  is  then  that 

I  -  -  I  ,  .  (4-100) 

P  d 

The  driven  current  is  determined  by  the  applied  electric  field,  and  propagating 
currents  are  then  led  into  the  ends  to  cancel  the  driven  current  at  thos'- 
points.  The  propagating  currents  propagate  along  the  wire,  modifying  the 
total  current  and  fields  as  they  go.  They  eventually  reach  the  opposite 
end  from  their  origin.  At  that  time,  additional  propagating  currents  are 
fed  into  the  ends  to  cancel  the  outgoing  propagating  currents.  The 
analysis  here  is  the  same  as  in  normal  transmission  line  analysis. 

This  analysis  is  not  exact.  While  the  propagating  solution  can 
cancel  the  wire  current  of  the  driven  solution,  it  does  not  cancel  the  fields 
of  the  latter  in  detail.  The  values  of  k,  and  therefore  of  r.,,  are 
different  for  the  driven  and  propagating  solutions  for  a  given  frequency. 

Hence  the  radial  distribution  of  fields  is  not  the  same.  Cancellation  of 
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the  currents  means  that  the  fields,  which  are  given  for  both  solutions 

by  liquation  4  39  near  t ho  wire  (out  to  about  one  skin  depth),  will  also 
cancel  approximately.  However,  as  will  be  seen  shortly,  the  fields  li  are 
substantially  different  for  the  driven  and  propagated  solutions,  and  comic 
nowhere  near  canceling.  It  is  reasonable  to  expect  an  adjustment  of  the 
radial  electric  field  over  a  region  of  the  size  of  one  skin  depth  near  the 
end  of  the  wire. 

liquation  4-9  relates  11  ^  to  .  Since  Bq  is  approximately  the 

same  for  driven  and  propagated  solutions,  given  liquation  4-100,  out  to  about 

one  skin  depth,  the  difference  in  li  comes  from  the  difference  in  k.  For 

the  driven  solution,  is  given  by  l-.quution  4-4.  For  the  propagating 

solutions,  k  is  given  by  liquation  4-87,  or  for  lower  frequencies  by  the 

approximate  liquation  4-92.  It  is  seen  that  k  is  typically  several  times 

k,.  Thus  li  is  several  times  b  ,. 
d  rp  ru 

The  adjustment  of  I;  involves  propagating  solutions  with  higher 
radial  modes.  The  k  for  these  modes  is  again  found  by  requiring  I  -  0, 
but  this  time  the  exact  liquation  4-25  must  he  used  instead  of  the  approximate 
liquation  4-28  resulting  from  the  small  argument  expansion.  These  solutions 
decay  rapidly  with  distance  away  from  the  wire  end,  and  carry  little  current 
so  that  liquation  4-100  remains  approximately  correct.  Their  role  is, 
roughly,  to  remove  li  from  the  soil  and  increase  li^  in  the  insulator 
near  the  end  of  the  wire.  They  arc  important,  therefore,  in  considerations 
of  insulator  breakdown  at  the  wire  end.  The  total  solution  depends  also 
on  the  structure  of  wire  and  insulation  at  the  end.  The  field  across  the 
insulation  can  be  reduced  by  connecting  the  end  of  the  wire  to  a  larger 
conducting  sphere  which  is  also  insulated  from  the  soil.  \'o  detailed 
soluions  were  available  at  the  time  of  writing  of  this  report. 

The  radial  field  in  the  insulator  is  at  least  as  large  as  the 
values  indicated  by  liquations  1-9<>  or  4-97.  With  the  approximation  liquation 
4-92  for  k,  these  equation:  become,  in  the  insulator, 
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r  2fTr  Vc' 


An (6  /a. ) 


An(a  /a  )  p 


C  l- 3  oi  ) 


A  measure  of  the  importance  of  the  fields  in  the  soil  to  those  in 


the  insulator  is  given  by  the  ratio  of  the  radial  voltage  drop  V  =  f\:.  dr 


in  the  two  regions.  On  taking  Bfi  ~  1/r,  liquation  4-9  leads  to 


n1  An(5/a2) 


P2  An (a^/a^) 


14-1025 


This  ratio  was  shown  to  be  small  for  oj  ^  10  in  Section  4.7.  It  is  graphed 
as  a  function  of  Laplace  domain  frequency  in  figure  4-9  for  our  usual  example 
which  shows  that  it  is  always  considerably  less  than  unity.  It  is  to  be 


expected  that  the  increase  in  0^  in  the  insulator  due  to  adjustment  of 


fields  near  the  wire  end  is  by  a  factor  1  +  (V  .  ,/V.  ),  which  is  not  large 

J  soi 1  ins 


The  phase  speed  of  propagation  along  the  cable  is  fk  is  the  real 


part  of  k) 


i>)  a 

v  -  — -  or  -  . 

*  kr  r 
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This  speed  is  considerably  less  than  c,  as  shown  by  liquation  1-92.  The 


ratio  v./c  is  graphed  in  figure  4-8. 

y 


The  maximum  value  of  the  current  at  the  center  of  the  cable  will 


usually  occur,  for  11MF  drive,  just  before  the  propagated  signals  arrive 
from  the  two  ends.  Attenuation  and  dispersion  of  the  propagated  signals 
must  he  taken  into  account.  An  approximate  way  of  treating  this  is  developed 
in  the  following  section. 


Fi 'jure  4-9.  Ratio  of  the  radial  voltage  drops  in  the  soil 

and  in  the  insulator  for  the  propagating  solution 
in  the  Laplace  domain,  as  a  function  of  frequency 
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ATTENUATION  AND  DISPERSION 


The  factor  exp(-jkz)  determines  how  the  signal  propagates  along 
the  cable.  In  order  to  estimate  the  effects  of  attenuation  and  dispersion, 
an  analytic  approximation  to  k(w)  is  needed.  Figure  4-8  shows  that  the 
real  part  k^  can  be  fitted  quite  well,  over  a  few  decades  of  oj  about 
any  point,  by  a  power  law.  In  choosing  a  fit,  care  must  be  taken  to  maintain 
the  reality  and  causality  conditions  (see  Section  2.3).  A  satisfactory 
approximation  is 

jk  =  a(joj)P  ,  (4-104) 

where  the  power  p  and  the  factor  a  are  both  real  constants.  For  then 
the  complex  conjugate  of  jk  is 

[jk(w)]*  =  jk (-oo)  ,  (4-105) 

and  reality  is  guaranteed.  Causality  requires  that  jk  be  analytic  in  the 
negative  imaginary  half  of  the  complex  ui  plane.  The  function  4-104  is 
analytic  in  any  region  that  does  not  enclose  the  origin.  A  cut  along  the 
imaginary  axis  from  w  =  0  to  w  --  j°°  prevents  encircling  the  origin. 

The  path  of  integration  in  any  Fourier  inversions  must  then  pass  below  the 
origin,  as  in  Figure  4-10a. 

Figure  4-8  shows  that  the  exponent  p  is  just  slightly  less  than 
unity.  Thus  q,  defined  as 

q  ~  1  -  p  ,  (4-100) 

is  a  small  number,  of  the  order  0.05  to  0.1. 

liquation  4-104  can  also  be  written,  for  w  on  the  positive  real 

axis , 


iL  <  AS'Jlil  A  C '!>• J.  3*- 1-1  J_  1 
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Figure  4-10.  Cut  and  integration  path  in  complex 
oj  plane. 

k  =  aj-V  -  a^expC-j  fl 

=  | cos (^-)  -  j sinful  |  (4-107) 

«  ao)l’|  1  -  j  ^j'j  (for  q  <<  1)  . 

This  sliows  the  relation  of  the  real  and  imaginary  parts  of  k.  The  closer 
p  is  to  unity,  the  smaller  is  the  ratio  of  k.  to  k ^ ,  which  explains  the 
shape  of  the  k.  curve  in  Figure  4-8.  If  we  choose  to  fit  the  k  curve 
in  the  vicinity  of  some  frequency  u>(),  where  k  =  k() ,  then  the  constant  a 
i  s 

a  «  k  /wjj  .  (4-108) 

Let  an  impulse  current  he  injected  at  t  =  l)  into  the  end  of  the 
cable,  which  is  at  z  =  0.  The  Fourier  transform  of  the  impulse  function 
is  unity.  The  signal  at  t,z  is  then 


i m.  ■  •  >  .rii< 


(4-109) 


Iim(t.z)  =  /  exp [jut  -  az(ju)‘]du  . 

.OO 


The  integral  here  id  difficult  to  evaluate  exactly,  but  it  can  be  estimated 
by  the  method  of  stationary  phase  (or  saddle  point).  The  argument  of  the 
exponential  is  stationary  when 

3(jw)  arS  =  t  -  paz(ju)"q  =  0  ,  (4-110) 

or  at 


Jw  =  jU)st 


(4-111) 


This  point  is  on  the  negative  imaginary  axis  (oj  =  -j  x  real  number),  so  that 
the  integration  contour  can  be  deformed  to  pass  through  it.  The  value  of 
the  argument  at  the  stationary  point  is 


arg0  -  t(*f) 


i/q 


az(^f) 


p/q 


L  I  •  ^ 


( 1 -q) / q 


(4-11. 


The  second  derivative  of  the  argument  is 


0“  0  ,-(l+q)  ,, 

“2  11  rh  =  "  - 7  arg  =  -  qpac(ju))  '  (4-1  lo) 

Din  cl(jw)“ 


and  evaluation  at  the  stationary  point  gives 

t  (l+qi/q 

arKo  =  ‘  qPazfe'J 


(4-114) 


In  the  vicinity  of  the  stationary  point,  u>  =  u  +  din 


arg  »  a  rg()  + 


\  arg'J 


(4-115) 
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Thus  the  exponential  in  Equation  4-109  is  real  and  decaying  if  the  path  of 
integration  passes  horizontally  (6co  real)  through  the  stationary  point. 


The  estimate  of 

I. 

im 

is  then 

I.  (t,z) 
nn  *  1 

1 

rv  - 

J~  (■ 

>  qpaz  v 

rO 

2TT 

(l+q)/2q  (l-q)/q 


expt-qazC2—) 


]  .  (4-116) 


Since  (l-q)/q  is  a  large  number  («  10  to  20),  the  exponential 
here  makes  I  very  small  until  t  is  large  enough  to  make  the  argument 
of  the  exponential  near  unity.  For  larger  t,  the  factor  ~ C J- ^ q) 2q 
I  decrease  rapidly  again.  Thus  the  original  impulse  function  is  spread 
out  over  a  short  time  about  that  time  that  makes  the  argument  of  the 
exponential  equal  to  unity. 


The  form  of  Equation  4-116  can  be  simplified  by  calling  the  argu¬ 


ment  of  the  exponential  -u  ,  i.e., 


u(t.z)  s  (H5.) 


(l-q)/2q 


(4-117) 


Note  then  that 


3u  1  1  f.paz..  ^'1 

2  V  qaz  '  t 

Comparison  of  this  expression  with  Equation  4-116  shows  that 

i  -  /  2  3u  -u2 

im  V  irp  Dt  ° 


(4-118) 


(4-119) 


Now,  the  response  I  of  the  cable  to  u  unit  stop  function  is  related  to 

I .  by 
i  m 


(t,z) 


0 


-U  , 

e  du  . 


(4-120) 


u(t,z) 


Note  that  t  = 


corresponds  to  u  =  0,  so  that 


W"'*1  ■  /e‘u  du  *  vk  ■  (4-m> 

0 

Actually,  Ig^(°°,z)  should  equal  unity.  The  factor  1/V"2p  w  0.73  is  not 
far  from  unity,  but  our  estimate  is  not  entirely  accurate.  It  will  be 
shown  below  that  the  error  comes  at  late  times  after  most  of  the  pulse  has 
arrived  at  the  point  z.  For  the  present  we  proceed  with  the  formulae  as 
they  stand. 


The  integral  in  liquation  4-120  attains  half  of  its  final  value  at 
about  u  =  1/2,  and  it  can  be  said  that  the  signal  arrives  at  the  correspond¬ 
ing  time.  This  time  can  be  evaluated  from  liquation  4-117,  which  gives  the 
arrival  time 

t  =  p(4q)q/p(az)1/p  .  (4-122) 

Ll 

If  a  is  evaluated  from  the  reference  values  and  liquation  4-108, 

this  result  becomes 


v„t 
0  a 


'(4qk()z)q/p  , 


(4-123) 


where  v^  is  the  phase  speed  at  the  reference  values, 


V0  “  Vko 


(4-124) 


liquation  4-123  shows  that  the  arrival  time  increases  faster  than  z/v^  with 
increasing  z;  this  comes  from  attenuation  of  the  higher  frequencies  with 
increasing  z. 


The  rate  of  rise  of  the  signal  at  the  arrival  time  can  be  esti¬ 
mated  from  liquation  4-119: 


111 


ii  u,„ 


'  "T,1  'VffV  '»»-< 


rise  rate  «,  JX  (-  |l)e-!/4  ,  JI  „  c->/4 

\  7Tp  v  3t'  \TTp  t,  a 


0.62 


.VjL 


(4-125) 


The  rise  time  t  is  the  ratio  of  the  final  value  to  the;  rise  rate, 


t  pst  1 . 14  ^  t 
r  pa 


Thus  t  is  a  small  fraction  of  t  . 
r  a 


(4-126) 


Equations  4-128  and  4-125  have  a  very  simple  and  useful  interpreta¬ 
tion,  from  which  they  could  have  been  foreseen.  Firs.,  noco  thJt  nu  specific 
choice  of  and  has  been  made,  except  that  they  go  together.  Thus 

in  Equation  4-125,  k()  can  be  set  equal  to  any  value  of  ky  (the  real  part 
of  k)  if  vfi  is  set  equal  to  the  phase  speed  v^  going  with  k^.  The 
equation  can  then  he  written 


t  ■  (pP^Jqk  I)1'7'’  . 
a  v,  *  1  r 

0 


(4-127) 


Now, 


,p/q 


tP/q 


p‘  '  ‘  =  (l-q)l''M  ^  e”*1  (q  «  1)  . 

For  0.9  <  p  <  0.95,  4c”'1  is  quite  close  to  ir/2  (within  3  percent).  Next, 
note  that  according  to  Equation  4-107, 


qk  w  k.  . 

2  r  l 

Thus  liquation  4-127  is  equivalent  to 

a  v ,  i 


(4-128) 


(4-129) 


The  interpretation  is  now  clear:  find  that  u^,  from  Figure  4-8  for 
example,  for  which 


k^WjJz  =  1  ; 


(4-130) 
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These  results  hold  for  a  step  function  input  current,  figure  4-7 
shows  that  short  I2M  pulses  produce  step-like  currents,  but  with  finite  rise 
rates.  The  rise  rate  at  z  cannot  be  faster  than  that  at  the  input.  Lot'.g 
r.M  pulses  produce  ramp- like  currents.  The  response  of  the  cable  to  a  ramp 
is  the  time  integral  of  the  step-function  response.  The  arrival  time  for  a 
ramp  is  essentially  the  same  as  for  a  .step  function.  The  rise  time  of  the 
step  corresponds  co  the  time  for  the  ramp  to  acquire  its  final  slope. 


The  failure  of  I  to  reach  unity,  liquation  4-121,  comes  from 
the  fact  that  the  expansion  of  the  argument,  liquation  4-115,  is  not  accurate 
at  late  times.  For  t  >  0,  the  integral  in  Equation  4-109  can  be  evaluated 
by  folding  the  contour  of  integration  about  the  cut  as  in  Figure  4- 10b. 

The  result  is 
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(4-134 


I. 

xm 


(t,z) 


l  J'e-3V°s<q*>»«1’,,. in[si„(^)a«P]d3  . 
0 


For  large  t,  most  of  the  contribution  to  this  integral  comes  from  small 
s,  so  that  the  functions  of  s^3  can  be  replaced  by  the  first  terms  in 
their  power  series  expansions.  Thus  for  large  t, 

00 

lim(t,z)  ^  i  fe“St sin(qir)azspds 
C 


n  qaz(p!)/t^+^  (t  large)  .  (4-135' 

This  result  shows  that  the  impulse  response  falls  to  zero  somewhat  more 
slowly  at  late  times  than  was  indicated  by  liquation  4-116.  This  explains 
why  the  time  integral  of  liquation  4-116  did  not  quite  roach  unity.  Tho 
behuvior  of  I  t  is  sketched  in  Figure  4-11.  Tho  current  rises  rapidly 
at  t  near  t  ,  but  not  quite  to  unity.  The  final  rise  to  unity  takes 
a  fow  tu's. 


Figure  4-11.  Shape  of  the  step  function  response. 
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CHAPTER  5 

COUPLING  TO  LONG  BURIED  CABLES;  AN  EXAMPLE 


5.1  INTRODUCTION 

In  this  chapter  the  theory  of  Chapter  4  will  be  extended  to  buried 
wires  that  arc  long  enough  that  the  amplitude  of  the  BMP  changes  appreciably 
over  the  length  of  the  wire.  In  most  of  this  chapter,  attention  is  given 
to  a  particular  example,  namely  that  of  a  buried  power  line  which  terminates 
at  a  buried  facility.  In  order  to  maximize  the  coupling,  the  nuclear  burst 
is  assumed  to  occur  directly  on  the  power  line,  at  u  distance  of  1  kilometer 
from  tlio  facility.  The  geometry  is  sketched  in  figure  5-1. 


The  air  in  the  fireball  is  very  hot,  witli  temperatures  in  the  range 
1  to  10  eV,  and  is  therefore  thermally  ionized.  The  electrical  conductivity  in 
the  fireball  is  in  excess  of  lo‘  mho/meter.  The  fireball  is  a  very  good 
conductor  compared  with  the  soil  and  with  the  air  outside  the  fireball.  The 
radius  R  of  the  fireball  increases  with  time  as 

2/5  1/5 

R (meters)  1300  t  '  \  '  ,  (5-1  ) 

where  t  is  the  time  in  seconds  and  V  is  the  yield  in  megatons.  This 

1  /3 

formula  is  valid  for  t  -  0.1  V  '  second. 

buried  power  lines  typically  have  a  central  "hot"  wire  surrounded 
by  insulation,  then  wrapped  with  several  return  conductors  which  are 
approximately  at  ground  potential.  Often  a  partially  conducting  plastic 
sheath  protects  the  return  wires  from  corrosion  by  the  soil.  IV e  shall  assume 
that  the  return  wires  are  in  eloctiieal  contact  with  the  soil,  and  take  the 
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Fireball 


Air 

Ground 

Facility 

Figure  5-1.  Geometry  of  burst,  power  line  and  facility. 

resistance  of  the  return  wires  to  be 

()  ^  0.3  milliohm/meter  .  ( - 2 ) 

The  integral  of  the  radial  liMP  electric  field  in  the  ground  between 
the  fireball  and  the  facility  would  be  of  the  order  o'  1  to  10  mogavo Its,  if 
the  power  line  were  not  present.  After  a  period  of  iniuctivc  limitation  of 
current  in  the  wire,  this  voltage  appears  partly  between  the  wire  and  the 
lirobaJl.  at  one  end,  and  partly  between  the  wire  and  t he  facility  at  the 
other  end.  It  is  likely  that  these  voltages  arc  large  enough  to  drive  arcs 
between  the  wire  and  the  fireball  and  between  the  wire  and  t he  facility 
walls.  It  is  assumed  1  ere  that  the  latter  arc  reinforced  concrete,  and 
have  a  low  impedance  to  distant  ground.  The  fireball,  which  is  in  contact 
with  the  soil,  also  has  a  low  impedance  to  distant  ground.  The  current  in 
the  wire  is  limited,  after  the  inductive  phase,  by  these  two  impedances.  It 
is  likely  also  that  the  soil  will  break  down  in  the  vicinity  of  the  wire 
along  its  length,  reducing  to  some  extent  the  voltage  and  current  delivered 
to  the  facility. 

The  question  as  to  whether  ground  shock  destroys  the  power  line  is 
immaterial.  There  is  no  shock  wave  outside  the  fireball  (at  the  times  of 
interest  here),  so  the  power  line  must  be  intact  at  the  fireball  radius.  If  the 
power  line  is  opened  at  some  point  underneath  the  fireball,  the  arc  can  still 
strike  near  its  edge. 


116 


5.2 


THE  DRIVING  ELECT  RIC  FIELD 


Figure  5-2  shows  a  crude  representation  of  the  time  and  radius 
dependence  of  the  radial  electric  field  in  the  ground  within  a  few  meters  of 
the  surface.  The  times  T  labeling  the  curves  are  retarued  times, 

f  =  t  -  -  .  (5-.Vi 

Actual  calculated  fields  do  r.ot  fall  precisely  exponentially  with  distance, 
but  the  representation  shown  is  not  a  bad  one.  'lhe  fields  given  are  most 
appropriate  for  a  few  megaton  explosion  over  soil  of  conductivity  (at  low 
frequency' 

0^  ~  lO**5  mho  /meter.  (5-ii 

The  ends  of  the  curves  are  placed  at  the  fireball  radius  obtained  from 
liquation  5-1  with  Y  =  S  iMT. 


The  behavior  of  the  field  can  be  understood  from  the  theory  pre¬ 
sented  in  Chapter  3.  At  early  tin.es  the  air  conductivity  is  larger  than  the 
soil  conductivity,  ,:nd  the  radial  electric  field  is  Approximately  equal  to 
the  saturated  field  .J./o(air),  which  varies  little  with  distance  out  to  1 
kilometer.  At  late  times,  in  the  quasistatic  phase,  the  Compton  current 
passing  outwards  through  the  hemisphere  of  radius  r  in  the  air  returns  as 
conduction  current  through  the  hemisphere  of  radius  r  in  the  ground,  so 
c  ha  t 

-  r  /  X  1 

H  ,  f-  -  .1  ./'■>  -  e  /_  /r“  .  (5-5') 

t  .S  U 


This  formula  indicates  tli.it  F  should  fall  In  about  a  factor  of  10  between 

r 

r  -  300  and  000  meters  P  500  meters),  .in  agreement  with  the  curve  at 


'!'  -  10  “  second.  At  intermed  iale  times,  the  curves  can  be  understood  by 
assuming  tiiat  the  Compton  current  within  one  skin  depth  in  the  air  returns 
as  conduction  current  within  one  skin  depth  m  the  group  i. 
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The  fields  of  Figure  5-2  are  represented  analyticaly  by  the  formula 


Er(r,t)  =  &{ t  -  ^)e'6r  .  (5-6) 

Here  S  is  a  function  only  of  retarded  time.  The  parameter  3  is  not 
independent  of  time,  because  the  slope  cf  the  curves  in  Figure  5-2  varies. 
From  the  slopes  we  obtain  the  following  table: 

-6  -5  -4  -5  -2 

t  =  10  10  -  10  10  10 

3  =  0.007  0.0012  0.0034  0.0054  0.0079 

It  is  seen  that  3  changes  only  slowly  with  time. 


set 

-1 

m 


(5-7) 


If  6  is  regarded  as  independent  of  time,  then  the  Fouriei  expansion 
of  Er(rit)  contains  terms  of  the  form 

E2~ei(“t-kOi  (5-81 

where 

k  =  £  -  j3  .  (5-9) 

Note  that  the  horizontal  coordinate  (along  the  power  line)  is  now  called 

z,  in  order  not  to  confuse  it  with  the  cylindrical  coordinate  r  of  Chapter  4. 


5.3  THE  WIRE  IMPEDANCE 


The  impedance  of  the  power  1  i no  to  driving  fields  of  t lie  type  of 

Equation  5-8  is  given  by  Equation  4-28.  The  absence  of  an  insulating  layer 

in  the  present  case  can  be  accounted  for  by  putting  a^  =  The  impedance 

Z  is  then 
w 


Z  (u>) 
w 


2ir 


•  -2-  (£n  -  4)] 

n,  Y>c2a2  2 


(5-10) 


)  1 9 


I  {#* 


The  quantity  is  given  by  Equation  4-11,  which  with  Equation  5-9  becomes 

->  Jwru 


(t  -  WY  > 


(5-11) 


where  the  relative  soil  admittance  is 


,  .  'JJt. 

'o°  ■*  T 


(5-12) 


In  deriving  the  time-varying  inductance  model  in  Section  4. fa,  we 
? 

showed  that  the  second  (k*-)  term  on  the  right  in  Equation  5-11,  with  B  =  0, 
was  small  compared  with  the  first  term  and  could  be  neglected,  inis  approx! 
mat  ion  must  now  be  reexamined.  Writing  cut  Equation  5-11  gives 


K 


Uj  ,  ,  .  .  U)  ,  _  ,,,  , 

— -r  (c - 1)  *■  j  -  (3()o  -  2d)  -  ts 

c  ~ 


(5-15) 


The  contributions  of  the  k"  term  are  evident.  Eirst.  there  is  the  term 
-1  in  the  factor  i  -  1,  which  can  be  neglected  as  in  Section  4.0.  Second, 


there  is  the  term 


in  the  factor  Z^o  -  2f>.  Now 


‘0 


o  ~  0.577  while 


2f>  <  0.010;  thus  20  can  also  be  neglected,  finally,  there  is  the  term 


-  r< 


i’h  i  s  t  e  rm  is  negligible  i  f 
I’" 


0) 


c  0 


<u' ,  s  i  nee  ...  , '  c  -  i-i  . ,  i 
O  0 


■  Vw 


O(0 


skin  depth  . 


(5-14) 


thus  the  time-van.  ing  inductance  model  is  still  valid  if  the  distance  of  dif¬ 
fusion  of  fields  along  the  wire  is  small  compared  with  the  distance  1/(5 
in  which  the  amplitude  of  the  driving  field  vanes  appreciably.  In  the  time 
domain,  the  inductance  mode)  is  valid  for  times 


V 


!  0 


sec  in  present  example 


(5-15) 


120 


>  JvVO.l.lVjilMl-, 


,.j  i  o.i,  ...I  .I  t  .;  j, 


ofdlM  iVifriVi  I 


1  n  order  to  find  an  approximate  evaluation  of  Z  (,  iv.  is  necessary 

2  W 

to  evaluate  noth  the  factor  and  the  logarithm  in  Equation  5-10. 

For  the  log  term,  the  variation  of  which  is  not  sensitive,  we  take 


2  _  _5_ 

V | k2  j  a2  a2  * 


(S-lh) 


where  6  is  given  by  Equations  4-74  or  4-81  if  Equation  5-15  is  satisfied,  and 


1.12 

6" 


for 


t 


0_ 

2 


(5-17) 


For  the  factor  k 2/n  2 »  notc  that  the  neglect  ion,  in  Equation  5-13,  of  -1 
compared  with  c  and  of  -28  compared  with  Z^o  is  equivalent  to  writing 
Equation  5-11  as 


jwn. 
2  c 


P2  , 


or 


k:  r2  „2 

_ 2  r4  :  w  .  b_  w  .  w  _ 

n2  J  c  n7  J  c  "  zoo 


(5-18) 


The  second  form  here  recognizes  the  dominance  of  Z  0  over  ao/c  at  late 

2  1 
times  when  the  3  term  is  significant. 


With  these  results  and  with  the  inclusion  of  the  resistance  E  ,, 

w() 

of  the  wire,  the  impedance  of  the  wire  can  be  written 


Z  (■■))  =  jojL  +  R  n  -  R  ,  (5-19; 

w  wO  3 


wlie  re 


►  -  n  5. 

L  •-  -  £ti(-(— )  Menry/mctcr  ,  (5-20) 


and 


JL 
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R,  =  ■=—  JLn(; — )  ohm/meter  . 

P  Z  HU  £i  ^ 


(S-21) 


As  in  Section  4.6,  the  time  domain  equation  for  the  wire  current  1  is 


[L(t)Ift)]  -  E(t)  -  [Rw0-R3(t)]I 


(5-22) 


Note  that  the  resistance  R^  represents  a  negative  resistance,  which 


accounts  for  the  increase  in  current  at  a  given  z  due  to  larger  driving 
field  at  smaller  z.  The  e-folding  time  of  the  current  due  to  this  negative 


resistance  is 


L  V 


R  '  2  * 

0  3 


(5-23) 


-4  -4 

For  the  example  considered  in  this  chapter,  t  10  second  at  t  =  10  second. 


Thus  diffusion  along  the  wire  accounts  for  most  of  the  increase  in  current 


after  t  =  10  seconu. 


The  solution  of  liquation  5-22  can  be  written  in  terms  of  integrals. 


(5-24) 


Then  liquation  5-22  can  be  written  as 


[e"g(t)L(t)I(t)]  -  e“8U(t)  , 


and  the  solution  of  this  equation  .is 


I(t)  = 


ii(t 1  )dt 1 


(5-25) 


.1  .. .  - 


v-  '  v 


- magma 


1 1 r  j a*. ■* i.-,* ".  ''i )  e  '!  ■* . r 


5.4 


DIFFUSION  ALONG  THE  W IKE 


The  physical  origin  of  the  term  R^I  is  made  clearer  by  writing 


V 


LI  = 


-  _3 


V 


LI 


(5-26) 


Now  it  was  assumed  in  Equation  5-6  that  the  dependence  of  the  driving  field 

on  z  and  t  is  <g’(t  -  z/c)e”^Z,  and  it.  was  shown  above  that  the  dependence 

of  on  z  produces  negligible  effects  when,  as  in  the  case  of  interest 

here,  the  wire  is  in  electrical  contact  with  the  soil.  Thus  the  essential 

"•  R  2 

dependence  of  E  on  z  is  e  ,  and  I  will  have  the  same  dependence, 
apart  from  end  effects  which  are  considered  below.  Therefore, 


(5-27) 


so  that  Equation  5-22  is  equivalent  to 


,  2 

[LI]  =  E  -  R  ,nT  + 

Of  w 0  |J ^0  g „ 


[LI] 


(5-28) 


This  diffusion  equation  shows  explicitly  that  the  magnetic  flux  LJ  per 
unit  length  of  wire  diffuses  along  the  wire,  and  is  valid  for  any  dependence 
of  E  on  z,  not  just  exponentials. 


ation  5-28  can  be  derived  directly  from  Maxwell's  equations  for 
the  present  case.  These  are,  in  the  cylindrical  coordinates  of  Figure  A-2, 


9Bn  31i  DE 
_0  z  r 

at  ' 


'  0  r 


Sr  3z  » 

3B0 
3z  ’ 


..  1  3  „ 

pnai,  =  -  v  ‘  i'Bn  . 

0  z  r  3r  0 


(5-29) 

(5-50) 

(5-51) 
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In  the  last  two  equations  the  displacement  current  has  been  neglected,  a 
good  approximation  in  most  soils  except  at  the  highest  frequencies  or 
very  early  times.  Using  liquation  5-30  in  liquation  5-29  gives 


8Be  3E,  .  i  r»o 

-3F  ■  IF  V  -,-p- 


(5-32) 


Let  the  fields  indicated  here  be  those  due  to  the  current  in  the  wire,  not 
including  the  incident  fields  applied  to  the  wire.  Then  Bq  and  B  vanish 
at  sufficiently  large  r.  Integrating  liquation  5-32  over  r  from  the  radius 
a^  of  the  wire  to  large  r  yields 


,  .  1  D“(|) 

'z  ^U2  +  jj  o  2  ’ 

U  <JZ 


(5-33) 


where  i|j  is  the  magnetic  flux  per  unit  length  of  wire.  If  the  wire  were 
perfectly  conducting,  then  11,  (a. ,)  =  -  ii  (the  applied  field);  for  resistive 
wire, 


I!  fa,)  *  -  li  ♦  Rnl 


w()  * 


(5-3  .n 


The  flux  (|)  is  estimated  by  assuming  that  is  the  field  of  the  current 

I  out  to  the  skin  depth  <5 , 


*  "  27f  1  /  T  =  27,  ■  1 


(5  35) 


When  liquations  5-34  and  5-35  are  used  in  liquation  5-33,  the  result  is  liquation 
5-28.  This  derivation  docs  not  give  the  formula  for  d,  which  comes  only 
from  the  complete  solution  of  the  equations  developed  in  previous  sections. 

The  complete  solution  also  makes  clearer  the  effects  of  approximations  made. 

5.5  TERMINATION  CONDITIONS 

IT.  order  to  determine  a  solution  of  liquation  5-28  for  a  wire  of 
finite  length,  an  additional  equation  is  needed  at  each  end  of  the  wire.  This 


*  i  l  '^|* 1  |'H  *  ’( J  i  12.  5 : ‘J  C,  t  .O-  r  /  .  L-.Q  lA  tO 


equation  can  be  obtained  by  integrating  liquation  5-30  over  r  from  an  to 
6.  The  I'esult  is 


where 


Vv  -  -  5 


■  /v>  ■ 


c)(J)  _  3  r 


[LI]  , 


(5-36) 


(5-37) 


is  the  voltage  between  the  end  of  the  wire  and  distant  ground  associated  with 
the  current  I  in  the  wire;  i.c.,  V  does  not  include  the  incident  electric 
field.  Usually,  V  is  related  to  1  by  a  simple  impedance,  which  at  all 
but  the  earliest  times  is  well  approximated  by  a  resistance.  If  the  facility 
in  figure  5-1  is  approximated  as  a  conducting  sphere  of  radius  a^.,  its 
resistance  to  distant  ground  is 


w  2  f  ..--.V...  re  ..A.  .. 
!  J  4m  r2o 


(u  »  a.,J 


(5-38) 


Tnc  factor  2  hero  comes  from  the  fact  that  the  facility  is  located  near  the 
surface  of  the  semi- in finite  soil  medium.  The  resistance  of  the  fireball  to 


■  stant  ground  is  that  of  a  disc  of  ratlins  a...  , 


R...  «  -  .  (d  »  a,..  ) 

i  I)  'la...  o  '  lb 

in 


( 5-3‘J) 


11'  the  load  resistance  is  designated  by  R,  (=  R  or  ) ,  then  the  relation 
between  V  and  1  is 


V  =  irl  ’ 


(5-40) 


where  I  is  the  current  flowing  ou_t_  of  the  wire  into  the  soil.  Combining 
this  equation  with  liquation  5-36  yields  the  termination  condition 


l  [Ml  - 


Kn(~-) 

tl  ^ 


(5-41) 
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5.6 


METHOD  OF  SOLUTION  OF  THE  EQUATIONS 


The  Flux  function  has  been  defined  (by  Equation  5-35)  as 

<|>(2,t)  L(t)  I  (z , t)  .  (5-42) 

Here  L  is  given  by  Equation  5-20,  with  6(t)  given  by  Equation  4-81  or 
Equation  5-17 ;  the  t  in  <5(t)  is  really  the  retarded  time,  or  time 
after  the  arrival  of  the  driving  field  at  the  position  z.  The  difference 
between  using  real  time  or  retarded  time  in  the  rest  of  the  equations  was 
shown  in  Section  5.3  to  be  negligible.  We  therefore  think  of  t  as  retarded 
time  in  the  remainder  of  this  chapter. 


These  equations  can  be  solved  quite  readily  by  finite-difference 
methods.  However,  approximate  solutions  can  also  be  found  analytically, 
and  these  are  useful  for  providing  understanding  and  checks  on  the  finite— 
d.i  ffcmiec  ivsul  ts . 


The  first  step  in  the  analytical  solution  is  to  eliminate  the 


term  R  n(f>/i  by  defining 


m  r  Rwo 
’1^  '  J  L(t* 


f^lf  , 


(5-46) 


and  letting 


(f>  =  e 


"8i  (t) 


(5-47) 


Then  Equation  5-43  becomes 


3<l,1  8,(t)  ,  A, 

■*r-°  li  +  v77* 


(5-48) 


Equations  5-44  are  left  unchanged  except  that  <|>  is  replaced  by  <J)^ . 

A  goner i.l  method  of  solving  Equation  5-48,  subject  to  the  end 
conditions  is:  first,  find  a  particular  solution  of  Equation  5-48  ignoring 
the  cud  conditions;  second,  find  solutions  to  the  homogeneous  equation  obtained 
by  setting  E  *  0,  again  ignoring  the  end  conditions;  third,  choose  a  linear 
combination  of  the  particular  and  homogeneous  solutions  which  satisfies  the 
end  conditions. 


=  r£"(t)  e"  , 


(5 -4D) 


then  a  pai  joe lav  solution  of  Equation  5-48  can  be  found  by  assuming 


<J>I  U.t)  =  <l>2  (t)  e 


(5-50) 


•Substitution  of  this  form  leads  to 


d'k  8,  ft)  p2 

~r  =  0  1  t)  +  <|._  . 

dt  p()0  2 


(5-51) 


Strictly  speaking,  this  procedure  is  valid  only  if  fi  is  a  constant.  How¬ 
ever,  liquation  5-51  is  approximately  correct  if  H  changes  slowly  with  time, 
which  we  assume.  Then  if  g  (t)  is  defined  by 


M1'  “./  a*’  ■ 


C  5  -  5  2 ) 


the  solution  of  liuuntion  5 - 5 1  is 


BiCt.)  r 

e  ”  J  explicit 


')  "  g.M'J  |/'"(t’)dl  '  .  f 5 - 5 3 ) 


This  completes  the  particular  solution  of  liquation  5  I M . 


5.7 


HOMOGENEOUS  SOLUTIONS 


In  fliulin)'  solutions  of  the  homogeneous 
convenient  to  redefine  the  space  variable  as 


form  of  liqua  t  i  on  !>  •  IS  ,  it  i  s 


y  :=  \fj^  . 

Then  the  equation  becomes  simply 


ih  I) 


h'l 


Uv‘ 


where  we  have  des  i  gnu  l  evl  the  homoj'.oneous  solution  hy  T  in  orde 
distinguish  it  I'roii  the  part  ieti  I  fir  solution  <f*j .  A  very  simple 
of  liquation  5 - 5 5  is 

at  •  <J7  v 

I  i  •  /i  /»  ' 


c  V 


where  s  is  an  arbitrary  positive  constant .  This  would  bo  a  us 
homogeneous  solution  to  add  to  the  particular  solution  -fj  if 
increased  exponentially  with  time,  liquation  f>-T>o  indicates  that 
<|)j  will  increase  approximately  exponentially  with  time  a i  late 


( 5  -  55) 

r  t  o 

solut ion 

( 5-5d I 
efu  I 

,i. 

1  1 

•I',  and 
*  i  mes . 
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At  each  end  of  the  wire,  we  should  choose  that  solution  that 
decays  with  Lnevensinjj  distance  into  the  wire.  Thus  at  the  ends  of  the 
wire  for  the  solution  yoiny  as  eSt,  we  have  the  relation  between  p  and 
D p/  D  z . 


Dip 

Dz 


V^i'oo  P 


at  rip.ht-hand  end  , 
at  left-hand  end  . 


(5-P7) 


it  will  turn  out  that  only  the  relation  between  p  and  Dp/D;’,  at  the  end 
of  the  wire  is  needed  to  calculate  the  total  current  at  the  end  of  the  wire, 
nt  those  times  before  diffusion  can  occur  over  the  entire  length  of  the  wire. 


At  early  times,  liquation  f>-f>.'5  indicates  that  p ,  will  vary  more 
like  a  power  of  the  time  than  exponent  dally .  Hence  It  would  he  useful  to 
find  homoj'.onuous  solutions  p  such  that 


PCt.ynO)  t" 


T»-r,8) 


where  11  is  a  positive  constant.  Such  solutions  are  conveniently  found 
l»v  Laplace  transform  of  the  time  variable  in  liquation  li-55.  It  is  easily 
shown  that  the  peneral  solution,  appropriate  to  the  rij'.ht-hand  end 
of  the  wire,  in  the  Lap  lace  domain  is 


P  (  s  ,  y )  ~  I’ (s')  e  '  , 

where  f(s)  is  an  arbitrary  function  of  the  Laplace  variable  s 
function  is  fixed  by  liquation  S-58.  At  y  ~  0,  f(s)  must  be  the 
transform  of  tn, 


fft-51)) 

Till  s 
Lnplace 


fj 


(T> -(>()) 


The  Laplace  transform  of  Dp/by  can  he  found  by  takinj;  the  derivative  of 
liquation  !>-!>;).  At  y  ■-  t), 
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Comparing  Equations  5-60  and  5-61  makes  it  clear  that  in  the  time  domain 


dM±>2)  \ 

*  I y=o 


_ n!.. _ t  n  - 1  /  2 

(n  -  |) ! 


(5-62) 


Tlius  the  relation  between  ip  and  d\p/dz  at  the  right-hand  end,  for  power 
law  time  dependence,  is 

'I  *  T~  VV77  f  '  (5-"5) 

(n  -  -p  • 

At  the  left-hand  end,  a  minus  sign  should  be  inserted  in  this  equation, 
figure  5-3  contains  a  graph  of  the  ratio  of  factorials  that  occurs  in 
liquation  5-63. 


The  time  dependence  of  <|>2 ,  Equation  5-53,  will  rarely  lie  purely 
exponential  or  purely  power  law.  However,  the  relation  between  iji  and 
dip/dz  is  not  very  sensitive  to  the  precise  form  of  the  time  dependence, 
for  example,  the  exponential  funccion  e'st  is  tangent  to  varying  power 
laws  at  varying  times.  Since 


s  (t  +cSt) 


w  l+sfit)  , 


(5-64) 


and 


i 

(t+iSt) "  ~  tn(l  +  ~  fit)  , 


(5-65) 


it  is  clear  that  at  time  t  the  exponential  is  tangent  to  a  power  law  with 


n  -  st  or  s  -  n/t  .  (5-66) 

If  in  the  end  relation  5-57  for  exponentials  we  make  the  replacement  s  =  n, 
we  obtain  the  approximate  power  law  result 
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The  crosses  plotted  in  figure  5-3  ;ire  which  is  to  be  compared  with  the 
ratio  of  factorials  of  liquation  5-63.  It  is  seen  that  treating  the 
exponential  case  by  power  law  fits  at  various  times  leads  to  only  a  small 
error . 

We  shall  write  the  relation  between  ig  and  0-.,' ,/ :» is 


J  @  -  *  A(t)  . 


(5-68) 


where  the  plus  sign  is  for  the  right-hand  end,  the  minus  sign  is  for  the 
left.  The  formulae  for  A  are 


st 


A  =  y  s;i  o  for  h  ~  o'  , 


(5-6D) 


. — -f —  v  li,.0'/t  for  <l>  ~  tn  . 

(n  -  ! 


Note  that  1/A  is  approx i mat  el v  the  distance  di Mused  along  the  wire. 


5.8  SOLUTION  OF  THE  TERMINATION  GONDII  ION 

FOR  EARLY  TIMES 


At  times  sufficiently  early  that  the  diffusion  distance  1/A  is 
small  compared  with  the  length  of  the  wire  the  contribution  of  it  ,  the 
homogeneous  solution  associated  with  the  left-hand  end,  to  the  flux  at  the 
right -hand  end  is  small  and  vice  versa  for  ij;  .  The  solution  of  the  tormina 
tion  condition  liquation  5-45  is  then  relatively  simple.  Having  evaluated 
liquation  5-55  for  the  part  r.lar  solution  ^(t)  at  the  right -band  end,  we 
find  the  power  aft)  appropriate  for  each  value  of  t  .  Let  ft)  lie 
the  value  of  the  light -hand  homogeneous  solution,  including  an  arbitrary 
constant  multiplier  which  is  absorbed  into  ifi  .  We  then  unsider  the  total 
so  1 ut i on 

13? 


<!>T  =  <f>2  +  liJr  ■ 


(5-70) 


exceeds  the  decay  length  1/(5  of  the  driving  electric  field.  No  change 
of  sign  occurs  because  A  approaches  3  at  these  times.  Equations  5-53 
and  5-52  show  that 


,  I  \ 

*2  exI’UT0V  t; 


(5-79) 


when  fTt/y^o  >>  1,  or  when  the  diffusion  distance  \  t/y^u  >>  1/fi.  Thus 
at  those  times  is  approximately  exponential  in  t  with  s  6“/y  o, 
and  from  Equation  5-69, 


A  ^  3 


(5-80) 


The  factor  A-3  in  Equation  5-78  becomes  small  as  <t>^  becomes  large 
(exponentially),  and  the  result  for  d)  is  net  clear. 

in  order  to  resolve  the  uncertainty,  it  is  necessary  to  examine 
Equation  5-78  more  carefully  in  the  Laplace  domain,  since  d>0  is  not  exactly 
exponential  in  time.  In  the  Laplace  domain,  the  equation 


\j  Uras  -  B 

4-  (s)  =  —  - <M*0  , 

Vo05  -  r-i 


(5-81) 


is  correct  provided  i 8  properly  evaluated.  This  can  be  achieved 

from  the  Laplace  transform  of  Equation  5-51,  which  is 


8<fS  =  «:(s) 

"V  “ 


(5-82) 


where  E(s)  is  the  Laplace  transform  of  the  function 


Fit)  F  e 


glflV(t)  . 


(5-83) 


The  solution  of  Equation  5-82  is 


>VJ 

4>,(s)  = - y  l  (s)  , 

y^os  -  3 


(5-84) 


and  Equation  5-Si  becomes 


The  inversion  of  Equation  5-85  to  the  time  domain  can  be  done 
approximately  for  those  times  at  which  the  diffusion  distance  is  larger 
than  1/3,  where  Equation  5-78  is  insufficient.  At  these  times 
can  be  neglected  compared  with  8.  In  practical  cases  is  comparable 

to  8  or  ’arger,  so  that  s  cun  also  be  neglected  compared  with 

Then 

V*)  “  ^  "is)  , 

which  leads  immediately  to 

iJ0o  gut) 

^At)e  1  .  (5-86) 

The  current  at  the  left-hand  end  if 

-8i(t) 

i(t)  =  c  4y'h 

»<ff(t)/3RL,  (5-87) 

where  R  is  the  load  resistance  at  the  left-hand  end.  According  to  Equation 
5-49,  <?(t)/3  is  the  integral  of  the  electric  field  along  the  wire,  or  total 
voltage.  Most  of  the  voltage  drop  occurs  near  the  left-hand  end.  The  current 
is  approximately  equal  co  the  voltage  drop  divided  by  the  load  resistance 
in  the  time  frame  assumed  in  this  paragraph. 


At  the  right-hand  end,  Equation  5-76  contains  no  such  cancellations. 
This  equation  is  correct  until  diffusion  from  the  left-hand  end  can  reach  the 
right-hand  end. 
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5.9 


QUASI STATIC  SOLUTION  AT  LATE  TIMES 


At  Late  times  the  end  effects  generated  at  each  end  can  diffuse 
to  the  other  cna ,  and  this  interaction  of  the  end  effects  must  be  taken 
into  account.  Ts  do  this  we  return  to  the  general  equations  for  <J>, 
liquations  5-43  and  5-44,  and  go  to  the  limit  in  which  diffusion  is  rapid 
compared  with  the  time  t.  Since  li  varies  more  and  more  slowly  as  time 
increases,  the  term  3c|>/9t  in  liquation  5-43  becomes  small  compared  with  other 
terms,  and  can  be  neglected  in  first  order.  This  procedure  yields  the 
static  solution  appropriate  to  the  li  and  L  at  each  time.  We  therefore 
need  to  solve  the  equation 


D  “  cf> 


-> 

7“<l> 


-  ,;oal! 


(S-88) 


whore  the  quantity  y  is  defined  by 

T  ■  V'V'V- 


/2'frok 


/  4SIIUK  ,*  \ 

- m 


,1/2 


\ 


6\l 


(5-89) 


We  again  use  the  form  of  liquation  5-19  for  I!,  and  choose  the  origin  of  the 
z  coordinate  to  be  at  the  left-hand  end  of  the  wire,  bet  d  be  the  distance 
from  the  fireball  edge  to  the  facility,  or  length  of  wire  exposed  to  the 
field  H. 

The  general  solution  of  liquation  5-88  is 


>V’‘f 

"n  ’  o 

r  - 


~!')Z  -  y.,  ,, 

■e  +  (.  j  e  -  Ue 


,  +yz 

n 


(5-90) 


where  and  C-,  are  arbitrary  constants  which  can  be  chosen  to  satisfy 

the  end  conditions,  liquation  5-44.  These  end  conditions  are 
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6  -  YCj  +  yC2  =  r,t  |-1  +  Cl  +  C2|  , 


Wl! 


-Ixl 


-Yd,. 


■Jd, 


-  YC  ,~C1  +  yc  (’2 


r,  f-c  W  +  ( 
r 


I 

^  '  VYdl-j 


(5-91) 


These  equations  can  bo  solved  for  Cj  and  C0 ,  with  the  results 

C, 


'1 


[(Ca+3)(Cr+Y)cYd  -  (C£-Y)(Cr-3)e“(idj/b  , 


C2  “  t(Ca+Y)((,r-[i)e“W  -  (r,a+3)(?Jr-Y)o"Yd]/l)  , 
u  =  (ya+y)(^,+y)cy‘i  -  (c.£-y)  (ci,-y)o”Yc1  • 


1 

( 


(5-92) 


When  and  have  been  evaluated,  the  flux  <|i  can  be  calculated  from 

liquation  b-90,  as  a  function  of  z.  The  fluxes  <]>„  and  <|> ,  at  the  left¬ 


'll,  -  T 

and  rij’ht-hand  ends  can  be  evaluated  directly,  with  the  results 


li„o<f; 


(|,  =  — i— — --a-  [(H-Y)(^,+Yi'.Yd  -  tii+Y)(Cr-Y)o"Yl1  +  2y(h -H)e“W]/l)  ,  (9-910 

V'  H‘  -  y2  1  1  1 


..OfS’e 


-W 


<l>  “  A - yl _  ( 3+Y )  ( Cj/+Y )  e^1  +  (3-Y)  (Cj-Y)o"Yd  +  2Y(r.y+H)oW]/l»  .  if. -91) 

■■  y“ 


Note  that  the  brackets  in  these  equations  transform  .into  each  other  when 
f,v  and  ^  are  interchanged  and  I'  is  replaced  by  as  is  required  by 

symmut  ry . 


The  expressions  derived  above  are  fairly  comp  1  i cited  and  the 
results  are  not  easy  to  visualise.  However,  considerable  simplification 
occurs  for  a  case  that  is  of  practical  importance.  This  ease  is 

(xl  »  i  ,  Yd  1  ,  Y  «  r,g  ,  r,v  .  (h-95) 

When  these  conditions  hold,  terms  contain  iny,  a  factor  e  can  bo  dropped 
and  the  approximation 
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**  1  ±  yz  , 


(5-96) 


can  be  made.  Then  c|)(z)  can  be  shown  to  be 


I'nV  r  a7  (1  -3X-)  (d+X  -z) 

qHz)  **-v-  -c +  — -* — 1 — 

S2  l  '  ‘  H  +  \  J 


(5-117) 


where  wo  have  defined  the  "extrapolation  lengths" 
x&  =  l/>r'Sl  »  xr  =  J/,"r  ‘ 


'5-98) 


The  boundury  conditions  (5-44)  are  equivalent  to  the  statement  that:  if  t}> 
is  extrapolated  with  constant  slope  to  a  distance  l/{,  ,  or  1  /r  beyond 
the  ends  of  the  wire,  the  extrapolated  <()  must  vanish.  Liquut i on  5-97 
satisfies  this  extrapolation  condition  approximately.  At  z  =  d  +  x  ,,  <|> 

_n  |  1* 

vanishes  in  the  approximation  that  e"  L  is  negligible.  At  z  =  -x  ,  the 

-  f;>  z 

extrapolated  value  of  o  is  1  +  Hx^,  so  that  the  bracket  vanishes 
exactly,1  A  sketch  of  the  geometrical  relation  of  the  exponential  and  linear 
terms  in  the  bracket  is  shown  in  figure  5-4,  along  with  the  shape  of  Hz). 
The  value  of  the  current  at  the  right-hand  end  is 


rj)  „  (Kill  ,, _ giracJf . B  .  .<?.  .  1  *  Hx 


,.2,t  ,  i)  d  V“x’  +'  x 
H  f.n  (.,  -  •)  2.  r 

'a 


(ilij  '  [T(T+x7+xry  ’ 


(5-99) 


where  Uj  is  the  termination  resistance  of  the  right-hand  end.  Note  again 
that  rP/p>  is  the  total,  voltage  applied  along  the  wire. 

5.10  APPLICATION  Or  FORMULAE  TO  EXAMPLE 

The  approximate  theory  developed  above  is  applied  in  this  section 
to  the  example  defined  in  this  chapter.  The  parameters  of  this  example  are, 
in  summary: 


/: .  xj  Aa^Lik'fjau 


Figure  5-4.  Sketch  showing  relation  of  exponential  and 
linear  terms  in  Equation  5-97,  extrapolation 
lengths  and  flux  (|> ( 2 ) . 


wire  resistance  ~  It  ^  -  0.3  ohms/km; 

-  7\ 

ground  conductivity  -  0  =  o()  «  10  '  mho/rn  ; 
ground  permittivity  «  cu)  «  10  ; 
radius  of  wire  =  a.,  -  0.01  m  ; 

resistance,  facility  to  distant  g round  «  it j,  ■■■  10  ohms 
distance  from  burst,  point,  to  facility  =  1000  m; 


(!>•■  1  ut) ) 


driving  electric  field  given  by  Figure  5-2 


From  these  input  parameters  the  derived  parameters  calculated  at  the  indicated 

times  are  listed  in  the  Table  5-101.  In  this  table,  (3  is  calculated  from 

the  slope  of  the  curves  in  Figure  5-2.  The  skin  depth  6^  comes  from 

liquation  4-81,  which  takes  into  account  approximately  the  frequency  variation 

of  0  and  t:  for  typical  soils,  liquation  5-17  determines  5.,.  The  smaller 

of  wj  find  is  used  in  liquation  5-20  to  determine  the  inductance  L. 

The  resistance  between  the  fireball  and  distant  ground  comes  from 

liquation  5-30,  with  the  fireball  radius  a^  read  off  Figure  5-2.  The 

logarithmic  derivatives  «'  anu  t,  come  from  liquation  5-45,  and  the 

y,  i’ 

extrapolation  lengths  am!  x  come  from  liquation  5-08. 
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1  520 
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124 

145 

164 

158 

152 

The'  next  r.  t  i.-p  is  to  calculate  the  exponent:  is]  a  r  oumon  t jtj(t)  fiiul 

C,  ,  ft),  defined  by  liquat  i  ons  5-46  and  5  •  5  2  It  is  clour  from  I  U|tui  t  ions  5-47, 

5-48  and  H  -  f»  .'•>  tli, st  these  arj’uments  need  to  lie  calculated  accurate  ly  only  when 

they  are  not  small  compared  with  unity.  If  we  take  I,  =  <1.0  x  K)  p||/m,  a 

v'alue  appropriate  to  tile  period  10  to  10  “  second,  liquat  ion  f* - d < >  qivcs 


Kj  ft)  -o  ir,0  t  , 


O.lf.  at  t  =  l(f‘\  1.5  at  t  ■•=  1  ( > “ see  . 


(5- 1 02) 


_  ') 

liquation  5-102  is  tiierel'ore  an  adequate  approximation  up  to  10  “  see,  From 
liquation  5 -11 2  it  can  be  seen  that  i;.,(t)  reaches  the  value  unity 
approx  i  inn  t  e  ly  when 


I  or  t  p  ,o/|-'.“  -  1.25  x  .it)  '  /(•'■“  . 


(5-104) 
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With  B  read  from  the  Table  5-101,  trial  and  error  leads  to  t  If' 

second  as  the  time  when  $>..  nears  unity.  A  simple  formula  representing 

2  z  -5-5 

B  within  20  percent  an  the  time  interval  from  4  x  10  to  10  '  second  is 

B2  »-  10 ayft  , 

from  which  the  result  follows, 

t  V2 

"  0.53  .  (5-104) 

10 

in  order  to  calculate  the  early-time  current  at.  the  right' hand  end 
of  the  wire.  <f»0( t)  must  next  he  evaluated  from  liquation  5-53.  The  >'■  in 
that  equation  is  to  he  read  from  figure  5-2  at  1  kilometer  from  the  luirst 
point.  The  points  read  off  are  graphed  .in  figure  5-5  and  a  smooth  eurvo  is 
drawn  through  them.  This  curve  is  then  multipl  ied  by  cxpfg . -g,,)  ,  and  the 
result  is  graphed,  and  integrated  numerically,  yielding  the  integral 
curve  shown,  filial  Iv  this  curve  is  multiplied  by  oxp(g.,)  to  yield  <!>)(t). 
All  of  these  operations  can  be  done  graphically  and  with  the  aid  of  a  pocket 
calculator  in  a  half  hour  or  so. 

The  next,  step  is  to  apply  the  end  condition,  liquation  5-7(>,  to 
obtain  '|'-j. •  To  determine  the  logarithmic  derivative  A  we  use  the  power 
law  approx i ma i  ion  in  liqtiai  ion  5 -(d).  from  the  slope  of  ■I'-,  in  figure  i:-5, 

the  values  of  the  power  u  are  determined;  A  is  then  calculated  from 
figure  5-3  and  liquation  15 -(it) .  Next,  the  factor  (A+Bl/ (Ait,  i  is  evaluated 
and  <|>,|,  is  determined  from  liquation  5-7(>.  finally,  the  current.  I  is 
calculated  from  liquation  5-77.  The  numbers  obtained  in  those  operations 
a  i  e : 
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t ,  see 
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o 

X 

to 
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The  early-time  current  computed  here  is  graphed  in  Figure  5-6.  It 
begins  to  rise  exponentially,  due  to  diffusion  along  the  wire,  at  about  10" 
second,  as  expected. 


The  late-time,  or  quasistatic  approximation  for  the  current  is 
given  by  liquation  5-99.  In  that  equation,  is  the  electric  field  at  the 
fireball  end  of  the  wire,  and  d  is  the  distance  from  the  fireball  edge  to 
the  facility,  which  is  read  from  Figure  5-2.  All  of  the  other  parameters 
In  liquation  5-99  have  been  calculated  above.  The  numbers  are: 


t ,  sec 
,  V/m 
d ,  m 
1 ,  Amps 


A 

1.7  x  If) 
960 

I.D7  x  10 5 


8.8  x  K)-1 

you 

7.9  X  jo'1 


10-2 

6.5  x  1()“ 
745 

5.0  x  K) 


3 


(5-1 00 1 


This  qnasistat ic  current  is  also  graphed  in  Figure  5-0,  where  it 


is  seen  that  the  quasi stat  ic  current  is  less  than  the  early-time  (inductively 
limited)  current  after  t  =  1.3  ><  !()  *  seconds,  at  about  the  same  time 


t  =  1  .  .8  x  |  o 

that  the  oarlv-time  current  becomes  exponent ia 1 ,  This  result  may  seem 


surprising.  The  time  to  go  into  the  exponential  phase  is  the  time  to  dif¬ 
fuse  a  distance  1/f  -  300  meters  (at  10  '  second),  and  this  time  is 
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1*11*1'.  >2.  ii  ^i'K<  1 


I  (Amp 


(5-107) 


t  ft*  PqO/3^  ~  10  ^  second  . 

The  time  to  diffuse  a  distance  d  is 

t  *»  pgad2  ~  10"3  second  ,  (5-108) 

_  3 

so  that  one  might  have  expected  the  quasistatic  phase  to  start  only  at  10 
second.  However,  because  the  driving  field  is  much  larger  at  the  left-hand 
end,  only  a  small  part  of  the  left-hand  end  effect  needs  to  diffuse  to  the 
right-hand  end  to  make  a  noticeable  effect  there. 

The  actual  current  is  estimated  by  joining  the  early-time  and 
quasistatic  currents  smoothly,  as  indicated  by  the  dashed  curve  in  figure 
5-0.  Note,  however,  that  in  the  decade  between  10  ^  and  10  3  second, 
where  the  peak  current  occurs,  the  quasistatic  approximation  has  nut  been 
shown  to  be  reliable.  A  better  treatment  of  the  diffusion  of  the  left-hand 
end  effect  is  needed.  Such  a  treatment  could  be  devised. 

Note  that  the  energy  delivered  into  the  assumed  10-ohm  load  at  the 

8 

facility  is  of  the  order  of  10  Joules.  This  explains  the  extensive 
electrical  damage  that  occurred  in  bunkers  that  had  long  wires  going  into 
them  in  the  early  days  of  nuclear  testing. 
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CHAPTER  6 

COUPLING  TO  OVERHEAD  LINES 


6.1  INTRODUCTION 


Ground-based  systems  often  employ  overhead  lines  for  power  instead 
of  buried  lines  as  they  are  substantially  cheaper  to  build  than  buried 
lines.  These  lines  are  not  intended  to  operate  after  an  attack,  as  they  arc 
relatively  easily  damaged  by  air  blast.  However,  their  presence  can  result  in 
large  EMP  signals  which  are  generated  along  their  length  and  transmitted 
to  equipment  which  does  have  a  post-attack  survivability  requirement. 

Overhead  power  lines  often  consist  of  a  set  of  several  wires  carrying  power 
and  a  neutral  wire  elevated  above  the  rest  which  is  periodically  grounded 
for  protection  against  lightning — for  simplicity  in  this  section  we  will 
discuss  coupling  to  a  single  line  located  at  a  height  of  10  meters  above 
the  ground. 


Two  features  complicate  the  theory  of  coupling  to  overhead  lines 
compared  to  buried  lines. 

1.  The  conductivity  of  the  medium  surrounding  the  wire  varies 
as  a  function  of  both  time  and  distance  from  the  burst. 

2.  'I  he  boundaries  of  the  wire  and  the  ground- air  interface  do 
not  fit  as  coordinate  surfaces  in  a  system  where  the  Helmholtz 
equat ion 

fV2  +  kJ)<|>  0  ,  fo-1) 

is  separable.  The  Laplace  equation  in  the  two  transverse 
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dimensions  (where  the  z  axis  is  para  Mel  to  the  wire) 


+  --)$  -  0  ,  (6-2) 
W  3y  ■ 


is  separable  in  bipolar  coordinates  and  we  will  exploit 
this  feature  in  deriving  transmission  line  equations  suitable 
for  late-time  calculations. 


»> 


The  discussion  of  current  on  the  overhead  line  system  breaks 
naturally  into  two  physical  regimes  —  the  first  of  these  encompasses  early 
times  when  the  skin  depth  in  the  air  is  less  than  the  height  of  the  wire  over 
the  ground.  In  this  regime  we  can  calculate  the  electromagnetic  fields 
about  the  wire  in  cylindrical  coordinates  centered  on  the  wire  axis  and 
ignore  the  effects  of  the  ground.  This  regime  is  further  subdivided  into 
two  phases  depending  on  whether  the  displacement  current  is  greater  than 
tiie  conduction  current  in  the  air  or  vicc-versa.  I11  analogy  to  the  discus¬ 
sion  of  surface-burst  hMP ,  we  will  call  the  first  of  these  the  wave  phase 
and  the  second  the  early  diffusion  phase.  In  the  second  regime,  which  we 
will  call  the  late -diffusion  phase,  the  skin  depth  in  the  air  is  larger  than 
the  height  of  the  line  and  it  is  possible  to  derive  a  set  of  transmission  line 
equations  for  the  current  on  the  overhead  lino.  These  three  phases  will 
be  treated  in  separate  sections  of  this  chapter. 

6.2  WAVE  PHASE 

In  this  section  we  derive  the  current  on  the  overhead  wire  at  early 
times  when  the  displacement  current  is  much  greater  than  the  conduction 
current  in  the  air.  Neglecting  the  field  dependence  of  the  air  conductivity, 
we  can  separate  the  electromagnetic  fields  into  incident  and  scattered 
parts.  For  example  the  electric  field  parallel  to  the  wire,  oriented  for 
convenience  along  the  c  axis,  is 


14? 


total  rinc  ^scattered 
r  -•  fi  +  E 

U..7  *7  ry 


(6-3) 


The  incident  field  is  calculated  without  including  the  presence  of  the  wire 
(but  including  the  effect  of  the  air-ground  interface).  The  scattered  field 
obeys  homogenous  Maxwell's  equations  (without  the  Compton  current).  The 
effect  of  the  wire  is  incorporated  by  setting  the  sum  of  the  scattered  and 
incident  electric  fields  (the  total  field)  equal  to  zero  at  the  surface  of 
the  wire.  We  will  only  be  concerned  with  the  response  of  the  wire  to  that 
portion  of  the  incident  electric  field  parallel  to  the  wire  axis;  the 
portion  which  lies  in  the  plane  perpendicular  to  the  wire  axis  results  in 
a  polarization  of  the  wire  across  its  width  which  is  inconsequential  for 
system  survivability,  We  will  also  ignore  the  variation  in  the  parallel 
component  of  the  incident  electric  Held  across  the  wire  as  the  width  of 
the  wire  is  much  smaller  than  the  spatial  variation  of  this  field.  In  the 
wave  phase,  where 


—  >>  z  a  , 
c  0  ’ 


(6-4) 


Equations  1-17  and  1-18  for  the  scattered  field  around  the  wire  are 

1  OH 

-  -  cV  x  B  .  (6-0) 


The  component  of  the  incident  electric  field  parallel  to  the  wire  .it  a  height 
10  meters  above  the  ground  varies  at  early  times  as 


..ino  ...  f  .  u(t-zcosx/c) 

E,  =  cosxi;0(z)c  A  , 


(6-7) 


where  X  is  the  angle  between  the  radial  from  the  burst  and  the  cable,  as 

shown  in  Figure  4-1.  H^(z)  varies  slowly  as  a  function  of  distance  along 

a 

the  line  —  the  variation  results  from  attenuation  of  gammas  and  the  l/r“ 
decrease  from  a  ‘point  source.  If  we  assume  that  the  scattered  fields  have 
the  same  variation  in  z  and  t  and  ignore  the  slow  variation  of  P  in 


ts — >'•! i -■i.'.-i ■>.' i ■  i  ..CA  -l'ii'-v.  •Aif.uaV.h 


z,  Maxwell's  equations  become 


ar. 


aB  .  E  +  ,  Z  , 

0  c  r  9r 

(6-8) 

—  E^  =  acosxBg  , 

(6-9) 

a  „  c  a  .. 

—  E  =  -  ■=—  rB„  . 
c  z  r  dr  0 

(6-10) 

Eliminating  E^  and  B0  we  arrive  at  the  following  equation  for  E 


BE 


2  2 

«  sin  x  P  _  I  JL  _ _ 

2  z  r  3r  3r  * 
c 

for  which  the  solution  vanishing  at  large  r  is 

n.ii) 


Me 


(6-11) 


(6-12) 


where  Ky  is  a  modified  (hyperbolic)  Bessel  function  of  the  second  kind. 
Since  H  - 


, .  me 

I;  at  r  -  a 
z 


_  _  i„c  Hotmalnx/c) 
r,z  ”  'z  K  (aasiny/c) 


(6-13) 


As  behaves  for  large  argument  as 


K()(C) 


ir  -C 

2  F,  c 


we  associate  the  distance 
6 


c 


(6-14) 


(6-15) 


exsinx  ’ 

with  a  skin  depth  about  the  wire.  The  magnetic  field  can  be  obtained  from 
.  2  r 

asm  XB0=  +  Tr 


K.  (rotsinx/e) 
a  . .  ..me  P _ _ _ _ 

+  c  •sjnX|-z  K0(aasinx/c)  ’ 


(6-16) 
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The  current  on  the  wire  is  given  by 


important  and  currents  at  a  given  position  are  strongly  influenced  by  the 
stronger  electric  fields  nearer  the  source.  Consider  the  geometry  of 
Figure  6-1.  The  time  difference  of  arrival  between  currents  generated  at 
the  observer  location  0  and  currents  generated  at  C  and  propagating  to 
0  is 


At 


A* 

c 


(1-cosy)  , 


(6-24) 


so  that  when  the  conductivity  is  small  unu  there  is  little  attenuation  of 
propagating  fields  near  the  wire,  we  can  ignore  the  variation  of  FL  in  z 
as  long  as 


3d, 


a  ..  ..  1  () 

t  (1‘cosx)  >:>  F'-5T 


(6-J5) 


everywhere  along  the  line.  When  the  burst  is  sufficiently  close  to  the  line 
that  either  this  condition  or  the  condition  that  the  skin  depth  is  smaller 
than  the  height  of  the  line 


Figure  6-1.  Geometry  used  iri  calculation  of  limits  of 
validity  of  Equation  6-21. 
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imvMVff  mvr 


— r —  <•  h 

asiMX  ~ 


(6-26) 


is  violated,  the  situation  becomes  substantially  more  complicated.  The 
case  of  a  short  wire  high  above  the  ground  with  the  burst  located  on  the 
wire  axis  is  calculated  in  Reference  6-1;  for  a  long  wire  with  the  burst 
near  the  wire  .axis,  the  presence  of  the  finite  ground  conductivity  becomes 
important  and  this  case  has  not  been  calculated. 


6.3  EARLY  DIFFUSION  PHASE 


At  most  locations  close  to  the  burst,  the  rising  conductivity  will 
cause  the  conduct  ion  current  to  rise  above  the  displacement  curia  nt.  In 
this  section,  we  will  calculate  the  wire  current  when 

Z0o  »  ~  ,  fb-27) 

following  the  derivation  of  Reference  6-2.  The  current  on  t lie  wire  generated 
during  this  phase  should  simply  be  added  to  the  current  on  the  wire  at  the 
end  of  the  wave  phase  as  long  as  the  skin  depth  is  smaller  than  the  wave 
phase  skin  depth  (liquation  6-15).  Ignoring  the  displacement  current,  and 
assuming  that  the  fields  vary  only  as  a  function  of  retarded  time  from  the 
hurst 

t '  ~  t  -  zcosy/c  .  (6-28) 

Maxwell's  equations  become 


Dt1  Bo 

w  l ,  V  1  . 

cosy  r  ,  z 
c  3t '  Dr  * 

(6-29) 

Z0(Ilir  = 

.  . ..  315o 
cosX  t  t 

(6-50) 

Z0Oliz  " 

C  f)  .. 

f  OT  r  0  ’ 

(6-31) 
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wifi  j  *  ton  iC*i*  t  ti 


which  can  be  combined  into  a  single  equation  for  i3y 


3 

at' 


2  315  A 

c°s_£  _ P) 

1.qjc  at  '7 


_c_  a_ )  _a_ 

Z0o"  3r  f  Dr  r  0  ‘ 


(6-32) 


Dropping  the  second  term  in  parenthesis  on  tire  left-hand  side  of  (6-3?)  by 
virtue  of  (6-27),  using  the  fact  that  c/Z ^  -  1/Mq  and  expanding  the  right- 
hand  side  wo  obtain 


2 

31, 


Dr' 


+  A 

r 


313, 


v 


3B1 

at:1 


(6-33) 


Changing  to  the  sealed  variables 
R  -  r/a  , 


t ' 

t  ,  f  _i'" 

’  J 


2 

J,.u  . 

0  t 


:  t " )  ’ 


0 


where  t'  is  the  time  at  the  end  of  the  wave  phase,  wo  have 


d  B0  l  ^0  *A)  *B0 

2"  +  r  au  '  2  =  a'i-  ’ 

u  K  K 


subject  to  the  boundary  conditions. 


TO  "“o’ 


R  =  1 


(6-34) 

(6-35) 


(6-36) 


(6-37) 


lim  13  =  0  . 

R  ™  U  (6-38) 

We  first  determine  the  Croon's  function  G(R,t-t')  which  is  the  magnetic 
field  resulting  from  an  impulse  in  i;  at  T  =  t*.  Since  liquation  6-36 
possesses  translational  invariance  in  t,  we  exploit  this  fact  by  writing 
C  as  a  Fourier  transform. 
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00 


GOV 


,, .  / 


„  .  10)T  . 

G(R,w)e  dw  , 


(6-39) 


and  the  impulse  function  can  be  written 

6(t) 


1  f  iwr 

2'it  J  ° 


dia  . 


The  equation  for  G(R,w)  is 


32G 

SR2 


1  3G  /  1  ,  \  „  n 

R  3R  "  ^2  1W)G  "  0  * 


which  has  the  solution 


G(it,w)  «  G1(w)I1(  ¥lo)R)  +  c;2 (w)  K1  (  V!wR) 


(6-40) 


(641) 


(6-42) 


where  1'^  and  arc  modified  Bessel  functions  of  the  first  and  second 

kinds,  in  the  Fourier  domain,  the  boundary  conditions  on  (i  imply 


(!  (w)I  (VTw)  -  G  (w)K  (VTw)  =  -  ~~~z  > 

2’ir  V  .(a) 


(6-43) 


1  i  in  Ci  (in)  I  ,  (VEjR)  +  fi  (ai)K  ( \TTwR)  »  0  , 

U  >  u.  1  1  i  i 


(6-44) 


The  asymptotic  limits  of  the  Bessel  functions  for  largo  arguments  are 

z 


•ii 


.  ii 

ary,  z  <  - 


a  z 


~  ■  T"  c 
2  •; 


I  ^  Sn 

ary  z  <  , 


K,  (z) 

so  if  we  choose  the  argument  of  VTui  so  that 
ary  ( V  i.w)  =  -  <n  <  0  , 


ary  (  VTui)  =  +  j 


u>  >  (i 


(6-45) 

(6-46) 


(6-47) 

((-48) 
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'!■ '  1  '  ‘if!  lillilltii 


(which  will  keep  the  integrand  on  the  same  sheet  when  we  transform  back  to 
a )  Gj  must  vanish  and 

G  (oi)  =  — — - - .  (6-49) 

2ir  ViaT  Kq  (  Y±u) 

Transforming  buck  to  the  t  domain 


-.00 


(h-no) 


The  integration  can  bo  transformed  into  one  over  positive  veal  values  by 
deforming  the  contour  in  the  complex  to  plane  as  shown  in  figure  6-2.  The 
integrand  has  no  singularities  except  a  brand)  point  at  the  origin  from 
which  we  may  run  the  branch  cut  along  the  +  :lw  axis.  As  the  integrand 
vanishes  exponentially  foi  \  >  0,  U  >  1  at  large  positive  imaginary  oi, 
we  may  ignore  the  contribution  from  the  axis  G.,  and  G(. .  Near  the  origin, 
the  integrand  behaves  as 


1w 


Figure  G-2.  Integration  path  for  Green's  function. 
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inmiwi— — i 


.i.um 


?,n  (Y  >  i'o) 


so  that  wo  may  ignore  tho  contribution  from  C , ,  and  the  tivoen's  tune t i on 


i  s  now 


<i(U,T) 


l  /  ^  cawt 
■,r  J  vt: 


K  C  V~:uo  10 


r'l  ^ 

S+S 


K()(  'Hu) 


(6-52) 


:i  0 


sotting  w  «  pc'  this  can  bo  written 


Cil'R.t) 


-i  r  'ip 


-1  |  dp 

2',r  J  V7T 


0 


■^(iVTU)  KjC-i  Vp"I0 


K(){i  VT>)  K0(-iV7>) 


CO-  53) 


Impressing  the  modified  liossoJ  functions  in  terms  of  llankci  functions  of  real 
arguments,  this  is  equal  to 


«UM) 


"  2V  / 


At.  0--n> 


ujl')  (VT  in  n{:,tVTiO 


L  0 


rn[V7)) 


(VT) 


(.6-54) 


Further  simplification  is  possible  if  we  look  at  h  =  1  which  will  give  us 
t he  magnetic  field  at  the  surface  of  tho  wire.  Using  tho  Wronskian 


l[15(a)!lj2)(z)  -  ll(()1)(z)il1(2)(-/n 


4  i 
•ir  z 


((»■  55) 


wo  obtain 


G(l,x)  = 


=  _2_  fdp 

■n2  J  P  II 


"TP 


<1)(Vp)hJ2)(Vp)  ■ 


(6-56) 


Setting  x  =  V”o  end  expressing  tlio  Ilankel  functions  in  terms  of  ordinary 
liessel  functions 


-•ix 


C(1  ,x) 


4_  f  dx  _  e 

ir 2  J  X  U2(x)  + 


0>-57) 


o(xl  +  Yo(x) 
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which  is  suitable  for  machine  computation .  It  is  found  that  a  fit  to  0(1, U) 
which  is  accurate  to  within  2  percent  for  all  T  is 


(6-58) 


(('-59) 


(G- GO") 


(6-61) 


We  have  made  a  comparison  between  this  accurate  calculation  and  two 
simplified  time-dependent  inductance  models  where  we  merely  set 


!■. 
i  ne¬ 


on  , 


(6-62) 


for  tlu  first  model  and 

|t-  L(t)I(t)  -  H,ln(.Ct)  ,  (6-631 

as  the  second  model.  The  second  model  is  analogous  to  that  used  in  Chapter  4 
(liquation  4-80)  and  Chapter  F  (liquation  5-33)  in  that,  in  the  absence  of  I5., 
the  current  on  the  wire  falls  us  the  magnetic  flux  <|)  diffuses  radially  away 
from  the  wire.  (This  effect  is  absent  in  the  first  model.)  The  simple, 
approximate  forms  of  the  air  conductivity  and  the  incident  electric  field  used 
in  the  comparison  arc 


WWfJH 


.1!  !LI  Ill'll  I  "UJP.l  I 


▼>  f  i.vttt 


a  =  1  cm  , 


-  at 

F.  =  1.57  x  10  — — 
me  .  3!. 

i  +  e 


v/m 


at 


a  =  ] .57 


1  +  e 


Bt 


inho/m 


where 


8  - 1 

a  =  2  x  10  see  , 


3  =  2.4  x  io°  see-1  , 


so  that  the  peak  electric  field  is  Ur  volts/m  and  the  peak  conductivity  is 
1  mho/m.  The  inductance  used  in  the  comparison  is 


1  0 

L  «  V.n 
2ir 


/  r*  i  ,  \ 

\'V  0'Tt)a/a/ 


t  <  t 


Pk  * 


(6-64) 
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4ir 
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dt. 


M()o(t)a; 


"pk 


1,0o(tpk)u,,‘ 


t  >  t 


pk 


(6-65) 


t  >  t  ^  is  the  square  of  the  ratio  of  the  skin  depth  to  t: lie  wire  radius 


where  t  ^  is  the  time  at  which  the  conductivity  peaks.  The  term  inside  the 
logarithm  for  t  <  is  the  ratio  of  the  skin  depth  in  the  air  to  the 

wire  radius  for  an  exponentially  rising  conductivity.  Inside  the  logarithm 
for 

as  defined  by  the  substitution  (6-35)  which  gave  rise  to  the  dimensionless 
t  in  the  Green's  functions  The  second  term  in  this  logarithm  ensures 
continuity  of  l,  at  t  The  comparison  shown' in  Figure  6-3  emphasises 
the  importance  of  including  the  reduction  in  the  current  due  to  the  diffusion  of 
magnetic  flux  away  from  the  wire.  (We  have  terminated  the  curves  in  ligure 
6-3  when  the  'kin  depth  equals  10  meters.) 


i  1  i  ■■  .i’.. 
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Equation  12 


Figure  6-3.  Comparison  between  Green's  function  and 
simple  models  for  early  diffusion  phase 


L 

iJili 


t.- 1  at  •  tj  ■  atM.-T. , 


At  late  times  when  the  skin  depth  becomes  larger  than  the  wave 
phase  skin  depth  the  wave  phase  current  contribution  to  the  total  current 
should  be  decreased  according  to 


I  (t)L(t)  =  1°  L° 
wave  wave 


(6-66) 


where  T  is  the  current  at  the  end  of  the  wave  phase  and 

wave  1 

inductance  calculated  with  the  wave  phase  skin  depth. 


I  is  the 


6.4  LATE  DIFFUSION  PHASE 


In  this  section  we  will  develop  a  model  for  the  lute-time  portion  of 

the  diffusion  phase  when  the  skin  depth  in  the  air  is  greater  than  the  height 

of  the  overhead  line.  Where  the  last  two  sections  of  this  chapter  have  dealt 

with  local  phenomena  where  the  current  on  the  line  was  only  a  function  of  the 

time  histories  of  the  incident  electric  field  and  conductivity  at  that,  point, 

the  current  on  the  line  at  a  given  point  in  the  late  diffusion  phase  involves 

the  time  histories  of  the  conductivity  and  incident  field  at  other  locations 

on  the  line.  We  will  assume  a  perfectly  conducting  ground  in  this  section — 

the  effects  of  finite  ground  conductivity  on  low-frequency  signals  on  overhead 

lines  was  first  investigated  by  .J.  R.  Carson  in.  1926  and  is  reviewed  in  Sonde's 

text  (Reference  6-3).  The  features  of  his  theory  which  are  relevant  to  us  arc 

the  two  modifications  to  the  inductance  of  an  overhead  line  which  result  from 

finite  soil  conductivity — the  first  of  these  is  thac  1_ he  inductance  is  increased 

by  the  skin  depth  in  the  ground.  This  modification  is  less  than  a  factor  of 

two  change  because  the  ratio  of  the  line  radius  to  twice  the  height  is  less 

than  tlu:  ratio  of  twice  the  height  to  the  skin  depth.  (the  reason  we  use 

twice  the  height,  will  be  apparent  shortly.)  The  second  modification  is  that 

there  is  a  scries  resistance  which  results  from  the  diffusion  of  energy  into 

-4 

the  ground.  At  a  tune  of  10  seconds,  tins  scries  resistance  is  less,  than 
10  ‘  ohms/m  and  falls  as  1/t  at  later  times.  With  this  value,  the  series 
resistance  of  the  line  is  usually  much  less  than  the  termination  resi  amice 
represented  by  facilities  (as  in  Chapter  S)  if  the  burst  to  facility  distance 
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is  loss  than  a  few  kilometers.  If  it  is  not,  numerical  solutions  of  the 
transmission  line  equations  wit  the  effects  of  finite  ground  conductivity 
ma)r  be  performed. 


The  transmission  line  equations  are  derived  by  evaluating  two 
sets  of  integrals  of  Maxwell's  equations.  Applying  Stoke's  theorem  to  (1-17) 
for  the  scattered  fields,  wc  obtain 


(6-67) 


where  J)  ds  is  the  path  enclosing  the  surface 


If  we  apply  this  to  the  path  shown  in  Figure  6-1  and  designate 


V(z+6z/2)  , 


( 6 - 68 ) 


n.  (Sz 


me 


(6-69) 


Finurp  6-4.  integration  path  for  derivation  of 
first  transmission  line  equation. 


V(z-6z/2)  , 


16-70) 


ft  ■  da  ■ 


L0l6z  , 


(6-71) 


so  that.  V  is  the  voltage  between  the  line  and  ground  and  this  defines  the 
inductance  I:.d s  is  zero  as  the  scattered  field  vanishes  at  the 

surface  of  the  ground.  As  ^z  approaches  zero 


V(z+6z/2)  -  V(z-6z/2)  ^  6z  , 


(6-72) 


and  we  obtain  the  transmission  line  equation 


,  3ri_  ,  3V 

^0  3t  3z  inc 


(6-73) 


The  second  transmission  line  equation  is  derived  by  applying  Stokes'  theorem 
to  (.1-18),  neglecting  the  displacement  current 


^0  j'  E  ■  da  -  c  ^  ds  •  B  . 


Using,  the  path  shown  in  Figure  6-5,  we  have 


(6-74) 


fd'i  .  $ 


Figure  6-5.  Integration  path  for  the  derivation  of 
second  transmission  "line  equation. 
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ft  •  ds  =  -  I (z+Sz/2)  , 


(6-75) 


— -  f  B  •  ds 

*/ 


=  I(z-<5z/2)  , 


(6-76) 


a  Tn  •  da  =  6zCV  , 

J 


(6-77) 


where  we  have  used  (1-19)  in  (6-75)  and  (6-76)  and  (6-77)  defines  G,  These 
and  (6-74)  provide  the  second  transmission  line  equation 

=  -  GV  .  (6-78) 


If  the  skin  depth  in  the  air  and  the  scale  length  over  which  variations  in 
I  and  V  occur  along  z  are  much  larger  than  the  height  of  the  wire  then 
the  derivatives  with  respect  to  x  and  y  are  much  larger  than  those  with 
respect  to  r.  and  ct  and  we  may  make  use  of  two-dimensional  electrostatic 
and  magnetostatic  models  to  determine  and  G.  The  electric  field  is 

defined  as  the  gradient  of  a  scalar  potential 


0  =  -  V(,’>  , 


(6-79) 


(.when  the  air  conductivity  is  uniform  in  the  plane  perpendicular  to  the  z 
axis)  and  the  magnetic  field  as  the  curl  of  the  z  component  of  a  vector 
potential . 


B  =  V  x  A 


(6-80) 


The  requirements  tiiat  h  tangential  to  conductors  and  B  perpendicular  to 
conductors  vanish  arn  set  if  (f  and  n^e  constant  on  each  conducting 

surface.  In  the  air 


v2(p  =  v2a„  =  0  , 


(o-8 i  ) 


wher the  Laplacians  are  in  the  two  dimensions  perpendicular  to  z.  These 
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conditions  can  be  mot  for  a  wire  of  radius  a  located  a  height  h  above  a 
ground  plane  by  using  bipolar  coordinates  (see  Figure  6-6)  (Reference  6-4). 

The  transformation  to  bipolar  coordinates  is  given  by 


_ li 1  sinh£ _ 

cosh£  +  cosO  * 


(6-82) 


_ h 1 sinO 

cosh$  +  cos0  * 


(6-88) 


and  Laplace's  equation  becomes 


92<J)  a  32(()  _  (coshC+eosO) 2  /  32cj)  A  92c|>\  _  . 

i  2  12  2 ~  ^  ' 

9x“  9y  h ’  \‘3C  30  / 


(6-84) 


The  line  £  =  0  is  x  =  0,  C,  =  is  a  circle  centered  at 

x  =  h  -  h'  coth  y  -  0  , 


(6-85) 


with  radius 


a  -  h'  cscli  f, 


(6-86) 


The  solution  to  our  electrostatics  jiroblem  is 


V;  , 


a„  -  (v:  . 


(6-87) 


From  ((>-71)  and  (6-u8) 


!'0l  L2^0  ’ 


v  =  c,r. 


(6-88) 


on  f  5  S 
v0 


The  norma)  derivative  of  4> 


Figure  6-6.  Coordinate  system  used  for  calculation  of 
inductance  and  conductance  terms  in  trans¬ 
mission  line  equations. 


coshf^  +  cosO 


and  the  normal  derivative  of  A7  is  the  same  with  'C,  replaced  by  C 
GV  can  be  found  by  integrating  this  around  the  wire 


J.  I! 

C  Jt|)  oh'd( 
=  "  J  9n  eoshS.,  +' 


=  +  2noC.  , 


so  that 


cosh  ^  fh/a) 


A  similar  development  of  L  yields 


'()  .1 

L0  =  2-1 7  cosh  • 


m 


If  h/a  is  large,  we  can  use 


cosh  1  (h/a)  £n  (— )  . 

3. 

The  two  transmission  line  equations  may  be  combined  to  obtain 

T  ii  _3_  A  LA.  - 

L0  at  ‘  3z  G  3z  inc  * 


(6-93) 


(6-94) 


We  will  assume  a  simple  form  for  the  late-time  shunt  conductance,  that  it  is 
separable  in  space  and  time  and  decreases  exponentially  in  distance 

G  =  f(t)e"az  ,  (6-95) 

and  we  will  examine  an  infinite  line  running  from  -°°  to  in  z. 

To  determine  the  Green's  function  for  the  current  on  the  line,  we 
examine  the  current  which  results  from  an  incident  electric  field  of  the 
form 


Einc  =  <Hz-7.0)6(t-t0)  , 


(6-96) 


or 


r  3#  3  Otz  3 rS  r,  +  +■ 

L0t(t)  37  '  37  e  37  =  f(t)5(z-z0)6(t-t0) 


(6-97) 


This  equation  may  be  simplified  by  the  substitution 

r" 

T 

t 


f  dtj_ 

J  f(t‘) 


Pk 


(6-98) 


_3  _  d  r  9  _ 
3t  dt 

_  __1_  _D_ 
f(t)  3t 


(6-99) 


■  fctr6(T'V  ■ 


1 66 


(6-100) 


..  ii  -•  ■  i  n  i J.  •*fi~irl?>ivf'-Vri,iit/,MtLit-*  h‘  .  fciCuiW 


Equation  6-97  becomes 


.  32?  3  az  3 (S  x,  . 

L«  ---  -  —  e  .  =  6(z-z0)6(t-t0J  , 


J0  3t  3z 


0 J 


(6-101) 


which  is  invariant  under  translations  in  r,  so  that  a  frequency  domain 
technique  such  as  a  Laplace  transform  is  a  useful  technique.  The  Laplace 
transform  of  Equation  6-101  is 


/  .  3  az  3  \  . 

\sI'0  "  3z  e  3z/g  =  * 


(6-102) 


where  g(z,s)  is  the  Laplace  transform  of  77(z,t).  The  last  equation  can  be 
reduced  to  a  variation  of  Bessel's  equation  by  the  substitution 


y  r 

dv 
d  z 


-az/'< 


ay 

3 


and  it  becomes j  for  z  f  z 


0 


(6-103) 

(6-104) 


.2 


y 


3  g  3g 
'  '  y  37  ' 


4sL 


0 


y“g  =  0  . 


3y  cT 

The  solution  (in  terms  of  z)  of  this  is,  from  Reference  6-5, 


ti 


Cic 


-az/2 


i  (~  Vsi; 

1  \a 


—  -az/2' 

'I  r\  ’ 


+  C2C 


-az/2 


Kl( 


0 


l  VsIT  o-*z/2 

a  0 


)  ■ 


(6-105) 


(6-106) 


where  1^  and  are  modified  Bessel  functions.  We  now  need  to  determine 

the  values  of  Cj  and  C ^  to  use  for  z  \  z ^  to  produce  the  discontinuity 
in  the  spatial  derivative  of  g  and  satisfy  physical  boundary  conditions 
as  z  -►  ±  As  z  ->  -lo  the  argument  of  the  modified  Bessel  functions 
become  large.  If  y  is  the  argument,  the  Bessel  functions  behave  as 
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J. 


*  jyUH"  .JiL'IMHi1 '  ?  T ui.;  i*j 


V+V9  U  .{If  ^PPPW 


wu,  u  A  m  "l  ^  i 


I 


i 


V^x 


Kl" 


f  TT 

Si  2x 


-x 


((i-10'7) 


((i-lOT) 


for  large  x  s°  we  use  the  solution  that  causes  small  currents  for  z  ->■  - 
or 


g 


C^e 


-az/2 


('(i-lOD) 


for  z  <  z  .  As  z  +  <*>,  the  argument  of  the  modified  Bessel  functions 
becomes  small,  and  they  behave  as 


Ij,  -»■  |  >  (6-110) 

K1  *  X  ’  (6-1)1) 

where  x  is  again  the  argument.  At  the  present  point  in  the  derivation,  we 
have  no  a  priori  justification  for  choosing  the  combination  of  C,  and  C, 
for  ::  >  .  Ultimately,  we  will  find  that  the  clioice  is  determined  by  the 

termination  resistance  at  large  z.  for  the  present,  we  set  g  for  z  :•  z 
equal  to 


g;. 


-az/2 


-az/2\ 
c  )  , 


( d-1121 


and,  after  transforming  back  to  the  time  domain,  we  will  find  that  this  cho.ce 
is  appropriate  for  an  infinite  termination  resistance  at  large  z.  We  will 
then  determine  the  modification  necessary  for  finite  termination  resistances, 
which  will  involve  a  term  proportional  to  .  The  solution  of  liquation  6-162 
at  z  =  z()  can  be  determined  by  making  the  first  derivative  of  g  discontinuous 
at  that  point  so  that 

gj..  =  g.  I ,  ,  (6-113) 

'  “0  "0 
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*r&4l2i*j  >  Vr-.c.v  '.i .*  .r  v.'  j. 


iidiii&u... 


dz'z0  dz  'z0 


(6-1  1 4 ) 


The  derivatives  of  g  are  easily  found  with  the  use  of  the  relations 


ax  [xli(x)1  *  :<lo(x)  • 


(6- US) 


IXK^X)  I  =  -  XK0(X)  » 


(6-1161 


to  be 


,  'azr  /2  V7T*  -«z/2\ 

'1  ^1)  <:  Ma  )  ' 


(6-117) 


„  ■xr~r~  ,-azf  / 2  . —  -az/2\ 

jz  =  C2™^0  ~  koia^Lo  u  )  » 


(6-118) 


and  the  solutions  for  gs  and  g<  which  satisfy  liquation  6-102  are 


j7  _ 

sLo  Ki(x())i()(X())  4  K0(x0)iifxc) 


(6-1.10) 


Vs!,'  K  ( 


KL (x) 1 i  Cx0) 


-az/2 


CXq")  I()(Xy)  4  ^()CXqD  (Xf)) 


(6-120) 


where 


X  ~  —  V  sL~  o 

A  a  0 


(6-121 ) 


This  can  be  simplified  by  using  the  Wronskian 


k1(x)i0Cx)  +  i'n(x)i1cx)  =  1/X  . 


(6-122) 


and  the  notation 


if  z  > 


'>  ~  2  Z<  ~  =0  ’ 


if  z  < 


>  0 


so  that 


>t  r,ii.i  t  «*«>.« Mvi  itOtJOCijillkiSi  £L 


^iL  mm a m. wrt&tdiii&'SkMS&i!*  >  ■- vl  a m^ma& &£Mk&iist 


( 6-132) 


The  quantity  on  the  left-hand  side  can  be  interpreted  as  the  skin  depth  at 
the  point  zQ  at  the  time  t  and  shows  that  the  energy  is  trapped  by  the 
air  conductivity  until  the  skin  depth  becomes  greater  than  the  conductivity 
attcnuat ion  length  1/a.  The  time  of  the  peak  current  can  be  calculated 
directly  from  liquation  0-129,  by  setting  the  derivative  equal  to  zero,  and 
occurs  when 


-azn 

.  ,  0.  -az. 

Ln(e  +e  ) 


,  f  df  V" 

“  /  f(t  ')  "  2 

a 


0 

The  voltage  on  the  power  line  may  be  determined  by 

-a(z+zQ)/2 


V  - 


1_  IM 
G  3z 


L„  az 
0  e  e 

7  fCt)' 


a 


(  -az  'a“0  ) 
l  c  +  e  1 
exp/- - r - / 


r  t - =j2  - 1 - ( 

f  J*L  f  J-tL 

U  V  J  f(f) 


('6-133') 


.  ,  -az  -cC2+s(,>/2  , 

-<  W  Tj  '  c  l0l  , 


(6-1 34a ) 


where  the  arguments  of  the  Bessel  functions  arc  the  same  as  in  liquation  6-126. 
When  we  Lake  the  small  argument  limit  of  the  Bessel  functions  and  let  z  -> 
the  term  proportional  to  1  dominates  as  it  is  constant  in  z  while  the 


0 

-f;  2,  j  ' 

other  term  approaches  zero  as  e  "  . 


The  resulting  limit  for  V 


1  s 


(6- 134b) 
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provided  f(t)  is  monotonically  decreasing.  To  summarize  the  results  of 
our  investigation  so  far,  we  have  derived  the  Green's  function  appropriate 
for  an  infinite  line  with  an  infinite  termination  impedance  at  i  -  +  » 
because  the  voltage  (Equation  6-134b)  approaches  a  constant  value  as  z  -+  m 
and  the  current  approaches  zero.  We  have  discovered  the  interesting  fact 
that  the  energy  is  trapped  on  the  line  whenever  the  local  skin  depth  is  smaller 
than  the  conductivity  attenuation  length  ]/a.  When  the  skin  depth  approaches 
1/a,  the  incident  electric  field  appears  as  a  voltage  across  the  line  at  large 
values  of  z  with  a  relatively  narrow  pulse  and  a  time  integral  equal  to 

fv dt  =  J'dzJ'dt  Einc  .  (6-140') 


If  we  look  at  the  voltage  and  current  ut  a  large  but  finite  value  of  z  wc 
find  that  the  ratio  between  the  two  is 


1_ 

V 


1  -az 

-  ftc 

a 


/ 


G  ( z  '  ,  t )  d  z  ' 


(6-141) 


We  now  want  to  determine  the  modifications  needed  to  the  Green's  function 
to  represent  a  finite  termination  resistance  at  large  z.  This  may  bo 
determined  by  examining  the  response  of  Equation  6-94  to  a  step  current  at 
large  z.  Using  the  term  proportional  to  C.7  in  Equation  6-106,  the 
solution  which  results  in  a  step  function  in  time  at  t()  current,  at  z(1 
and  vanishes  as  z  -»■  -  00  at  finite  s  is 


l(z,s)  =  -■  - 


-az/2,,  ,2  -r—t—  -ccz/2. 

e  K,  C—  VsL  e  ) 


5 

e  k 


;  (- VslT 
1  \a  0 


-azo/2'  ’ 


(6-1 12) 


where  we  have  used  the  face  that  the  Laplace  transform  of  a  step  function  is 
1/s.  If  we  look  at  t  for  z  <  z^,  we  can  let  z()  -k  "  and  expand  the 
Bessel  function  in  t lie  denominator  as 
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fur  whinh  the  .  ^ce  inverse  x 


]  ( z  ,  t )  —  exp / -  — 


-az  % 


— - - \ 

j  f(f)' 


(6-145) 


Taking  the  spatial  derivative  of  T 


r, -  -az..  1 2  -j-r—  -az/2 

/  III  e  Vu^Y)e _ 

V  s  -o.z,J2  „  -a.2„ 


h(5^o 


“V2 )  ’ 


(6-146) 


if  we  now  again  let  z  ->•  +  c-°  so  we  take  the  small  argument  limit  of  K. 


J0  -az,.  1 2 


-az,.  i  t  v* —  -az/2\ 
e  k,.(— Vsn  e  ' 

O' a  0  / 


(6-147) 


The  Laplace  inverse  of  this  is 


,  -az  .  T  -az 

?J  T  L0e  L()C 

v - - - exp/ - — - 


)  2  r1  at  - 

aJ  t’Ct’T  *  a  J  f(f) 


(6-148) 


where  t  is  the  time  of  the  applied  current  and  the  voltage  resulting  from 

J  ' 7 

the  step  current  is 


V  =  - 


Lo _ 0SJ_  Hf  "xz  ) 


(6-149) 


The  phy.f .  i  ltorpretation  of  the  early  time  response  is  clarified  by 
examining  the  time  integral  of  V(z,t)  from  0  to  t.  This  integral  is 
simplified  by  the  substitution 


\  = 


L0e 


-az 


rdf 

J  flt'V 


(6-150) 


dy  - 


and  we  obtain 


L0e 


-az 


dt 


a 


dt' 

J  f(f)| 


1 2  f(t)  ’ 


to 


where 


fv dt'  =  -  —  f  -A-  e"X 

tJ  a  /  X 


'a  W  ’ 


(6-151) 


(6-152) 


L0C 


-az 


C  dt ' 

J  f(t') 


(6-155) 


and  i; ^  ( ■■■'  )  is  the  exponent  Lai  inte,  a i  .  For  small  xL)  (large  z)  the  limit 
of  the  exponential  integral  is 


!•:,  (X())  ~  "  ^(YX0)  Y  =  1 .781  , 


(6-154  i 


so 
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Vdt  ~  --  in 
a 


f  dt' 

J  f(tr 


(6-155) 


=  -  Ln(z-zJ  +  0.577  Ln/a  , 


where  z ^  is  defined  by  the  edge  of  the  source  region  where  the  skin  depth 
equals  the  exponential  conductivity  scale  length 


r _ dt  ■ 

J  V^o’1 


(6-156) 


The  late-time  response  can  also  be  understood  in  the  same  manner  —  for  late 
times  and  f  decreasing  sufficiently  rapidly  that 


t  t,-, 

[JIL.  »  f  Jll L 

J  f(f)  J  f(t!)  * 

Li  L  , 


(6-157) 


we  are  looking  into  a  resistance  with  the  value 


af  (t) 


f  dt’ 

J  Fit')' 


(6-158) 


which  is 


R  =  -  L 


0  dt  ’ 


(6-159) 


and  z  ,  is  defined  by  liquation  6-156.  To  summarize  the  outer  boundary 
conditions,  an  observer  at  z  looking  back  towards  the  source  region  sees  the 
inductance  of  the  line  between  z  and  and  a  resistance  at  the  edge 

of  the  source  region  given  by  liquation  6-158. 
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A  couple  of  useful  cases  for  which  the  time  integrals  of  ,1/f  can 
be  analytically  evaluated  are  when  f  behaves  as  an  exponential  function  or 
as  a  power  of  t.  If 


C  -  GQe 


-az-3t 


(6-160] 


f  dt'  1  ,  6t 

J  fCT'T  ’  O  (e  > 


t»t0  , 


If  we  define  z  =  0  as  the  point  where 


(6-161) 


j  UAz,t) 


(6-162) 


so  that,  at  c  =  0,  z  ~  0  Is  the  edge  of  the  source  region  then 


o '  r0F  ■ 


(6-163) 


The  Green's  function  has  a  simple  limit  when  z  <<  0,  t.  <  0,  z  >>  0,  t  >>  t 


,,,  a 

(!z  - —  e 


-a(z+z^)-26t  (  -az^-ptj 

exp,-  c  ( 


(6-164) 


which  peaks  at 


[it  =  -  azfl  -  J?.n2 


(6-165) 


The  voltage  associated  with  rff  is,  in  the  same  limit 


-az()-3t  j  -az0-htj 
V  ~  e  exp | -  e  , 


(6-166) 


. . . . 


which  peaks  at 


3t  =  -  azQ  , 


(6-167) 


which  is  the  time  that  the  edge  of  the  source  region  defined  by  Equation 
6-162  passes  z^.  The  edge  of  the  source  region  moves  in  the  -  z  direction 
at  a  uniform  velocity 


v  = 


6 

a  ’ 


(6-168) 


and  the  late-time  resistance  looking  into  the  source  region  defined  by  Equation 
6-158  is 


3L 


R  = 


0 


a 


(6-169) 


The  second  case  is  where  the  conductivity  decreases  as  a  power  of  t 


G  =  Gorue'az 


so  that 


/  -  — . —  (t 

J  f(t')  (n+l)Gn  ^ 


n+1  n+1 

"t0  J 


0 

Tn+1)(50‘  r(> 


(6-170) 


(6-171  ) 


The  edge  of  the  source  region,  defined  by 

t.  li/: 


f  Jt  _  1 

J  b()G(z,t)  ~  a 


(6-172) 


is 


i  /  2  n+i 

.  L  in  i_±± _ 

a  \  ( n+1  )  )-0f->0 


(6- 172) 
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The  limit  of  the  Green's  function  for  z  »  z^,  t  >>  t^  and  (z^t^)  deep 
within  the  source  region  is 


C"+i)\c  o 

a3t2"*2 


-a(zQ+z) 


exp 


-az. 


(n+l)G0LQe 

27n+T 
a  t 


! 

j ' 


(6-174) 


which  peaks  when 


-az 


1 


/(n+l)Lf)G0e 
\  2a2 


0  ,  n+1 


(6-175) 


The  voltage  associated  with  V  is 


-az. 


(n+1)  L0G()  -az  j  (n+l)G0L0e  ) 
- 2uiT2—  e  '  exP  \ - 2jH+T - 

at. 


I 


a  t 


which  peaks  at 


(6-176) 


t  = 


lr  ,,2,  ~aZG,  n+1 
(ii+l)  LQC>0e  \ 

(  (n+2)a2  / 


(6-177) 


The  edge  of  the  source  region  moves  in  the  -  z  direction  at  a  decreasing 
velocity 


v  = 


n  +  1 
at 


(6-178) 


and,  consequently,  the  resistance  represented  by  the  source  region  is 
decreasing 


The  results 
model  shown 
I  'm  i  n  t  z 


R 


(n+1 )  I 


at 


0 


of  this  section 
in  Figure  6-7, 
where  the  local 


(6-179) 

and  Section  6.5  can  be  summarized  in  the  eiicuit 
The  edge  of  the  source  region  is  defined  by  the 
skin  depth  equals  the  conductivity  attenuation 
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Jmi. 


length.  To  the  left  of  this  point,  the  current  on  the  line  is  proportional 
to  the  time  integral  of  the  electric  field  in  the  diffusion  phase  (plus  the 
wave  phase  current) .  The  external  circuit  shown  is  composed  of  the  voltage 
source  (Vs)  determined  from  the  Green's  function,  the  resistance  (R  )  determined 
from  the  solution  looking  back  into  the  source  region,  and  the  external  circuit 
parameters  Lp,  V  and  R  .  L  is  the  line  inductance  between  zs  and  the 
termination,  V  is  the  applied  electric  field  integrated  over  the  same 

region  and  R^  is  the  line  termination  resistance. 

The  assumption  that  the  local  current  deep  within  the  source  region 
is  determined  solely  by  the  time  integral  of  the  electric  field  is  modified 
slightly  by  the  expansion  of  the  fireball.  As  the  fireball  grows,  it 
pushes  most  of  the  magnetic  field  on  its  powerline  ahead  of  it.  The 
practical  effect  of  this  is  small  because  the  magnetic  relaxation  distance 
in  the  ionized  air  ahead  of  the  fireball  is  substantially  larger  than  the 
fireball  radius. 

6.5  APPLICATION  OF  FORMULAE  TO  EXAMPLE 

In  this  section  we  apply  the  approximate  theory  of  Section  6.4  to 
determine  the  late-time  currents  at  a  shelter.  The  parameters  of  the  example 
arc  the  same  as  that  of  Section  5.10  except  that  the  line  is  at  a  height  of 
10  meters  above  the  surface  of  the  earth.  The  small  line  resistance  (0.3 
milliohms/m)  is  ignored.  The  air  conductivity  as  a  function  of  range  at 
various  times  is  shown  in  Figure  0-8.  This  conductivity  was  calculated 
with  gamma  fluxes  of  Chapter  3  scaled  up  to  represent  a  3  MT  burst — the 
conductivity  resulting  from  device  X  rays  and  ionization  causcu  by  the 
elastic  scattering  of  neutrons  off  air  nuclei  was  ignored.  The  first  step 
in  the  calculation  was  to  determine  the  point  where 
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Air  Conduc 


Range  (m)  ► 

Figure  6-B.  Air  conductivity  as  a  function  of  range  and  fiii.e.  The  asterisks 
indicate  the  edge  of  the  source  region  and  the  vertical  bars 
indicate  the  fireball  radius. 
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t 

f  dt  =  J_ 
ju^tz.l)  "  a2  ' 
0 


(6-180) 


which  is  (6-172)  with  the  fact  that 

L0G(z,t)  =  U0a(z,t)  .  (6-18.1) 

To  do  this  we  assure  that  a  varies  locally  as 


a 


t 


-n  -az 
e 


(6-182) 


n  is  evaluated  as  a  function  of  time  from  Figure  6-8,  and  we  first  assume 
a  ~  200  meters  a(z,t)  is  then  evaluated  from 


o(z,t) 


cc2t 


(n+l)u 


0 


(6-188) 


Since  a  is  a  function  of  distance  this  formula  is  iterated  until  we  find 

the  point  where  (6-183)  is  satisfied.  One  iteration  suffices  as  a  is  a 

weak  function  of  distance,  if  we  use  the  a  from  the  last  time  step  as  the 

initial  value  in  (6-183).  The  velocity  of  the  edge  of  the  source  region  is 

evaluated  from  (6-178).  After  the  edge  of  the  source  region  reaches  the 

fireball  radius,  V  is  zero  and  R  is  the  sum  of  R  or  R_.  and  the 
s  s  s  lb 

10  ohm  termination  resistance.  The  time  integral  of  the  incident  electric 
field  is  evaluated  from  the  curves  at  the  appropriate  ranges  in  Figure 
6-9.  and 


t 


V[:  is  the  spatial  integral  of  the  electric  field  outside  the  source  rep,  ion 
anil  V.j  is  the  sum  of  and  V^.,  L  is  the  inductance  outside  the 

source  region;  we  have  included  the  effects  of  this  inductance  in  the  cal- 


(volts 


Time  (sec)  -> 

Figure  6-9.  Early  time  radial  electric 
fields  at  various  ranges. 


culation  of  the  current  in  an  iterative  manner.  An  initial  guess  at  the 
shelter  current  is  made  by 


h  =  VR 


(6-185) 


(we  have  assumed  that  the  overhead  line  extends  sufficiently  far  into  the 
fireball  to  result  in  a  low  resistance  between  the  line  and  the  fireball) 
is  calculated  and  added  to  R_  and  the  current  is  recalculated  as  I^.  As 
the  difference  between  I-j  and  1 2  is  always  less  than  about  30  percent, 
it  is  not  necessary  to  iterate  further.  These  operations  are  summarized 
in  Table  6-1  and  the  resulting  current  is  plotted  in  Figure  6-10.  The  peak 
current  occurs  somewhat  earlier  and  is  larger  than  that  shown  in  figure  5-6. 
This  results  primarily  from  the  air  conductivity  being  smaller  than  the 
ground  conductivity  at  most  ranges  so  that  the  large  electric  fields  near 
the  source  are  seen  at  an  earlier  time. 
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CHAPTER  7 


COUPLING  TO  SHORT  VERTICAL  CONDUCTORS 


7.1  INTRODUCTION 

Systems  that  are  hardened  to  source-region  EMP  from  surface  or 
near-surface  nuclear  explosions,  and  to  the  accompanying  blast  effects,  are 
unlikely  to  use  above-ground  vertical  antennas  for  mission-critical 
functions.  However,  such  antennas  and  other  vertical  conductors  may  be 
piesent  as  part  of  non-survivable  subsystems  that  serve  peacetime  functions. 
Examples  are  communications  antennas  for  maintenance  and  security  operations, 
light  poles,  etc.  It  i.s  necessary  to  know  what  currents  will  be  collected 
by  such  structures  under  EMP  conditions  and  whether  these  currents  can  get 
into  the  parts  of  the  system  that  are  required  to  survive. 

It  is  likely  that  such  structures  will  be  no  more  than  about  10 
meters  in  height  and  no  more  than  0.3  meters  in  diameter.  Thus  the  incident 
vertical  electric  field  can  be  taken  as  uniform  around  the  circumference. 
Usually,  the  conductor  will  be  thin  to  gamma  rays,  but  we  shall  include 
the  case  in  which  it  is  not.  The  conductor  may  or  may  not  be  grounded  at 
its  base,  but  since  it  is  not  likely  to  be  hardened  to  EMP,  an  arc  to  ground 
may  foim  anyway,  for  the  purposes  of  this  chapter  we  shall  assume  that  the 
base  of  the  structure  is  An  electrical  contact  with  soil;  the  impedance  of 
this  connection  is  included  in  the  analysis.  One  would  hope  that  a  wire 
does  not  run  from  the  structure  into  any  shielded  enclosure  containing  mission 
critical  electronic  equipment.  The  current  flowing  into  such  wires  could 
be  calculated,  usually  as  a  perturbation,  if  the  specifics  of  their  connection 
were  known. 
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DRIVING  FIELD  AND  SKIN  DEPTHS 


We  shall  use  the  vertical  field  Eg  :  uphei.'  _gure  3-5  as  an 

example  in  this  chapter.  That  field  is  appropriate  at  the  ground-air 
surface.  We  need  to  discuss  how  Eg  depends  on  height  above  the  surface. 

This  discussion  uses  ideas  contained  in  Sections  3.B  and  3.9. 

At  the  onset  of  a-saturation  at  the  location  of  our  conductor, 

Eg  extends  into  the  air  to  a  height  (Equation  3-49) 

<S„  c/a  ^1.5  meters  (example)  .  (7-1) 

C 

This  is  the  time  when  Efi  reaches  its  peak  value.  Previous  to  this  time 

ext  /  2 

(in  the  period  of  A-saturation)  Eg  has  been  rising  as  e  ,  and  6p 
has  been  decreasing.  The  current  in  the  conductor  at  ground  level  cannot 
be  appreciably  affected  by  the  field  at  heights  greater  than  2c/a.  Thus, 
for  times  before  the  peak  of  E0,  it  is  a  modest  overestimate  to  assume  that 
the  driving  field  is  independent  of  height  up  to  <5^. 

After  the  onset  of  a-saturation.  Eg  extends  up  to  the  conductivity- 

determined  skin  depth  5,  in  the  air.  In  this  (diffusion)  phase  also,  the 

conductor  current  at  ground  level  cannot  depend  on  the  field  at  heights 

greater  than  6  .  Thus  again  only  a  modest  overestimate  is  made  by  regarding, 
n 

Eg  as  independent  of  height  up  to  o  .  Eventually,  the  incident  Eg  is 
ine'eed  independent  o'  height,  at  least  up  to  10  meters.  We  shall  regard 
it  as  constant  in  height  at  all  times  over  that  range  of  heights  which  can 
affect  the  current  at  the  base. 


A  composite  formula  for  this  height  is 

£  =  c/a  when  a  <  c,,a  , 

a  0 


f  SBWJ5^FnfrWVTTr*»«IT,r‘ ' 


The  first  line  here  is  appropriate  up  to  the  time  of  t lie  peak  in  Eg,  the 
phase  in  which  displacement  current  exceeds  conduction  current.  The  second 
line  .is  appropriate  from  this  time  up  to  the  time  of  the  peak  in  o.  The 
third  line  is  appropriate  after  o  has  peaked,  and  is  based  on  the  assump¬ 
tion  that  a  falls  as  t  n.  The  power  n  is  determined  by 

n  -  [£n(o  /o) ] / [£n(t/t^) ]  ,  (7-3) 


where  a  -  a(t) 


and  a 

P 


aft  ).  The  value  of  n  varies  somewhat  with  t. 


It  can  be  seen  that  Equation  7-2  provides  continuous  values  of 

6  .  For  the  conductivity  of  Figure  3-3,  the  <5,  computed  from  Equation  7-2 
a  a  8  ~  i 

is  graphed  in  Figure  7-1.  For  this  case,  a  =  2  x  10‘  sec 

We  sec  from  Figure  7-1  that  reaches  10  meters  at  t  =  2  x  10 

u 

second.  As  will  be  seen,  the  coupling  to  a  conductor  of  height  h  depends 

on  whether  6  is  smaller  or  greater  than  h. 
a 


7.3  INDUCTIVELY  LIMITED  CURRENT 


The  current  in  the  conductor  is  limited,  in  varying  degrees  at  various 
times,  by  both  inductive  and  capacitive  reactance  and  by  resistance  in  both 
air  and  ground.  We  shall  first  calculate  the  current  based  on  inductive 
limitation  alone.  This  calculation  assumes  that  h  »  $  ,  that  the  ground 

tl 

is  perfectly  conducting,  and  that  the  driving  field  Eg  is  independent  of 
height.  The  current  I  in  the  conductor  is  determined  by  the  equation 

Jp  (LI)  =  Eg (t )  .  (7-4) 


The  inductance  I,  per  unit  length  is 


Mn  6„ 

in(l  +  ~ ) 


0  „  ,,  “a.  henry 
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(H  =  2  x  10 
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Figure  7-1.  Skin  depths  da  in  air  and  6 
Inductances  L(2)  and  L('IO)  of 
radius  conductors. 
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Here  a  is  the  radius  of  the  conductor.  The  argument  of  the  logarithm  has 
been  written  so  as  to  cover  the  possibility  that  6  may  be  less  than  a. 

3. 

We  shall  treat  two  examples,  with  a  =  2  cm  and  10  cm.  The  inductance  for 
these  two  cases  is  graphed  in  Figure  7-1. 

Before  a-saturation,  L  is  constant  and  E  rises  exponentially 
as  eat/^.  In  this  phase, 


I(t)  =  ~  Ea(t)  . 
v  J  aL  J 


For  the  two  conductors  chosen  as  examples, 
aL 


(7-6) 


=87  ohms/m 


=  55  ohms/m 


(7-7) 


At  the  time  of  a- saturation,  E0  reaches  its  peak  value  and  is  no  longer 
rising  exponentially.  Use  of  the  peak  value  of  Eg  from  Figure  3-5  would 
give  an  underestimate  of  the  current  at  this  time.  Use  of  the  exponentially 
extrapolated  value  of  Eg,  i.e.,  Eg  =  cB^,  will  give  an  overestimate.  We 


choose  the  latter,  setting 

E0  (peak)  =  2  x  105  V/m  . 

Thus  the  currents  at  this  time  are 

I (a-. sat.)  =  2.3  x  io^  Amp 

,3 


(7-8) 


2  cm 


(7-9) 


=  3.6  x  1Q‘  Amp  ,  a  =  10  cm  . 


.-8 


In  Figure  3-5,  a-saturation  occurs  at  1.8  x  lo  second.  Accord- 

Oft 

lag  to  Figure  3-3,  o  continues  to  rise  exponentially  (as  e  )  until  about 
-8 

5  x  lo  second,  then  changes  much  more  slowly.  The  increase  in  a  causes 


in  Equation  7-4,  as  we  shall  now  explain. 
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The  magnetic  flux  0  circling  the  conductor  is  determined  by 
integrating  the  Maxwell  equation 


3B  _  3 _E 
3t  3r 


(7-10) 


over  r  from  the  radius  a  of  the  conductor  to  large  distances.  (In  this 
equation  E  is  the  electric  field  parallel  to  the  conductor  and  B  is  the 
magnetic  field  circling  the  conductor).  Integration  over  r  yields 

It'  4 /Bdr  '  B~-EM  ■  Ee  •  (7‘n) 


where  E^  is  the  electric  field  at  large  distances  and  E(a)  =0  is  the 
field  at  the  conductor  surface  (resistance  of  conductor  neglected) .  E^  is 
the  Eg  of  the  burst  coordinate  system.  Equation  7-11  would  be  the 
same  as  Equation  7-4  if  we  define  the  inductance  by  the  Equation  $  =  LI. 

The  inductance  defined  by  Equation  7-5  assumes  that  B  ~  I/r  out  to 
r  =  a  +  6  and  then  falls  rapidly.  That  assumption  is  not  correct  when  the 
conductivity  increases  rapidly  with  time.  In  this  case,  magnetic  flux  pro¬ 
duced  in  a  given  time  interval  1/a  is  frozen  in  place  shortly  afterward 
by  the  increasing  conductivity.  The  additional  flux  produced  in  the  next, 
time  interval  is  distributed  only  over  the  decreased  skin  depth  given  by 
the  second  line  of  Equation  7-2.  Thus  the  inductance  defined  using  this  skin 
depth  applies  only  to  the  increment  in  current  in  the  next  time  interval, 
and  all  previously  established  current  is  frozen  in,  i.e.,  does  not  change 
appreciably.  Therefore,  in  the  time  period  from  a-saturation  to  the  time 

__  O 

5  x  10  1  second  when  a  stops  increasing  exponentially,  Equation  7-4  should 
be  replaced  by 


dl  _  LJ) 

dt  "  L 


(7-12) 


In  the  time  period  indicated  here,  the  ratio  Eg/I,  does  not 
change  much  for  either  of  our  example  conductors.  Average  values  are 
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l!y/L  =  1.5  x  10  1  Amp/ see  ,  a  =  2  cm  , 

=  2.9  x  10^  Amp/scc  ,  a  =  10  cm  . 

Tin.  Increase  In  current  in  this  time  period  is  therefore 
A I  =  4.8  x  1  i)'^  Amp  ,  a  =  2  cm  , 

=  9.3  x  l()'5  Amp  ,  a  =  10  cm  . 


(7  -13) 


(7-14) 


liquation  -4  would  have  given  a  larger  increase  in  current  in  this 
time  period,  since  that  equation  can  be  written  as 


1 

L  dt*  ‘ 


.7-15) 


During  the  period  after  o  peaks,  in  which  L  increases,  this  equation 
would  yield  smaller  currents  than  liquation  7-12.  Which  equation  should  we 
use  in  this  time  period?  We  argue  as  follows.  The  flux  going  with  the  current 
at  the  time  of  a-saturat ion,  liquation  7-9,  is  frozen  in  over  radii  up  to 

_7 

<5  -  1.5  meters  and  will  not  diffuse  appreciably  until  t  =  2  *  10  second 

when  A  again  reaches  that  value;  that  part  of  the  current  will  remain 
constant.  We  add  to  this  constant  current  one  half  of  A],  liquation  7-14. 

The  other  half  of  Al  is  associated  (we  say)  with  flux  distributed  only  up 

to  (S  (mi.it)  =  0.09  meters,  and  this  flux  diffuses  immediately,  together  with 

'  -8 
additional,  flux  produced  by  li^  after  t  =  5  x  10  second,  for  this 

part  of  the  current  liquation  7-4  Is  appropriate.  Thus  the  current  has 

two  parts,  for  which,  until  t  =  1  x  10  1  second, 


1|  =  constant  =  4.7  x  10‘  Am])  ,  a  -  2  cm  , 

=  8.2  x  nr’’  Amp  ,  a  =  10  cm  , 


(7-lh) 


dT  ai2>  =  ,:o 


(7-17) 


where ,  at  t  =  5  y  10  second. 


>•  vi-  a --i'1.', ..uAiiJ;..., .  c. f,** i*y I’-iir  f>lc. 


9 


2  cm 


(7-18) 


U  = 


=  2.4  x  K)‘  Amp  ,  a  =  2  cm  ,  j 

(7-18) 

-  4.6  x  io"  Amp  ,  a  =  10  cm  .  ) 

Hand  integration  of  liquation  7-17  in  three  steps  yields  for  the  total  current 


5  x  io"8 

CO 

1 

i — < 

X 

r>. 

1  x  IQ'7 

2  x  io"7  sec 

) 

7.1 

8.5 

8.4 

1 

9.6  KA,  a  =  2  cm  /  (7-19) 

12.8 

15.9 

14.5 

15.5  KA,  a  =10  cm  > 

I 

The  currents  are  graphed  as  a  function  of  time  in  Figure  7-2. 


After  t  -  2  x  io  7  second,  all  of  the  flux  diffuses,  and  Equation 
7-2  is  appropriate.  Since  Eq  is  approximately  constant  in  this  period, 
that  equation  yields 

LI  =  (LI).  +  (t-t.)H0  .  (7-20) 

where  the  subscript  i  indicates  evaluation  at  t  =  2  x  10 second. 
Currents  calculated  from  this  equation  are  used  to  extend  the  curves  in 
Figure  7-2  out  to  10  5  second. 


7.4  FFFFCT  OF  GROUND  TERMINATION 


In  considering  the  effect  of  a  finitely  conducting  ground  termina¬ 
tion  of  the  vertical  conductor,  it  is  important  to  understand  that  the 
current  carried  by  the  conductor  is  exactly  equal  to  the  current  (conduction 
and  displacement)  removed  from  the  surrounding  medium  due  to  the  presence 
of  the  conductor.  This  follows  from  the  fact  that  the  magnetic  field  at 
distances  appreciably  larger  than  t he  skin  depth  is  unaffected  by  the 
presence  of  the  conductor;  hence  the  net  change  in  current  over  an  area 
comparable  with  iuS~  must  vanish  by  Stokes'  theorem. 

The  flow  of  (tiie  change  in)  current  in  the  ground  must  therefore 
be  as  sketched  in  Figures  7-5a  and  7 -3b  for  the  cases  in  which  the  air 


Current  (amperes) 


Time  (seconds) 


Figure  7-2.  Currents  at  base  of  conductor.  I[_,  inductively  limited 

current;  Ir,  resistively  limited  current;  IR0,  resistively 
limited  current  with  perfectly  conducting  ground;  Is, 
Compton  current  collected  by  aluminum  pipe  conductor. 
Number  in  parentheses  is  radius  of  conductor  in  cm. 
Conductor  (2)  is  3  m  high,  conductor  (10)  is  10  m  high. 
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In  the  time  period  before  u-saturat ion ,  the  conductor  current  rises 

.  , ,  at/2  T  .  .  .  . 

exponentially  as  e  .  In  this  period, 


C/U 


’  1 1>  \/  i.i  .o 


17-2  1  ) 


O  L> 


SI  s 

With  a  -  2  r  1()‘  sec"  and  a  =  1.6  *  10  '  miio/m  (10  percent  water  soil 


at  frequency  l.o  *  !()'  lie,  see  Chapter  2)  we  obtain 


6  =  1.5  m  ,  6  -  0.7  m  , 

a  8 


(7-25) 


liquation  7-21  qives 

=  01  ohms  ,  a  -  2  cm  ,  | 

=  30  ohms  j  a  =  10  cm  .  ^ 


("-26) 


i  A 


Now  according  to  liquation  7-7,  the  inductive  reactance  ot  a  length  *5  =  1.5 
meters  of  the  conductor  Is  130  and  82  ohms  in  the  two  cases.  Thus  the 
current  before  u-saturat i on  will  be  reduced  bv  a  factor  of  about 


1 30  o  ■  <82  _ 

130"+  61  •  ~  82  +  3 ~9  ’ 


(7-27) 


in  both  cases  because  of  the  termination  resistance. 

-  8 

By  t  -  5  <  10  second,  the  rise  time  i  ( e- 1  old inq  time)  of  the 
conductor  current  has  increased  to  about 


T  !-  2.6  x  10  >S  second  . 


(7-28) 


At.  frequency  !/2ai  -  (>  *  10  lie,  Chapter  2  qives  •’i  =  1 . 3  x  10  "  mhu/m, 

and 


fce.,  1, 


"i;  "  yJ\T>i 

(i  r 


1  ,  3  ill  . 


(7-2(1) 


As  described  in  Section  7.3,  the  magnetic  flux  in  the  air  is  distributed  over 
a  ranj’.e  of  radii  (^'s)  from  0.0P  to  1.5  meters.  The  geometric  mean  yives 
A  -  0.37  rioter.  liquation  7-22  then  pives  the  termination  resistance 
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R  -  39  ohms  ,  a  =  2  cm 


(7-30) 


=  24  ohms  ,  a  =  10  cm  . 

At  6  =0.37  meter,  the  inductance  of  a  length  6  of  the  conductor  is 

a  h  a 

LS_  =  0.22  Jill  ,  a  :  2  cm  , 

=  0.11  pi  1  ,  a  =  10  cm  . 

The  relaxation  time  of  the  current  into  the  ground  is 

L6u  R 

=  0.56  x  io  sec  ,  a  =  2  cm  , 

t  (7-321 

=  0.46  x  10  S  sec  ,  a  =  10  cm  . 

because  this  relaxation  time  is  short  compared  with  the  rise  time  t  of 
the  conductor  current,  most  of  the  conductor  current  indicated  in  Tigure  7- 
wi.ll  flow  in  the  air  just  above  the  ground  rather  than  in  the  ground.  The 
fraction  of  the  current  flowing  in  the  ground  is  about 


(7-31) 


Vg  _  (0.013)  (1.3) 

a  6  +  a  6  "  0.0169  +  (0.4)  (0.37) 

g  g  a  a 


0.10  . 


f  7 -33) 


Thus  the  ground  termination  resistance  substantially  reduces  the  current  in 

_ 

the  base  of  the  conductor  at  times  in  the  neighborhood  of  5  x  10  second. 
If  the  ground  surface  is  covered  by  a  conducting  slice t  or  counter  poise, 
the  current  in  Figure  7-2  is  correct. 


7.5  RES ISTIVELY  LIMITED  CURRENT 


Section  7.4  has  shown  that  resistance  of  the  ground  termination 
affects  the  current  at  the  base  of  the  conductor  at  quite  early  times. 

As  the  air  conductivity  falls  (Figure  3-3),  resistance  in  the  air  will  also 
limit  tiie  current.  In  this  section,  we  shall  ignore  inductive  effects  and 
calculate  the  current  as  in  a  static  problem.  It  is  helpful  to  distinguish 


between  two  phases,  according  to  whether  the  skin  depth  6  in  the  air  is 

cl 

less  or  larger  than  the  height  h  of  the  conductor. 


First  Phase:  6  <  h. 

a 


The  vertical  electric  field  Eg  is  not  just  a  static  field;  the 
curl  of  E  does  not  vanish  at  early  times  since,  as  shown  by  Figure  3-5, 

_7 

B.  changes  rapidly  with  time  before  t  =  10  second.  However,  that  does 
not  matter  for  the  conductor  response,  which  depends  only  on  the  vertical 
electric  field.  The  same  Eg  at  the  position  of  the  (thin)  conductor,  even 
if  it  were  derivable  from  a  potential,  would  produce  the  same  current  in 
the  conductor.  We  can  therefore,  for  the  convenience  of  familiarity,  think 
in  terms  of  a  voltage  V(z), 


V(z)  =  j HgU'jdt'  , 


(7-34) 


where  z  is  the  distance  above  the  ground  surface. 

In  the  diffusion  phase  (which  begins  at  a-saturat ion) ,  Eg  extends 
only  up  to  the  skin  depth  6.  above  the  ground;  above  that  height  the 
electric  field  is  approximately  radial  from  the  burst  point.  Thus  V(z) 
has  the  z-dependence  indicated  in  Figure  7-4a.  The  maximum  voltage  V 
i  s  about 


V  «  Un6  , 
m  0  a 


(7-35) 


where  Ff1  is  the  field  iust  above  the  ground. 


If  we  imagine  the  conductor  to  be  opened  just  above  the  ground,  as 
in  Figure  7-4b,  to  what  voltage  will  the  conductor  come?  Remember  that 
the  conductor  is  in  a  conducting  medium.  Because  currents  in  the  conductor 
are  limited  by  diffusion  in  this  medium,  the  lower  end  of  the  opened 
conductor  cannot  be  affected  by  conditions  existing  at  heights  much  greater 


Figure  7-4.  Potential  function  and  open  circuited 
conductor. 

than  6,  .  We  conclude  therefore  that  the  open-circuit  voltage  of  the 

tl 

conductor  is  about 

Vn  »  V  /2  *  Li,,5  /2  .  (7-36) 

0  m  0  a  v 

If  we  now  reconnect  the  conductor  to  its  ground  end,  a  current 
will  flow  across  the  junction.  The  magnitude  of  thi  a  current  will  he 

I  -  V0/(R.+Rt)  ,  (7-37) 

where  R.  is  the  "internal"  resistance  of  the  source  of  V,.  and  R.  is 
i  0  t 

ground  termination  resistar.ee.  R  has  been  estimated  by  liquations  7-21, 
and  23.  An  estimate  of  IT  is 

b  -  shr 1,1(1  *  t’  •  (7'38) 

a  a 

This  is  the  resistance  between  a  conductor  of  radius  a,  length  <5^  and 
distant  points  in  a  medium  of  conductivity  ojf 

.ce:ond  Phase:  6  >  h. 

a 


In  this  phase  the  open-circuit  voltage  is  about 


(7-39) 


V  w  [•'  U/~> 

V()  LyM/-  , 

and  the  reconnected  current  is 


1  =  V0/fR*+Rp  . 


(7-40) 


The  internal  resistance  here  is  modified  by  replacing  6  with  h ,  i.c., 


R!  =  -  £n  (--)  , 

l  2  it  a  h  hi  * 

a 


(7-41) 


anil  the  termination  resistance  is 


R*  =  ;--1— r-  £n(h-)  ,  (h  >  d  >  6  ) 

t  duo  A  hi  ’  1  g 

S  l! 


(7-42) 


=  r  | fcn (— )  +  1  -  — |  +  fi-n (4^—)  ,  (h  >  A  >  d)  , 

2ira  d  1  hr  A  1  2iro  A  '■O  ’  g 

R  R  R  R  S 

(7-43) 


=  [  ftn  -  +  1  -  '’  I  ,  (iS  >  h  >  d) 

2'iro  d  1  a  h 1  ’  g 

g 


(7-44) 


The  estimates  given  here  are  continuous  between  the  three  regimes  of  A  . 

bet  us  calculate  the  resistively  limited  current  for  our  two 
examples,  choosing  the  heights  and  depths 


it  -  3  m  ,  d  -  J  in  for  a  =  2  cm  , 

=  10  ill  ,  d  =  1  m  for  a  =  10  cm  . 


(7-43) 


The  skin  depth  6  in  the  air  has  been  graphed  in  figures  7-1.  for  the  soil 

we  take  o  =  0.01  mho/meter  and 

g 


/  2  t 1 

1(0.7  m)  +  -  ----  , 

0  g 


(7-40) 


where 


f  =  t  -  2  x  10  sec  . 


(7-47) 


This  estimate  of  <5^  is  continuous  with  the  value  in  the  exponential  phase 
given  by  liquation  7-25.  6  is  also  graphed  in  Figure  7-1. 

it 

The  first  phase,  6^  <  h,  ends  at 


t^  =  3.7  x  10  sec  for  a  =  2  cm  , 


=  2.0  x  io  sec  for  a  =  10  cm 


(7-48) 


In  this  phase  V  is  given  by  Equation  7-36  and  lb  by  Equation  7-38.  Thus 
quantities  are  graphed  in  Figure  7-5,  where  they  are  also  extended  into  the 
second  phase,  6  >  h,  by  use  of  Equations  7-39  and  7-41. 

For  the  termination  resistance,  we  sec  that  5^  is  greater  than 

the  assumed  d  at  almost  all  times  of  interest,  and  that  6  >6  in  most 

g  a 

of  the  first  phase  (when  d  <  h) .  Since  liquations  7-21,  22,  23  do  not  apply 

cl 

in  tills  case,  we  need  another  estimate  of  R.  for  the  case  h  >  <S  ,  iS  >  6  , 
'  t  a  g  a  ’ 

3  >  d.  This  is 


*t  ■  SF 51  *»<5>  ■  <5,  <  J1  ' 

& 


(7-4  9) 


L_  [ g.n (---)  +  1  -  ,  (6  >  d)  . 

2'iro  d  1  \n  iS  1  ’  a 


(7-50) 


From  Figure  7-1  we  see  that  cS  exceeds  d  after  t  =  1.6  x  K)  second. 

Before  this  time,  Equation  7-19  gives  R  .  From  this  time  until  tj 
(liquation  7-48),  R  is  given  by  Equation  7-50.  After  t ^ ,  R^  is  given  by 
Equation  7-44.  R  is  also  graphed  in  Figure  7-5  for  the  two  conductors. 


In  computing  R.  at  times  before  u-saturat ion  in  Figure  7-5,  we 
have  replaced  o  by  u  yx  in  liquations  7-38  and  7-41,  since  the  displace¬ 
ment  current  is  larger  than  the  conduction  current  in  that  phase.  Hence 
the  current  is  limited  by  capacitive  reactance  rather  than  by  air  resistance 
in  that  phase.  Note  that,  the  capacitive  reactance  is  real  For  exponentially 
rising,  field,  and  is  approximately  constant. 


Time  (second) 

Figure  7-5.  Open  circuit  voltage  V q  of  two  vertical  conductors: 

(2)  2  cm  radius,  3  in  high;  (10)  10  cm  radius,  10  m  high. 
Source  resistance  Ri  and  termination  resistance 
for  the  two  conductors. 
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The  resistively  (and  capaeitvely)  limited  currents  calculated 
from  Equations  7-37  and  7-40  are  graphed  in  Figure  7-2,  for  the  two  conductor 
radii  and  for  both  perfectly  and  imperfectly  conducting  ground.  It  is  seen 
that  these  currents  are  substantially  smaller  at  all  times  than  the 
inductively  limited  current.  For  the  examples  considered,  inductive 
reactance  is  small  at  all  times  compared  with  resistance  or  capacitive 
reactance.  The  correct  current  T  can  be  estimated  from  the  equation 

1  *  Vl^V'V  *  (7-51) 

where  1^  and  Ij,  are  the  resistively  and  inductively  limited  currents. 

I  is  only  a  little  less  than  1^  in  our  examples.  The  sign  of  the  current 
is  such  that  electrons  flow  down  the  conductor  into  the  ground. 

7.6  COMPTON  CURRENT  COLLECTION 

If  the  conductor  provides  appreciable  attenuation  of  t lie  gamma  rays, 
the  Compton  current  emerging  from  its  back  side  will  be  less  than  that  entering 
its  front  side.  Thus  the  conductor  collects  negative  charge  due  to  gamma 
attenuation,  (lamina  attenuation  lengths  are  of  the  order  of  30  grams/em“. 

If  tiie  conductor  is  made  of  high-utomic-number  material,  then  the 
ratio  of  Compton  electron  flux  to  gamma  flux  coming  out  of  the  conductor  is 

smaller  than  the  ratio  going  in,  provided  the  conductor  is  an  electron  range 

.  .  4  2 
m  thickness.  Compton  electron  ranges  arc  of  the  order  of  0.3  grams/cm 

in  air  and  aluminum,  but  arc  smaller  (due  to  nuclear  scattering)  in  high  Z 

materials.  Most  conductors  will  be  thicker  than  an  electron  range.  An  iron 

conductor  collects  about  30  percent  of  the  Compton  current  striking  it,  even 

without  the  gamma  attenuation  effect. 

Let  us  assume  that  our  10  cm  radius  conductor  is  an  aluminum  pipe 

3 

(density  p  =  2,7  gram/cm')  with  wall  thickness  D  -  1/4  inch  0.0  cm.  Its 
average  projected  thickness  in  grams/cm^  is  then 


m 


=  'irDp  =  5.1  gram/cm 


(7-52) 


t-  M  ci  U  (J 
~2 ~ 


The  fraction  of  Compton  current  collected,  due  to  gamma  attenuation,  is 

f  =  1  ~  exp (-m/ 30)  ^  0.16  .  (7-53) 

The  total  Compton  current  collected  in  a  height  of  the  pipe  equal  to  the 
smaller  (6r  ,h)  of  6  and  h  is 

rl  ll 

1  -  2  (1  a (6  , h)  ,  (7-54) 

S  b  il 


where  is  the  Compton  current  density,  which  we  take  from  figure  5-2. 

The  current  computed  from  liquation  7-54  is  also  graphed  in  Figure  7-2.  Note 
that  it  is  small  compared  with  the  other  currents,  except  for  a  short  period  of 
time  for  the  resistvely  Limited  currents  into  imperfectly  conducting  ground. 
Actually,  the  current  I  affects  the  resistvely  limited  current.  Since 
the  impedance  of  the  Compton  current  source  is  very  large,  the  current  into 
the  ground  is  given  by 


I  = 


I  R. 

R  i  ♦  Rt"  +  R;p-R; 


R. 


s  R.  +  R 
i  t 


(7-55) 


where  1^  is  the  resist ively  limited  current  calculated  in  Section  7.5  Iv’e 
see  from  Figure  7-5  that  the  ratio  R./(R,+R  )  is  small  (T°  0.1)  in  the 
time  period  in  which  I  exceeds  i^.  Thus  the  collected  Compton  current 
makes  little  difference  for  our  aluminum  conductor.  If  the  pipe  were  made 
of  iron,  f  and  I  would  lie  about  four  times  larger  than  for  aluminum. 


7.7  NUCLEAR  LIGHTNING 


It  is  likely  that  a  discharge  would  form  in  the  air  at  the  upper  end 
of  the  conductors  in  our  examples,  and  grow  upwards.  Such  discharges  were 


observed  in  several  large  yj <  Id  nuclear  tests.  The  theory  of  this  "nuclear 
lightning"  is  currently  under  development  at  MRC  (by  J.  Gilbert,  It.  Gardner, 

M.  Frese,  and  C.  Longmirc) .  It  is  believed  that  the  currents  in  these 
discharges  reach  peak  values  of  several  times  1(/*  amperes,  much  larger 
than  the  resistively  limited  currents  of  Figure  7-2  because  of  the  increased 
height  of  the  discharges.  These  heights  were  observed  to  reach  several  hundred 
meters  in  the  millisecond  time  frame.  The  authors  hope  to  add  a  chapter  on 
nuclear  lightning  to  this  report  when  the  theory  is  Firmly  established. 
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